double lat(lat) ;
	double lon(lon) ;
	double nlev(nlev) ;
	double nlevp1(nlevp1) ;
	double time(time) ;
	double height_f(nlev, lat, lon) ;
	double pressure_f(nlev, lat, lon) ;
	double temp(nlev, lat, lon) ;
	double qv(nlev, lat, lon) ;
	double rh(time, nlev, lat, lon) ;
	double theta(time, nlev, lat, lon) ;
	double rv(time, nlev, lat, lon) ;
	double u(nlev, lat, lon) ;
	double v(nlev, lat, lon) ;
	double ug(time, nlev, lat, lon) ;
	double vg(time, nlev, lat, lon) ;
	double w(time, nlev, lat, lon) ;
	double uadv(time, nlev, lat, lon) ;
	double uadvh(time, nlev, lat, lon) ;
	double uadvv(time, nlev, lat, lon) ;
	double vadv(time, nlev, lat, lon) ;
	double vadvh(time, nlev, lat, lon) ;
	double vadvv(time, nlev, lat, lon) ;
	double tadv(time, nlev, lat, lon) ;
	double tadvh(time, nlev, lat, lon) ;
	double tadvv(time, nlev, lat, lon) ;
	double qadv(time, nlev, lat, lon) ;
	double qadvh(time, nlev, lat, lon) ;
	double qadvv(time, nlev, lat, lon) ;
	double thadv(time, nlev, lat, lon) ;
	double thadvh(time, nlev, lat, lon) ;
	double thadvv(time, nlev, lat, lon) ;
	double radv(time, nlev, lat, lon) ;
	double radvh(time, nlev, lat, lon) ;
	double radvv(time, nlev, lat, lon) ;
	double radcool(time, nlev, lat, lon) ;
	double sfc_sens_flx(time, lat, lon) ;
	double sfc_lat_flx(time, lat, lon) ;
	double ts(time, lat, lon) ;
	double ustar(time, lat, lon) ;
	double q1(time, nlev, lat, lon) ;
	double q2(time, nlev, lat, lon) ;
	double ustress(time, nlev, lat, lon) ;
	double vstress(time, nlev, lat, lon) ;
	double ps(time, lat, lon) ;
	double omega(time, nlev, lat, lon) ;
	double orog(lat, lon) ;
	double height_h(nlevp1, lat, lon) ;
	double pressure_h(nlevp1, lat, lon) ;
	double coor_par_a(nlevp1, lat, lon) ;
	double coor_par_b(nlevp1, lat, lon) ;