Changeset 2787 for LMDZ5/branches/testing/libf/phylmd/Dust
- Timestamp:
- Jan 30, 2017, 5:54:45 PM (8 years ago)
- Location:
- LMDZ5/branches/testing
- Files:
-
- 4 edited
- 4 copied
-
. (modified) (1 prop)
-
libf/phylmd/Dust/phys_output_write_spl_mod.F90 (modified) (44 diffs)
-
libf/phylmd/Dust/spla_output_dat.h (modified) (1 diff)
-
libf/phylmd/Dust/spla_output_write.h (modified) (7 diffs)
-
libf/phylmd/Dust/splaeropt_5wv_rrtm.F90 (copied) (copied from LMDZ5/trunk/libf/phylmd/Dust/splaeropt_5wv_rrtm.F90)
-
libf/phylmd/Dust/splaeropt_6bands_rrtm.F90 (copied) (copied from LMDZ5/trunk/libf/phylmd/Dust/splaeropt_6bands_rrtm.F90)
-
libf/phylmd/Dust/splaeropt_lw_rrtm.F90 (copied) (copied from LMDZ5/trunk/libf/phylmd/Dust/splaeropt_lw_rrtm.F90)
-
libf/phylmd/Dust/splaerosol_optic_rrtm.F90 (copied) (copied from LMDZ5/trunk/libf/phylmd/Dust/splaerosol_optic_rrtm.F90)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/branches/testing
- Property svn:mergeinfo changed
/LMDZ5/trunk merged: 2728,2730-2775,2777-2785
- Property svn:mergeinfo changed
-
LMDZ5/branches/testing/libf/phylmd/Dust/phys_output_write_spl_mod.F90
r2641 r2787 6 6 !JE20150620<< 7 7 !JE20150620>> 8 9 8 !JE20150620<< 10 9 11 USE time_phylmdz_mod, only: day_step_phy, start_time, itau_phy10 USE time_phylmdz_mod, ONLY: day_step_phy, start_time, itau_phy 12 11 13 12 USE phytracr_spl_mod, ONLY : ok_chimeredust, id_prec, id_fine, id_coss, & … … 93 92 ! defined and initialised in phys_output_mod.F90 94 93 95 USE dimphy, only: klon, klev, klevp1, nslay 96 USE control_mod, only: day_step, iphysiq 97 USE phys_output_ctrlout_mod, only: o_phis, o_aire, is_ter, is_lic, is_oce, & 94 USE dimphy, ONLY: klon, klev, klevp1 95 USE ocean_slab_mod, ONLY: nslay 96 USE control_mod, ONLY: day_step, iphysiq 97 USE phys_output_ctrlout_mod, ONLY: o_phis, o_aire, is_ter, is_lic, is_oce, & 98 98 is_ave, is_sic, o_contfracATM, o_contfracOR, & 99 99 o_aireTER, o_flat, o_slp, o_tsol, & … … 287 287 !JE20150620>> 288 288 289 USE phys_state_var_mod, only: pctsrf, paire_ter, rain_fall, snow_fall, &289 USE phys_state_var_mod, ONLY: pctsrf, paire_ter, rain_fall, snow_fall, & 290 290 qsol, z0m, z0h, fevap, agesno, & 291 291 nday_rain, rain_con, snow_con, & … … 314 314 vphiSTD, wTSTD, u2STD, v2STD, T2STD, missing_val_nf90 315 315 316 USE phys_local_var_mod, only: zxfluxlat, slp, zxtsol, zt2m, &316 USE phys_local_var_mod, ONLY: zxfluxlat, slp, zxtsol, zt2m, & 317 317 t2m_min_mon, t2m_max_mon, evap, & 318 318 zu10m, zv10m, zq2m, zustar, zxqsurf, & … … 369 369 d_q_ch4, pmfd, pmfu, ref_liq, ref_ice, rhwriteSTD 370 370 371 USE phys_output_var_mod, only: vars_defined, snow_o, zfra_o, bils_diss, &371 USE phys_output_var_mod, ONLY: vars_defined, snow_o, zfra_o, bils_diss, & 372 372 bils_ec,bils_ech, bils_tke, bils_kinetic, bils_latent, bils_enthalp, & 373 373 itau_con, nfiles, clef_files, nid_files, zvstr_gwd_rando 374 USE ocean_slab_mod, only: tslab, slab_bils, slab_bilg, tice, seaice375 USE pbl_surface_mod, only: snow376 USE indice_sol_mod, only: nbsrf377 USE infotrac, only: nqtot, nqo, nbtr, type_trac378 USE geometry_mod, only: cell_area379 USE surface_data, only: type_ocean, version_ocean, ok_veget, ok_snow380 ! USE aero_mod, only: naero_spc381 USE aero_mod, only: naero_tot, id_STRAT_phy382 USE ioipsl, only: histend, histsync383 USE iophy, only: set_itau_iophy, histwrite_phy384 USE netcdf, only: nf90_fill_real374 USE ocean_slab_mod, ONLY: tslab, slab_bils, slab_bilg, tice, seaice 375 USE pbl_surface_mod, ONLY: snow 376 USE indice_sol_mod, ONLY: nbsrf 377 USE infotrac, ONLY: nqtot, nqo, nbtr, type_trac 378 USE geometry_mod, ONLY: cell_area 379 USE surface_data, ONLY: type_ocean, version_ocean, ok_veget, ok_snow 380 ! USE aero_mod, ONLY: naero_spc 381 USE aero_mod, ONLY: naero_tot, id_STRAT_phy 382 USE ioipsl, ONLY: histend, histsync 383 USE iophy, ONLY: set_itau_iophy, histwrite_phy 384 USE netcdf, ONLY: nf90_fill_real 385 385 386 386 #ifdef CPP_XIOS 387 387 ! ug Pour les sorties XIOS 388 388 USE xios, ONLY: xios_update_calendar 389 USE wxios, only: wxios_closedef, missing_val389 USE wxios, ONLY: wxios_closedef, missing_val 390 390 #endif 391 USE phys_cal_mod, only : mth_len 392 391 USE phys_cal_mod, ONLY : mth_len 393 392 394 393 IMPLICIT NONE 395 396 394 397 395 ! INCLUDE "temps.h" … … 445 443 CALL set_itau_iophy(itau_w) 446 444 447 IF (.NOT.vars_defined) THEN445 IF (.NOT.vars_defined) THEN 448 446 iinitend = 2 449 447 ELSE … … 456 454 !$OMP MASTER 457 455 IF (vars_defined) THEN 458 if (prt_level >= 10) then456 IF (prt_level >= 10) THEN 459 457 write(lunout,*)"phys_output_write: call xios_update_calendar, itau_w=",itau_w 460 endif458 ENDIF 461 459 ! CALL xios_update_calendar(itau_w) 462 460 CALL xios_update_calendar(itap) 463 END 461 ENDIF 464 462 !$OMP END MASTER 465 463 !$OMP BARRIER … … 751 749 CALL histwrite_phy(o_uq, uq) 752 750 CALL histwrite_phy(o_vq, vq) 753 IF (iflag_con.GE.3) THEN ! sb751 IF (iflag_con.GE.3) THEN ! sb 754 752 CALL histwrite_phy(o_cape, cape) 755 753 CALL histwrite_phy(o_pbase, ema_pcb) 756 754 CALL histwrite_phy(o_ptop, ema_pct) 757 755 CALL histwrite_phy(o_fbase, ema_cbmf) 758 if (iflag_con /= 30) then756 IF (iflag_con /= 30) THEN 759 757 CALL histwrite_phy(o_plcl, plcl) 760 758 CALL histwrite_phy(o_plfc, plfc) 761 759 CALL histwrite_phy(o_wbeff, wbeff) 762 end if760 ENDIF 763 761 764 762 CALL histwrite_phy(o_cape_max, cape) … … 771 769 CALL histwrite_phy(o_ftime_con, zx_tmp_fi2d) 772 770 IF (vars_defined) THEN 773 IF (iflag_thermals>=1)THEN771 IF (iflag_thermals>=1)THEN 774 772 zx_tmp_fi3d=dnwd+dnwd0+upwd+fm_therm(:,1:klev) 775 773 ELSE … … 820 818 DO k=1, nlevSTD 821 819 bb2=clevSTD(k) 822 IF (bb2.EQ."850".OR.bb2.EQ."700".OR. &820 IF (bb2.EQ."850".OR.bb2.EQ."700".OR. & 823 821 bb2.EQ."500".OR.bb2.EQ."200".OR. & 824 822 bb2.EQ."100".OR. & … … 840 838 #endif 841 839 #endif 840 842 841 #ifdef CPP_XIOS 843 IF (ok_all_xml) THEN842 IF (ok_all_xml) THEN 844 843 !XIOS CALL xios_get_field_attr("u850",default_value=missing_val) 845 844 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! … … 847 846 DO k=1, nlevSTD 848 847 bb2=clevSTD(k) 849 IF (bb2.EQ."850".OR.bb2.EQ."700".OR. &848 IF (bb2.EQ."850".OR.bb2.EQ."700".OR. & 850 849 bb2.EQ."500".OR.bb2.EQ."200".OR. & 851 850 bb2.EQ."100".OR. & … … 940 939 ELSE 941 940 CALL histwrite_phy(o_tslab, tslab) 942 END 941 ENDIF 943 942 IF (version_ocean=='sicINT') THEN 944 943 CALL histwrite_phy(o_slab_bilg, slab_bilg) 945 944 CALL histwrite_phy(o_slab_tice, tice) 946 945 CALL histwrite_phy(o_slab_sic, seaice) 947 END 946 ENDIF 948 947 ENDIF !type_ocean == force/slab 949 948 CALL histwrite_phy(o_weakinv, weak_inversion) … … 1029 1028 CALL histwrite_phy(o_tausumaero(naero), & 1030 1029 tausum_aero(:,2,naero) ) 1031 END 1030 ENDDO 1032 1031 ENDIF 1033 1032 IF (flag_aerosol_strat>=1) THEN … … 1046 1045 CALL histwrite_phy(o_sollwad0, sollwad0_aero) 1047 1046 !====MS forcing diagnostics 1048 if (new_aod) then1047 IF (new_aod) THEN 1049 1048 CALL histwrite_phy(o_swtoaas_nat, topsw_aero(:,1)) 1050 1049 CALL histwrite_phy(o_swsrfas_nat, solsw_aero(:,1)) … … 1057 1056 CALL histwrite_phy(o_swsrfcs_ant, solsw0_aero(:,2)) 1058 1057 !cf 1059 if (.not. aerosol_couple) then1058 IF (.not. aerosol_couple) THEN 1060 1059 CALL histwrite_phy(o_swtoacf_nat, topswcf_aero(:,1)) 1061 1060 CALL histwrite_phy(o_swsrfcf_nat, solswcf_aero(:,1)) … … 1064 1063 CALL histwrite_phy(o_swtoacf_zero,topswcf_aero(:,3)) 1065 1064 CALL histwrite_phy(o_swsrfcf_zero,solswcf_aero(:,3)) 1066 endif1067 endif! new_aod1065 ENDIF 1066 ENDIF ! new_aod 1068 1067 !====MS forcing diagnostics 1069 1068 ENDIF … … 1085 1084 ENDIF 1086 1085 ! Champs 3D: 1087 IF (ok_ade .OR. ok_aie) then1086 IF (ok_ade .OR. ok_aie) THEN 1088 1087 CALL histwrite_phy(o_ec550aer, ec550aer) 1089 1088 ENDIF … … 1170 1169 CALL histwrite_phy(o_alb2, albsol2) 1171 1170 !FH Sorties pour la couche limite 1172 if (iflag_pbl>1) then1171 IF (iflag_pbl>1) THEN 1173 1172 zx_tmp_fi3d=0. 1174 1173 IF (vars_defined) THEN 1175 donsrf=1,nbsrf1176 dok=1,klev1174 DO nsrf=1,nbsrf 1175 DO k=1,klev 1177 1176 zx_tmp_fi3d(:,k)=zx_tmp_fi3d(:,k) & 1178 1177 +pctsrf(:,nsrf)*pbl_tke(:,k,nsrf) … … 1199 1198 ENDIF 1200 1199 CALL histwrite_phy(o_dtcon, zx_tmp_fi3d) 1201 if(iflag_thermals.eq.0)then1200 IF (iflag_thermals.eq.0)THEN 1202 1201 IF (vars_defined) THEN 1203 1202 zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys + & … … 1205 1204 ENDIF 1206 1205 CALL histwrite_phy(o_tntc, zx_tmp_fi3d) 1207 else if(iflag_thermals.ge.1.and.iflag_wake.EQ.1)then1206 ELSEIF (iflag_thermals.ge.1.and.iflag_wake.EQ.1)THEN 1208 1207 IF (vars_defined) THEN 1209 1208 zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys + & … … 1212 1211 ENDIF 1213 1212 CALL histwrite_phy(o_tntc, zx_tmp_fi3d) 1214 endif1213 ENDIF 1215 1214 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_u_con(1:klon,1:klev)/pdtphys 1216 1215 CALL histwrite_phy(o_ducon, zx_tmp_fi3d) … … 1220 1219 CALL histwrite_phy(o_dqcon, zx_tmp_fi3d) 1221 1220 1222 IF (iflag_thermals.EQ.0) THEN1221 IF (iflag_thermals.EQ.0) THEN 1223 1222 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys 1224 1223 CALL histwrite_phy(o_tnhusc, zx_tmp_fi3d) 1225 ELSE IF(iflag_thermals.GE.1.AND.iflag_wake.EQ.1) THEN1224 ELSEIF (iflag_thermals.GE.1.AND.iflag_wake.EQ.1) THEN 1226 1225 IF (vars_defined) THEN 1227 1226 zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys + & … … 1243 1242 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 1244 1243 ! Sorties specifiques a la separation thermiques/non thermiques 1245 if (iflag_thermals>=1) then1246 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lscth(1:klon,1:klev)/pdtphys1244 IF (iflag_thermals>=1) THEN 1245 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lscth(1:klon,1:klev)/pdtphys 1247 1246 CALL histwrite_phy(o_dtlscth, zx_tmp_fi3d) 1248 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lscst(1:klon,1:klev)/pdtphys1247 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lscst(1:klon,1:klev)/pdtphys 1249 1248 CALL histwrite_phy(o_dtlscst, zx_tmp_fi3d) 1250 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_lscth(1:klon,1:klev)/pdtphys1249 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_lscth(1:klon,1:klev)/pdtphys 1251 1250 CALL histwrite_phy(o_dqlscth, zx_tmp_fi3d) 1252 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_lscst(1:klon,1:klev)/pdtphys1251 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_lscst(1:klon,1:klev)/pdtphys 1253 1252 CALL histwrite_phy(o_dqlscst, zx_tmp_fi3d) 1254 1253 CALL histwrite_phy(o_plulth, plul_th) 1255 1254 CALL histwrite_phy(o_plulst, plul_st) 1256 1255 IF (vars_defined) THEN 1257 dok=1,klev1258 doi=1,klon1259 if (ptconvth(i,k)) then1256 DO k=1,klev 1257 DO i=1,klon 1258 IF (ptconvth(i,k)) THEN 1260 1259 zx_tmp_fi3d(i,k)=1. 1261 else1260 ELSE 1262 1261 zx_tmp_fi3d(i,k)=0. 1263 endif1262 ENDIF 1264 1263 enddo 1265 1264 enddo … … 1267 1266 CALL histwrite_phy(o_ptconvth, zx_tmp_fi3d) 1268 1267 IF (vars_defined) THEN 1269 doi=1,klon1268 DO i=1,klon 1270 1269 zx_tmp_fi2d(1:klon)=lmax_th(:) 1271 1270 enddo 1272 1271 ENDIF 1273 1272 CALL histwrite_phy(o_lmaxth, zx_tmp_fi2d) 1274 endif! iflag_thermals>=11273 ENDIF ! iflag_thermals>=1 1275 1274 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 1276 1275 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)/pdtphys 1277 1276 CALL histwrite_phy(o_dtvdf, zx_tmp_fi3d) 1278 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_diss(1:klon,1:klev)/pdtphys1277 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_diss(1:klon,1:klev)/pdtphys 1279 1278 CALL histwrite_phy(o_dtdis, zx_tmp_fi3d) 1280 1279 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf(1:klon,1:klev)/pdtphys … … 1320 1319 CALL histwrite_phy(o_dqthe, zx_tmp_fi3d) 1321 1320 ENDIF !iflag_thermals 1322 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_ajsb(1:klon,1:klev)/pdtphys1321 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_ajsb(1:klon,1:klev)/pdtphys 1323 1322 CALL histwrite_phy(o_dtajs, zx_tmp_fi3d) 1324 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_ajsb(1:klon,1:klev)/pdtphys1323 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_ajsb(1:klon,1:klev)/pdtphys 1325 1324 CALL histwrite_phy(o_dqajs, zx_tmp_fi3d) 1326 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_swr(1:klon,1:klev)/pdtphys1325 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_swr(1:klon,1:klev)/pdtphys 1327 1326 CALL histwrite_phy(o_dtswr, zx_tmp_fi3d) 1328 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_sw0(1:klon,1:klev)/pdtphys1327 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_sw0(1:klon,1:klev)/pdtphys 1329 1328 CALL histwrite_phy(o_dtsw0, zx_tmp_fi3d) 1330 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lwr(1:klon,1:klev)/pdtphys1329 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lwr(1:klon,1:klev)/pdtphys 1331 1330 CALL histwrite_phy(o_dtlwr, zx_tmp_fi3d) 1332 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lw0(1:klon,1:klev)/pdtphys1331 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lw0(1:klon,1:klev)/pdtphys 1333 1332 CALL histwrite_phy(o_dtlw0, zx_tmp_fi3d) 1334 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)/pdtphys1333 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)/pdtphys 1335 1334 CALL histwrite_phy(o_dtec, zx_tmp_fi3d) 1336 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_u_vdf(1:klon,1:klev)/pdtphys1335 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_u_vdf(1:klon,1:klev)/pdtphys 1337 1336 CALL histwrite_phy(o_duvdf, zx_tmp_fi3d) 1338 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_v_vdf(1:klon,1:klev)/pdtphys1337 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_v_vdf(1:klon,1:klev)/pdtphys 1339 1338 CALL histwrite_phy(o_dvvdf, zx_tmp_fi3d) 1340 1339 IF (ok_orodr) THEN 1341 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_u_oro(1:klon,1:klev)/pdtphys1340 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_u_oro(1:klon,1:klev)/pdtphys 1342 1341 CALL histwrite_phy(o_duoro, zx_tmp_fi3d) 1343 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_v_oro(1:klon,1:klev)/pdtphys1342 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_v_oro(1:klon,1:klev)/pdtphys 1344 1343 CALL histwrite_phy(o_dvoro, zx_tmp_fi3d) 1345 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_oro(1:klon,1:klev)/pdtphys1344 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_oro(1:klon,1:klev)/pdtphys 1346 1345 CALL histwrite_phy(o_dtoro, zx_tmp_fi3d) 1347 1346 ENDIF 1348 1347 IF (ok_orolf) THEN 1349 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_u_lif(1:klon,1:klev)/pdtphys1348 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_u_lIF (1:klon,1:klev)/pdtphys 1350 1349 CALL histwrite_phy(o_dulif, zx_tmp_fi3d) 1351 1350 1352 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_v_lif(1:klon,1:klev)/pdtphys1351 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_v_lIF (1:klon,1:klev)/pdtphys 1353 1352 CALL histwrite_phy(o_dvlif, zx_tmp_fi3d) 1354 1353 1355 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lif(1:klon,1:klev)/pdtphys1354 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lIF (1:klon,1:klev)/pdtphys 1356 1355 CALL histwrite_phy(o_dtlif, zx_tmp_fi3d) 1357 1356 ENDIF 1358 1357 1359 1358 ! IF (ok_hines) THEN 1360 ! IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_u_hin(1:klon,1:klev)/pdtphys1359 ! IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_u_hin(1:klon,1:klev)/pdtphys 1361 1360 ! CALL histwrite_phy(o_duhin, zx_tmp_fi3d) 1362 ! IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_v_hin(1:klon,1:klev)/pdtphys1361 ! IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_v_hin(1:klon,1:klev)/pdtphys 1363 1362 ! CALL histwrite_phy(o_dvhin, zx_tmp_fi3d) 1364 ! IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_hin(1:klon,1:klev)/pdtphys1363 ! IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_hin(1:klon,1:klev)/pdtphys 1365 1364 ! CALL histwrite_phy(o_dthin, zx_tmp_fi3d) 1366 1365 ! ENDIF 1367 1366 1368 ! IF (ok_gwd_rando) then1369 ! CALL histwrite_phy(o_du_gwd_rando, du_gwd_ran do/ pdtphys)1370 ! CALL histwrite_phy(o_dv_gwd_rando, dv_gwd_ran do/ pdtphys)1367 ! IF (ok_gwd_rando) THEN 1368 ! CALL histwrite_phy(o_du_gwd_rando, du_gwd_ranDO / pdtphys) 1369 ! CALL histwrite_phy(o_dv_gwd_rando, dv_gwd_ranDO / pdtphys) 1371 1370 ! CALL histwrite_phy(o_vstr_gwd_rando, zvstr_gwd_rando) 1372 ! endIF1373 1374 IF (ok_qch4) then1371 ! ENDIF 1372 1373 IF (ok_qch4) THEN 1375 1374 CALL histwrite_phy(o_dqch4, d_q_ch4 / pdtphys) 1376 1375 ENDIF … … 1384 1383 CALL histwrite_phy(o_rlucs, lwup0) 1385 1384 CALL histwrite_phy(o_rldcs, lwdn0) 1386 IF (vars_defined) THEN1385 IF (vars_defined) THEN 1387 1386 zx_tmp_fi3d(1:klon,1:klev)=d_t(1:klon,1:klev)+ & 1388 1387 d_t_dyn(1:klon,1:klev) 1389 1388 ENDIF 1390 1389 CALL histwrite_phy(o_tnt, zx_tmp_fi3d) 1391 IF (vars_defined) THEN1390 IF (vars_defined) THEN 1392 1391 zx_tmp_fi3d(1:klon,1:klev)=d_t_swr(1:klon,1:klev)/pdtphys + & 1393 1392 d_t_lwr(1:klon,1:klev)/pdtphys 1394 1393 ENDIF 1395 1394 CALL histwrite_phy(o_tntr, zx_tmp_fi3d) 1396 IF (vars_defined) THEN1395 IF (vars_defined) THEN 1397 1396 zx_tmp_fi3d(1:klon,1:klev)= (d_t_lsc(1:klon,1:klev)+ & 1398 1397 d_t_eva(1:klon,1:klev)+ & … … 1400 1399 ENDIF 1401 1400 CALL histwrite_phy(o_tntscpbl, zx_tmp_fi3d) 1402 IF (vars_defined) THEN1401 IF (vars_defined) THEN 1403 1402 zx_tmp_fi3d(1:klon,1:klev)=d_qx(1:klon,1:klev,ivap)+ & 1404 1403 d_q_dyn(1:klon,1:klev) 1405 1404 ENDIF 1406 1405 CALL histwrite_phy(o_tnhus, zx_tmp_fi3d) 1407 IF (vars_defined) THEN1406 IF (vars_defined) THEN 1408 1407 zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys+ & 1409 1408 d_q_eva(1:klon,1:klev)/pdtphys … … 1411 1410 CALL histwrite_phy(o_tnhusscpbl, zx_tmp_fi3d) 1412 1411 CALL histwrite_phy(o_evu, coefm(:,:,is_ave)) 1413 IF (vars_defined) THEN1412 IF (vars_defined) THEN 1414 1413 zx_tmp_fi3d(1:klon,1:klev)=q_seri(1:klon,1:klev)+ & 1415 1414 ql_seri(1:klon,1:klev) 1416 1415 ENDIF 1417 1416 CALL histwrite_phy(o_h2o, zx_tmp_fi3d) 1418 if (iflag_con >= 3) then1419 IF (vars_defined) THEN1417 IF (iflag_con >= 3) THEN 1418 IF (vars_defined) THEN 1420 1419 zx_tmp_fi3d(1:klon,1:klev)=-1 * (dnwd(1:klon,1:klev)+ & 1421 1420 dnwd0(1:klon,1:klev)) 1422 1421 ENDIF 1423 1422 CALL histwrite_phy(o_mcd, zx_tmp_fi3d) 1424 IF (vars_defined) THEN1423 IF (vars_defined) THEN 1425 1424 zx_tmp_fi3d(1:klon,1:klev)=upwd(1:klon,1:klev) + & 1426 1425 dnwd(1:klon,1:klev)+ dnwd0(1:klon,1:klev) 1427 1426 ENDIF 1428 1427 CALL histwrite_phy(o_dmc, zx_tmp_fi3d) 1429 else if (iflag_con == 2) then1428 ELSEIF (iflag_con == 2) THEN 1430 1429 CALL histwrite_phy(o_mcd, pmfd) 1431 1430 CALL histwrite_phy(o_dmc, pmfu + pmfd) 1432 end if1431 ENDIF 1433 1432 CALL histwrite_phy(o_ref_liq, ref_liq) 1434 1433 CALL histwrite_phy(o_ref_ice, ref_ice) 1435 if(RCO2_per.NE.RCO2_act.OR.RCH4_per.NE.RCH4_act.OR. &1434 IF (RCO2_per.NE.RCO2_act.OR.RCH4_per.NE.RCH4_act.OR. & 1436 1435 RN2O_per.NE.RN2O_act.OR.RCFC11_per.NE.RCFC11_act.OR. & 1437 1436 RCFC12_per.NE.RCFC12_act) THEN 1438 IF (vars_defined) zx_tmp_fi2d(1 : klon) = swupp ( 1 : klon, klevp1 )1437 IF (vars_defined) zx_tmp_fi2d(1 : klon) = swupp ( 1 : klon, klevp1 ) 1439 1438 CALL histwrite_phy(o_rsut4co2, zx_tmp_fi2d) 1440 IF (vars_defined) zx_tmp_fi2d(1 : klon) = lwupp ( 1 : klon, klevp1 )1439 IF (vars_defined) zx_tmp_fi2d(1 : klon) = lwupp ( 1 : klon, klevp1 ) 1441 1440 CALL histwrite_phy(o_rlut4co2, zx_tmp_fi2d) 1442 IF (vars_defined) zx_tmp_fi2d(1 : klon) = swup0p ( 1 : klon, klevp1 )1441 IF (vars_defined) zx_tmp_fi2d(1 : klon) = swup0p ( 1 : klon, klevp1 ) 1443 1442 CALL histwrite_phy(o_rsutcs4co2, zx_tmp_fi2d) 1444 IF (vars_defined) zx_tmp_fi2d(1 : klon) = lwup0p ( 1 : klon, klevp1 )1443 IF (vars_defined) zx_tmp_fi2d(1 : klon) = lwup0p ( 1 : klon, klevp1 ) 1445 1444 CALL histwrite_phy(o_rlutcs4co2, zx_tmp_fi2d) 1446 1445 CALL histwrite_phy(o_rsu4co2, swupp) … … 1460 1459 ! Champs interpolles sur des niveaux de pression 1461 1460 missing_val=missing_val_nf90 1462 DO iff=7, nfiles 1461 DO iff=7, nfiles-1 !--here we deal with files 7,8 and 9 1463 1462 1464 1463 CALL histwrite_phy(o_tnondef,tnondef(:,:,iff-6),iff) … … 1470 1469 CALL histwrite_phy(o_va,vwriteSTD(:,:,iff-6),iff) 1471 1470 CALL histwrite_phy(o_wap,wwriteSTD(:,:,iff-6),iff) 1472 IF (vars_defined) THEN1471 IF (vars_defined) THEN 1473 1472 DO k=1, nlevSTD 1474 1473 DO i=1, klon 1475 IF (tnondef(i,k,iff-6).NE.missing_val) THEN1476 IF (freq_outNMC(iff-6).LT.0) THEN1474 IF (tnondef(i,k,iff-6).NE.missing_val) THEN 1475 IF (freq_outNMC(iff-6).LT.0) THEN 1477 1476 freq_moyNMC(iff-6)=(mth_len*un_jour)/freq_calNMC(iff-6) 1478 1477 ELSE … … 1487 1486 ENDIF 1488 1487 CALL histwrite_phy(o_psbg,zx_tmp_fi3d_STD,iff) 1489 IF (vars_defined) THEN1488 IF (vars_defined) THEN 1490 1489 DO k=1, nlevSTD 1491 1490 DO i=1, klon 1492 IF (O3sumSTD(i,k,iff-6).NE.missing_val) THEN1491 IF (O3sumSTD(i,k,iff-6).NE.missing_val) THEN 1493 1492 zx_tmp_fi3d_STD(i,k) = O3sumSTD(i,k,iff-6) * 1.e+9 1494 1493 ELSE … … 1499 1498 ENDIF 1500 1499 CALL histwrite_phy(o_tro3,zx_tmp_fi3d_STD,iff) 1501 if(read_climoz == 2) THEN1502 IF (vars_defined) THEN1500 IF (read_climoz == 2) THEN 1501 IF (vars_defined) THEN 1503 1502 DO k=1, nlevSTD 1504 1503 DO i=1, klon 1505 IF (O3daysumSTD(i,k,iff-6).NE.missing_val) THEN1504 IF (O3daysumSTD(i,k,iff-6).NE.missing_val) THEN 1506 1505 zx_tmp_fi3d_STD(i,k) = O3daysumSTD(i,k,iff-6) * 1.e+9 1507 1506 ELSE … … 1512 1511 ENDIF 1513 1512 CALL histwrite_phy(o_tro3_daylight,zx_tmp_fi3d_STD,iff) 1514 endif1513 ENDIF 1515 1514 CALL histwrite_phy(o_uxv,uvsumSTD(:,:,iff-6),iff) 1516 1515 CALL histwrite_phy(o_vxq,vqsumSTD(:,:,iff-6),iff) … … 1527 1526 #endif 1528 1527 #ifdef CPP_XIOS 1529 IF (ok_all_xml) THEN1528 IF (ok_all_xml) THEN 1530 1529 ! DO iff=7, nfiles 1531 1530 … … 1538 1537 CALL histwrite_phy(o_va,vlevSTD(:,:)) 1539 1538 CALL histwrite_phy(o_wap,wlevSTD(:,:)) 1540 ! IF (vars_defined) THEN1539 ! IF (vars_defined) THEN 1541 1540 ! DO k=1, nlevSTD 1542 1541 ! DO i=1, klon 1543 ! IF (tnondef(i,k,3).NE.missing_val) THEN1544 ! IF (freq_outNMC(iff-6).LT.0) THEN1542 ! IF (tnondef(i,k,3).NE.missing_val) THEN 1543 ! IF (freq_outNMC(iff-6).LT.0) THEN 1545 1544 ! freq_moyNMC(iff-6)=(mth_len*un_jour)/freq_calNMC(iff-6) 1546 1545 ! ELSE … … 1555 1554 ! ENDIF 1556 1555 ! CALL histwrite_phy(o_psbg,zx_tmp_fi3d_STD) 1557 IF (vars_defined) THEN1556 IF (vars_defined) THEN 1558 1557 DO k=1, nlevSTD 1559 1558 DO i=1, klon 1560 IF (O3STD(i,k).NE.missing_val) THEN1559 IF (O3STD(i,k).NE.missing_val) THEN 1561 1560 zx_tmp_fi3d_STD(i,k) = O3STD(i,k) * 1.e+9 1562 1561 ELSE … … 1567 1566 ENDIF 1568 1567 CALL histwrite_phy(o_tro3,zx_tmp_fi3d_STD) 1569 if(read_climoz == 2) THEN1570 IF (vars_defined) THEN1568 IF (read_climoz == 2) THEN 1569 IF (vars_defined) THEN 1571 1570 DO k=1, nlevSTD 1572 1571 DO i=1, klon 1573 IF (O3daySTD(i,k).NE.missing_val) THEN1572 IF (O3daySTD(i,k).NE.missing_val) THEN 1574 1573 zx_tmp_fi3d_STD(i,k) = O3daySTD(i,k) * 1.e+9 1575 1574 ELSE … … 1580 1579 ENDIF 1581 1580 CALL histwrite_phy(o_tro3_daylight,zx_tmp_fi3d_STD) 1582 endif1581 ENDIF 1583 1582 CALL histwrite_phy(o_uxv,uvSTD(:,:)) 1584 1583 CALL histwrite_phy(o_vxq,vqSTD(:,:)) … … 1598 1597 IF (type_trac == 'lmdz' .OR. type_trac == 'repr') THEN 1599 1598 1600 !jyg<1601 !! CALL histwrite_phy(o_trac(iq-nqo), qx(:,:,iq))1602 1599 CALL histwrite_phy(o_trac(iq-nqo), tr_seri(:,:,iq-nqo)) 1603 !>jyg1604 1600 CALL histwrite_phy(o_dtr_vdf(iq-nqo),d_tr_cl(:,:,iq-nqo)) 1605 1601 CALL histwrite_phy(o_dtr_the(iq-nqo),d_tr_th(:,:,iq-nqo)) … … 1612 1608 CALL histwrite_phy(o_dtr_ls(iq-nqo),d_tr_ls(:,:,iq-nqo)) 1613 1609 ! CALL histwrite_phy(o_dtr_dyn(iq-nqo),d_tr_dyn(:,:,iq-nqo)) 1614 !!CALL histwrite_phy(o_dtr_cl(iq-nqo),d_tr_cl(:,:,iq-nqo))1610 ! CALL histwrite_phy(o_dtr_cl(iq-nqo),d_tr_cl(:,:,iq-nqo)) 1615 1611 CALL histwrite_phy(o_dtr_trsp(iq-nqo),d_tr_trsp(:,:,iq-nqo)) 1616 1612 CALL histwrite_phy(o_dtr_sscav(iq-nqo),d_tr_sscav(:,:,iq-nqo)) … … 1618 1614 CALL histwrite_phy(o_dtr_uscav(iq-nqo),d_tr_uscav(:,:,iq-nqo)) 1619 1615 zx_tmp_fi2d=0. 1620 IF (vars_defined) THEN1616 IF (vars_defined) THEN 1621 1617 DO k=1,klev 1622 !jyg<1623 !! zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+zmasse(:,k)*qx(:,k,iq)1624 1618 zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+zmasse(:,k)*tr_seri(:,k,iq-nqo) 1625 !>jyg1626 1619 ENDDO 1627 1620 ENDIF 1628 1621 CALL histwrite_phy(o_trac_cum(iq-nqo), zx_tmp_fi2d) 1629 endif1630 ENDDO 1631 ENDIF 1632 1633 IF (.NOT.vars_defined) THEN1622 ENDIF 1623 ENDDO 1624 ENDIF 1625 1626 IF (.NOT.vars_defined) THEN 1634 1627 !$OMP MASTER 1635 1628 #ifndef CPP_IOIPSL_NO_OUTPUT … … 1652 1645 vars_defined = .TRUE. 1653 1646 1654 END 1655 1656 END 1657 1658 IF (vars_defined) THEN1647 ENDIF 1648 1649 ENDDO 1650 1651 IF (vars_defined) THEN 1659 1652 ! On synchronise les fichiers pour IOIPSL 1660 1653 #ifndef CPP_IOIPSL_NO_OUTPUT … … 1664 1657 CALL histsync(nid_files(iff)) 1665 1658 ENDIF 1666 END 1659 ENDDO 1667 1660 !$OMP END MASTER 1668 1661 #endif -
LMDZ5/branches/testing/libf/phylmd/Dust/spla_output_dat.h
r2641 r2787 3 3 4 4 !Dust emission module 5 type(ctrl_out),save :: o_m1dflux = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &6 'm1dflux','m1dflux','', (/ ('', i=1, 9) /))7 type(ctrl_out),save :: o_m2dflux = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &8 'm2dflux','m2dflux','', (/ ('', i=1, 9) /))9 type(ctrl_out),save :: o_m3dflux = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &10 'm3dflux','m3dflux','', (/ ('', i=1, 9) /))5 type(ctrl_out),save :: o_m1dflux = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 6 'm1dflux','m1dflux','', (/ ('', i=1, 10) /)) 7 type(ctrl_out),save :: o_m2dflux = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 8 'm2dflux','m2dflux','', (/ ('', i=1, 10) /)) 9 type(ctrl_out),save :: o_m3dflux = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 10 'm3dflux','m3dflux','', (/ ('', i=1, 10) /)) 11 11 12 12 ! traceur_spl 13 type(ctrl_out),save :: o_taue550 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 14 'taue550','Tau ext 550','', (/ ('', i=1, 9) /)) 15 type(ctrl_out),save :: o_taue670 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 16 'taue670','Tau ext 670','', (/ ('', i=1, 9) /)) 17 type(ctrl_out),save :: o_taue865 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 18 'taue865','Tau ext 865','', (/ ('', i=1, 9) /)) 19 type(ctrl_out),save :: o_taue550_tr2 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 20 'taue550_tr2','Tau ext 550tr2','', (/ ('', i=1, 9) /)) 21 type(ctrl_out),save :: o_taue670_tr2 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 22 'taue670_tr2','Tau ext 670tr2','', (/ ('', i=1, 9) /)) 23 type(ctrl_out),save :: o_taue865_tr2 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 24 'taue865_tr2','Tau ext 865tr2','', (/ ('', i=1, 9) /)) 25 type(ctrl_out),save :: o_taue550_ss = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 26 'taue550_ss','Tau ext 550ss','', (/ ('', i=1, 9) /)) 27 type(ctrl_out),save :: o_taue670_ss = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 28 'taue670_ss','Tau ext 670ss','', (/ ('', i=1, 9) /)) 29 type(ctrl_out),save :: o_taue865_ss = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 30 'taue865_ss','Tau ext 865ss','', (/ ('', i=1, 9) /)) 31 type(ctrl_out),save :: o_taue550_dust = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 32 'taue550_dust','Tau ext 550dust','', (/ ('', i=1, 9) /)) 33 type(ctrl_out),save :: o_taue670_dust = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 34 'taue670_dust','Tau ext 670dust','', (/ ('', i=1, 9) /)) 35 type(ctrl_out),save :: o_taue865_dust = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 36 'taue865_dust','Tau ext 865dust','', (/ ('', i=1, 9) /)) 37 type(ctrl_out),save :: o_taue550_dustsco = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 38 'taue550_dustsco','Tau ext 550dustsco','', (/ ('', i=1, 9) /)) 39 type(ctrl_out),save :: o_taue670_dustsco = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 40 'taue670_dustsco','Tau ext 670dustsco','', (/ ('', i=1, 9) /)) 41 type(ctrl_out),save :: o_taue865_dustsco = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 42 'taue865_dustsco','Tau ext 865dustsco','', (/ ('', i=1, 9) /)) 43 44 45 type(ctrl_out),save :: o_taue550_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 13 type(ctrl_out),save :: o_taue550 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 14 'taue550','Tau ext 550','', (/ ('', i=1, 10) /)) 15 type(ctrl_out),save :: o_taue670 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 16 'taue670','Tau ext 670','', (/ ('', i=1, 10) /)) 17 type(ctrl_out),save :: o_taue865 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 18 'taue865','Tau ext 865','', (/ ('', i=1, 10) /)) 19 type(ctrl_out),save :: o_taue550_tr2 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 20 'taue550_tr2','Tau ext 550tr2','', (/ ('', i=1, 10) /)) 21 type(ctrl_out),save :: o_taue670_tr2 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 22 'taue670_tr2','Tau ext 670tr2','', (/ ('', i=1, 10) /)) 23 type(ctrl_out),save :: o_taue865_tr2 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 24 'taue865_tr2','Tau ext 865tr2','', (/ ('', i=1, 10) /)) 25 type(ctrl_out),save :: o_taue550_ss = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 26 'taue550_ss','Tau ext 550ss','', (/ ('', i=1, 10) /)) 27 type(ctrl_out),save :: o_taue670_ss = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 28 'taue670_ss','Tau ext 670ss','', (/ ('', i=1, 10) /)) 29 type(ctrl_out),save :: o_taue865_ss = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 30 'taue865_ss','Tau ext 865ss','', (/ ('', i=1, 10) /)) 31 type(ctrl_out),save :: o_taue550_dust = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 32 'taue550_dust','Tau ext 550dust','', (/ ('', i=1, 10) /)) 33 type(ctrl_out),save :: o_taue670_dust = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 34 'taue670_dust','Tau ext 670dust','', (/ ('', i=1, 10) /)) 35 type(ctrl_out),save :: o_taue865_dust = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 36 'taue865_dust','Tau ext 865dust','', (/ ('', i=1, 10) /)) 37 type(ctrl_out),save :: o_taue550_dustsco = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 38 'taue550_dustsco','Tau ext 550dustsco','', (/ ('', i=1, 10) /)) 39 type(ctrl_out),save :: o_taue670_dustsco = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 40 'taue670_dustsco','Tau ext 670dustsco','', (/ ('', i=1, 10) /)) 41 type(ctrl_out),save :: o_taue865_dustsco = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 42 'taue865_dustsco','Tau ext 865dustsco','', (/ ('', i=1, 10) /)) 43 44 type(ctrl_out),save :: o_taue550_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 46 45 'taue550_terra','Tau ext 550 terra','', & 47 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &48 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 49 type(ctrl_out),save :: o_taue550_fine_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &46 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 47 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 48 type(ctrl_out),save :: o_taue550_fine_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 50 49 'taue550_fine_terra','Tau ext fine 550 terra','', & 51 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &52 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 53 type(ctrl_out),save :: o_taue550_coss_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &50 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 51 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 52 type(ctrl_out),save :: o_taue550_coss_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 54 53 'taue550_coss_terra','Tau ext coss 550 terra','', & 55 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &56 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 57 type(ctrl_out),save :: o_taue550_codu_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &54 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 55 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 56 type(ctrl_out),save :: o_taue550_codu_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 58 57 'taue550_codu_terra','Tau ext codu 550 terra','', & 59 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &60 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 61 type(ctrl_out),save :: o_taue550_scdu_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &58 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 59 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 60 type(ctrl_out),save :: o_taue550_scdu_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 62 61 'taue550_scdu_terra','Tau ext scdu 550 terra','', & 63 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &64 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 65 66 type(ctrl_out),save :: o_taue670_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &62 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 63 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 64 65 type(ctrl_out),save :: o_taue670_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 67 66 'taue670_terra','Tau ext 670 terra','', & 68 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &69 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 70 type(ctrl_out),save :: o_taue670_fine_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &67 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 68 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 69 type(ctrl_out),save :: o_taue670_fine_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 71 70 'taue670_fine_terra','Tau ext fine 670 terra','', & 72 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &73 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 74 type(ctrl_out),save :: o_taue670_coss_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &71 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 72 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 73 type(ctrl_out),save :: o_taue670_coss_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 75 74 'taue670_coss_terra','Tau ext coss 670 terra','', & 76 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &77 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 78 type(ctrl_out),save :: o_taue670_codu_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &75 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 76 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 77 type(ctrl_out),save :: o_taue670_codu_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 79 78 'taue670_codu_terra','Tau ext codu 670 terra','', & 80 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &81 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 82 type(ctrl_out),save :: o_taue670_scdu_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &79 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 80 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 81 type(ctrl_out),save :: o_taue670_scdu_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 83 82 'taue670_scdu_terra','Tau ext scdu 670 terra','', & 84 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &85 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 86 87 type(ctrl_out),save :: o_taue865_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &83 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 84 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 85 86 type(ctrl_out),save :: o_taue865_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 88 87 'taue865_terra','Tau ext 865 terra','', & 89 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &90 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 91 type(ctrl_out),save :: o_taue865_fine_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &88 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 89 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 90 type(ctrl_out),save :: o_taue865_fine_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 92 91 'taue865_fine_terra','Tau ext fine 865 terra','', & 93 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &94 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 95 type(ctrl_out),save :: o_taue865_coss_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &92 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 93 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 94 type(ctrl_out),save :: o_taue865_coss_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 96 95 'taue865_coss_terra','Tau ext coss 865 terra','', & 97 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &98 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 99 type(ctrl_out),save :: o_taue865_codu_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &96 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 97 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 98 type(ctrl_out),save :: o_taue865_codu_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 100 99 'taue865_codu_terra','Tau ext codu 865 terra','', & 101 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &102 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 103 type(ctrl_out),save :: o_taue865_scdu_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &100 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 101 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 102 type(ctrl_out),save :: o_taue865_scdu_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 104 103 'taue865_scdu_terra','Tau ext scdu 865 terra','', & 105 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &106 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 107 108 type(ctrl_out),save :: o_taue550_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &104 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 105 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 106 107 type(ctrl_out),save :: o_taue550_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 109 108 'taue550_aqua','Tau ext 550 aqua','', & 110 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &111 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 112 type(ctrl_out),save :: o_taue550_fine_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &109 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 110 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 111 type(ctrl_out),save :: o_taue550_fine_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 113 112 'taue550_fine_aqua','Tau ext fine 550 aqua','', & 114 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &115 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 116 type(ctrl_out),save :: o_taue550_coss_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &113 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 114 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 115 type(ctrl_out),save :: o_taue550_coss_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 117 116 'taue550_coss_aqua','Tau ext coss 550 aqua','', & 118 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &119 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 120 type(ctrl_out),save :: o_taue550_codu_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &117 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 118 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 119 type(ctrl_out),save :: o_taue550_codu_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 121 120 'taue550_codu_aqua','Tau ext codu 550 aqua','', & 122 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &123 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 124 type(ctrl_out),save :: o_taue550_scdu_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &121 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 122 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 123 type(ctrl_out),save :: o_taue550_scdu_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 125 124 'taue550_scdu_aqua','Tau ext scdu 550 aqua','', & 126 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &127 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 128 129 type(ctrl_out),save :: o_taue670_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &125 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 126 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 127 128 type(ctrl_out),save :: o_taue670_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 130 129 'taue670_aqua','Tau ext 670 aqua','', & 131 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &132 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 133 type(ctrl_out),save :: o_taue670_fine_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &130 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 131 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 132 type(ctrl_out),save :: o_taue670_fine_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 134 133 'taue670_fine_aqua','Tau ext fine 670 aqua','', & 135 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &136 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 137 type(ctrl_out),save :: o_taue670_coss_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &134 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 135 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 136 type(ctrl_out),save :: o_taue670_coss_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 138 137 'taue670_coss_aqua','Tau ext coss 670 aqua','', & 139 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &140 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 141 type(ctrl_out),save :: o_taue670_codu_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &138 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 139 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 140 type(ctrl_out),save :: o_taue670_codu_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 142 141 'taue670_codu_aqua','Tau ext codu 670 aqua','', & 143 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &144 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 145 type(ctrl_out),save :: o_taue670_scdu_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &142 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 143 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 144 type(ctrl_out),save :: o_taue670_scdu_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 146 145 'taue670_scdu_aqua','Tau ext scdu 670 aqua','', & 147 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &148 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 149 150 type(ctrl_out),save :: o_taue865_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &146 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 147 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 148 149 type(ctrl_out),save :: o_taue865_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 151 150 'taue865_aqua','Tau ext 865 aqua','', & 152 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &153 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 154 type(ctrl_out),save :: o_taue865_fine_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &151 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 152 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 153 type(ctrl_out),save :: o_taue865_fine_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 155 154 'taue865_fine_aqua','Tau ext fine 865 aqua','', & 156 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &157 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 158 type(ctrl_out),save :: o_taue865_coss_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &155 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 156 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 157 type(ctrl_out),save :: o_taue865_coss_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 159 158 'taue865_coss_aqua','Tau ext coss 865 aqua','', & 160 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &161 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 162 type(ctrl_out),save :: o_taue865_codu_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &159 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 160 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 161 type(ctrl_out),save :: o_taue865_codu_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 163 162 'taue865_codu_aqua','Tau ext codu 865 aqua','', & 164 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &165 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 166 type(ctrl_out),save :: o_taue865_scdu_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &163 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 164 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 165 type(ctrl_out),save :: o_taue865_scdu_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 167 166 'taue865_scdu_aqua','Tau ext scdu 865 aqua','', & 168 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 169 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 170 171 172 173 ! type(ctrl_out),save :: o_taue550_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 167 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 168 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 169 170 ! type(ctrl_out),save :: o_taue550_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 174 171 ! 'taue550_terra','Tau ext 550 terra','', & 175 ! (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &172 ! (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 176 173 ! 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 177 174 ! 178 ! type(ctrl_out),save :: o_taue670_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &175 ! type(ctrl_out),save :: o_taue670_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 179 176 ! 'taue670_aqua','Tau ext 670 aqua','', & 180 ! (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &177 ! (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 181 178 ! 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 182 179 ! 183 ! type(ctrl_out),save :: o_taue670_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &180 ! type(ctrl_out),save :: o_taue670_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 184 181 ! 'taue670_terra','Tau ext 670 terra','', & 185 ! (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &182 ! (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 186 183 ! 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 187 184 ! 188 ! type(ctrl_out),save :: o_taue865_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &185 ! type(ctrl_out),save :: o_taue865_aqua = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 189 186 ! 'taue865_aqua','Tau ext 865 aqua','', & 190 ! (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &187 ! (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 191 188 ! 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 192 189 ! 193 ! type(ctrl_out),save :: o_taue865_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &190 ! type(ctrl_out),save :: o_taue865_terra = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 194 191 ! 'taue865_terra','Tau ext 865 terra','', & 195 ! (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &192 ! (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 196 193 ! 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 197 198 194 199 type(ctrl_out),save :: o_trm01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 200 'trm01','Burden PREC','', (/ ('', i=1, 9) /)) 201 type(ctrl_out),save :: o_trm02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 202 'trm02','Burden FINE','', (/ ('', i=1, 9) /)) 203 type(ctrl_out),save :: o_trm03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 204 'trm03','Burden COSS','', (/ ('', i=1, 9) /)) 205 type(ctrl_out),save :: o_trm04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 206 'trm04','Burden CODU','', (/ ('', i=1, 9) /)) 207 type(ctrl_out),save :: o_trm05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 208 'trm05','Burden SCDU','', (/ ('', i=1, 9) /)) 209 210 type(ctrl_out),save :: o_sconc01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 211 'sconc01','Surf. Conc. PREC','', (/ ('', i=1, 9) /)) 212 type(ctrl_out),save :: o_sconc02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 213 'sconc02','Surf. Conc. FINE','', (/ ('', i=1, 9) /)) 214 type(ctrl_out),save :: o_sconc03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 215 'sconc03','Surf. Conc. COSS','', (/ ('', i=1, 9) /)) 216 type(ctrl_out),save :: o_sconc04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 217 'sconc04','Surf. Conc. CODU','', (/ ('', i=1, 9) /)) 218 type(ctrl_out),save :: o_sconc05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 219 'sconc05','Surf. Conc. SCDU','', (/ ('', i=1, 9) /)) 220 221 222 195 type(ctrl_out),save :: o_trm01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 196 'trm01','Burden PREC','', (/ ('', i=1, 10) /)) 197 type(ctrl_out),save :: o_trm02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 198 'trm02','Burden FINE','', (/ ('', i=1, 10) /)) 199 type(ctrl_out),save :: o_trm03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 200 'trm03','Burden COSS','', (/ ('', i=1, 10) /)) 201 type(ctrl_out),save :: o_trm04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 202 'trm04','Burden CODU','', (/ ('', i=1, 10) /)) 203 type(ctrl_out),save :: o_trm05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 204 'trm05','Burden SCDU','', (/ ('', i=1, 10) /)) 205 206 type(ctrl_out),save :: o_sconc01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 207 'sconc01','Surf. Conc. PREC','', (/ ('', i=1, 10) /)) 208 type(ctrl_out),save :: o_sconc02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 209 'sconc02','Surf. Conc. FINE','', (/ ('', i=1, 10) /)) 210 type(ctrl_out),save :: o_sconc03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 211 'sconc03','Surf. Conc. COSS','', (/ ('', i=1, 10) /)) 212 type(ctrl_out),save :: o_sconc04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 213 'sconc04','Surf. Conc. CODU','', (/ ('', i=1, 10) /)) 214 type(ctrl_out),save :: o_sconc05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 215 'sconc05','Surf. Conc. SCDU','', (/ ('', i=1, 10) /)) 223 216 224 217 !lessivage 225 218 226 type(ctrl_out),save :: o_flux01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 227 'flux01','emission PREC','', (/ ('', i=1, 9) /)) 228 type(ctrl_out),save :: o_flux02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 229 'flux02','emission FINE','', (/ ('', i=1, 9) /)) 230 type(ctrl_out),save :: o_flux03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 231 'flux03','emission COSS','', (/ ('', i=1, 9) /)) 232 type(ctrl_out),save :: o_flux04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 233 'flux04','emission CODU','', (/ ('', i=1, 9) /)) 234 type(ctrl_out),save :: o_flux05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 235 'flux05','emission SCDU','', (/ ('', i=1, 9) /)) 236 237 238 type(ctrl_out),save :: o_ds01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 239 'ds01','Depot sec PREC','', (/ ('', i=1, 9) /)) 240 type(ctrl_out),save :: o_ds02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 241 'ds02','Depot sec FINE','', (/ ('', i=1, 9) /)) 242 type(ctrl_out),save :: o_ds03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 243 'ds03','Depot sec COSS','', (/ ('', i=1, 9) /)) 244 type(ctrl_out),save :: o_ds04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 245 'ds04','Depot sec CODU','', (/ ('', i=1, 9) /)) 246 type(ctrl_out),save :: o_ds05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 247 'ds05','Depot sec SCDU','', (/ ('', i=1, 9) /)) 248 249 type(ctrl_out),save :: o_dh01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 250 'dh01','Depot hum PREC','', (/ ('', i=1, 9) /)) 251 type(ctrl_out),save :: o_dh02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 252 'dh02','Depot hum FINE','', (/ ('', i=1, 9) /)) 253 type(ctrl_out),save :: o_dh03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 254 'dh03','Depot hum COSS','', (/ ('', i=1, 9) /)) 255 type(ctrl_out),save :: o_dh04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 256 'dh04','Depot hum CODU','', (/ ('', i=1, 9) /)) 257 type(ctrl_out),save :: o_dh05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 258 'dh05','Depot hum SCDU','', (/ ('', i=1, 9) /)) 259 260 type(ctrl_out),save :: o_dtrconv01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 261 'dtrconv01','Tiedke convective PREC','', (/ ('', i=1, 9) /)) 262 type(ctrl_out),save :: o_dtrconv02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 263 'dtrconv02','Tiedke convective FINE','', (/ ('', i=1, 9) /)) 264 type(ctrl_out),save :: o_dtrconv03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 265 'dtrconv03','Tiedke convective COSS','', (/ ('', i=1, 9) /)) 266 type(ctrl_out),save :: o_dtrconv04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 267 'dtrconv04','Tiedke convective CODU','', (/ ('', i=1, 9) /)) 268 type(ctrl_out),save :: o_dtrconv05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 269 'dtrconv05','Tiedke convective SCDU','', (/ ('', i=1, 9) /)) 270 271 type(ctrl_out),save :: o_dtherm01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 272 'dtherm01','Thermals dtracer PREC','', (/ ('', i=1, 9) /)) 273 type(ctrl_out),save :: o_dtherm02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 274 'dtherm02','Thermals dtracer FINE','', (/ ('', i=1, 9) /)) 275 type(ctrl_out),save :: o_dtherm03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 276 'dtherm03','Thermals dtracer COSS','', (/ ('', i=1, 9) /)) 277 type(ctrl_out),save :: o_dtherm04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 278 'dtherm04','Thermals dtracer CODU','', (/ ('', i=1, 9) /)) 279 type(ctrl_out),save :: o_dtherm05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 280 'dtherm05','Thermals dtracer SCDU','', (/ ('', i=1, 9) /)) 281 282 type(ctrl_out),save :: o_dhkecv01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 283 'dhkecv01','KE dep hum convective PREC','', (/ ('', i=1, 9) /)) 284 type(ctrl_out),save :: o_dhkecv02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 285 'dhkecv02','KE dep hum convective FINE','', (/ ('', i=1, 9) /)) 286 type(ctrl_out),save :: o_dhkecv03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 287 'dhkecv03','KE dep hum convective COSS','', (/ ('', i=1, 9) /)) 288 type(ctrl_out),save :: o_dhkecv04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 289 'dhkecv04','KE dep hum convective CODU','', (/ ('', i=1, 9) /)) 290 type(ctrl_out),save :: o_dhkecv05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 291 'dhkecv05','KE dep hum convective SCDU','', (/ ('', i=1, 9) /)) 219 type(ctrl_out),save :: o_flux01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 220 'flux01','emission PREC','', (/ ('', i=1, 10) /)) 221 type(ctrl_out),save :: o_flux02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 222 'flux02','emission FINE','', (/ ('', i=1, 10) /)) 223 type(ctrl_out),save :: o_flux03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 224 'flux03','emission COSS','', (/ ('', i=1, 10) /)) 225 type(ctrl_out),save :: o_flux04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 226 'flux04','emission CODU','', (/ ('', i=1, 10) /)) 227 type(ctrl_out),save :: o_flux05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 228 'flux05','emission SCDU','', (/ ('', i=1, 10) /)) 229 230 type(ctrl_out),save :: o_ds01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 231 'ds01','Depot sec PREC','', (/ ('', i=1, 10) /)) 232 type(ctrl_out),save :: o_ds02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 233 'ds02','Depot sec FINE','', (/ ('', i=1, 10) /)) 234 type(ctrl_out),save :: o_ds03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 235 'ds03','Depot sec COSS','', (/ ('', i=1, 10) /)) 236 type(ctrl_out),save :: o_ds04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 237 'ds04','Depot sec CODU','', (/ ('', i=1, 10) /)) 238 type(ctrl_out),save :: o_ds05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 239 'ds05','Depot sec SCDU','', (/ ('', i=1, 10) /)) 240 241 type(ctrl_out),save :: o_dh01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 242 'dh01','Depot hum PREC','', (/ ('', i=1, 10) /)) 243 type(ctrl_out),save :: o_dh02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 244 'dh02','Depot hum FINE','', (/ ('', i=1, 10) /)) 245 type(ctrl_out),save :: o_dh03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 246 'dh03','Depot hum COSS','', (/ ('', i=1, 10) /)) 247 type(ctrl_out),save :: o_dh04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 248 'dh04','Depot hum CODU','', (/ ('', i=1, 10) /)) 249 type(ctrl_out),save :: o_dh05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 250 'dh05','Depot hum SCDU','', (/ ('', i=1, 10) /)) 251 252 type(ctrl_out),save :: o_dtrconv01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 253 'dtrconv01','Tiedke convective PREC','', (/ ('', i=1, 10) /)) 254 type(ctrl_out),save :: o_dtrconv02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 255 'dtrconv02','Tiedke convective FINE','', (/ ('', i=1, 10) /)) 256 type(ctrl_out),save :: o_dtrconv03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 257 'dtrconv03','Tiedke convective COSS','', (/ ('', i=1, 10) /)) 258 type(ctrl_out),save :: o_dtrconv04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 259 'dtrconv04','Tiedke convective CODU','', (/ ('', i=1, 10) /)) 260 type(ctrl_out),save :: o_dtrconv05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 261 'dtrconv05','Tiedke convective SCDU','', (/ ('', i=1, 10) /)) 262 263 type(ctrl_out),save :: o_dtherm01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 264 'dtherm01','Thermals dtracer PREC','', (/ ('', i=1, 10) /)) 265 type(ctrl_out),save :: o_dtherm02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 266 'dtherm02','Thermals dtracer FINE','', (/ ('', i=1, 10) /)) 267 type(ctrl_out),save :: o_dtherm03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 268 'dtherm03','Thermals dtracer COSS','', (/ ('', i=1, 10) /)) 269 type(ctrl_out),save :: o_dtherm04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 270 'dtherm04','Thermals dtracer CODU','', (/ ('', i=1, 10) /)) 271 type(ctrl_out),save :: o_dtherm05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 272 'dtherm05','Thermals dtracer SCDU','', (/ ('', i=1, 10) /)) 273 274 type(ctrl_out),save :: o_dhkecv01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 275 'dhkecv01','KE dep hum convective PREC','', (/ ('', i=1, 10) /)) 276 type(ctrl_out),save :: o_dhkecv02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 277 'dhkecv02','KE dep hum convective FINE','', (/ ('', i=1, 10) /)) 278 type(ctrl_out),save :: o_dhkecv03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 279 'dhkecv03','KE dep hum convective COSS','', (/ ('', i=1, 10) /)) 280 type(ctrl_out),save :: o_dhkecv04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 281 'dhkecv04','KE dep hum convective CODU','', (/ ('', i=1, 10) /)) 282 type(ctrl_out),save :: o_dhkecv05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 283 'dhkecv05','KE dep hum convective SCDU','', (/ ('', i=1, 10) /)) 292 284 293 285 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 294 type(ctrl_out),save :: o_d_tr_ds01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &295 'd_tr_ds01','Tendance dep sec PREC','', (/ ('', i=1, 9) /))296 type(ctrl_out),save :: o_d_tr_ds02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &297 'd_tr_ds02','Tendance dep sec FINE','', (/ ('', i=1, 9) /))298 type(ctrl_out),save :: o_d_tr_ds03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &299 'd_tr_ds03','Tendance dep sec COSS','', (/ ('', i=1, 9) /))300 type(ctrl_out),save :: o_d_tr_ds04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &301 'd_tr_ds04','Tendance depot sec CODU','', (/ ('', i=1, 9) /))302 type(ctrl_out),save :: o_d_tr_ds05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &303 'd_tr_ds05','Tendance dep sec SCDU','', (/ ('', i=1, 9) /))286 type(ctrl_out),save :: o_d_tr_ds01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 287 'd_tr_ds01','Tendance dep sec PREC','', (/ ('', i=1, 10) /)) 288 type(ctrl_out),save :: o_d_tr_ds02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 289 'd_tr_ds02','Tendance dep sec FINE','', (/ ('', i=1, 10) /)) 290 type(ctrl_out),save :: o_d_tr_ds03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 291 'd_tr_ds03','Tendance dep sec COSS','', (/ ('', i=1, 10) /)) 292 type(ctrl_out),save :: o_d_tr_ds04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 293 'd_tr_ds04','Tendance depot sec CODU','', (/ ('', i=1, 10) /)) 294 type(ctrl_out),save :: o_d_tr_ds05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 295 'd_tr_ds05','Tendance dep sec SCDU','', (/ ('', i=1, 10) /)) 304 296 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 305 297 306 type(ctrl_out),save :: o_dhkelsc01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 307 'dhkelsc01','KE dep hum large scale PREC','', (/ ('', i=1, 9) /)) 308 type(ctrl_out),save :: o_dhkelsc02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 309 'dhkelsc02','KE dep hum large scale FINE','', (/ ('', i=1, 9) /)) 310 type(ctrl_out),save :: o_dhkelsc03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 311 'dhkelsc03','KE dep hum large scale COSS','', (/ ('', i=1, 9) /)) 312 type(ctrl_out),save :: o_dhkelsc04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 313 'dhkelsc04','KE dep hum large scale CODU','', (/ ('', i=1, 9) /)) 314 type(ctrl_out),save :: o_dhkelsc05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 315 'dhkelsc05','KE dep hum large scale SCDU','', (/ ('', i=1, 9) /)) 316 317 318 type(ctrl_out),save :: o_d_tr_cv01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 319 'd_tr_cv01','cvltr d_tr_cv PREC','', (/ ('', i=1, 9) /)) 320 type(ctrl_out),save :: o_d_tr_cv02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 321 'd_tr_cv02','cvltr d_tr_cv FINE','', (/ ('', i=1, 9) /)) 322 type(ctrl_out),save :: o_d_tr_cv03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 323 'd_tr_cv03','cvltr d_tr_cv COSS','', (/ ('', i=1, 9) /)) 324 type(ctrl_out),save :: o_d_tr_cv04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 325 'd_tr_cv04','cvltr d_tr_cv CODU','', (/ ('', i=1, 9) /)) 326 type(ctrl_out),save :: o_d_tr_cv05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 327 'd_tr_cv05','cvltr d_tr_cv SCDU','', (/ ('', i=1, 9) /)) 328 329 330 type(ctrl_out),save :: o_d_tr_trsp01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 331 'd_tr_trsp01','cvltr d_tr_trsp PREC','', (/ ('', i=1, 9) /)) 332 type(ctrl_out),save :: o_d_tr_trsp02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 333 'd_tr_trsp02','cvltr d_tr_trsp FINE','', (/ ('', i=1, 9) /)) 334 type(ctrl_out),save :: o_d_tr_trsp03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 335 'd_tr_trsp03','cvltr d_tr_trsp COSS','', (/ ('', i=1, 9) /)) 336 type(ctrl_out),save :: o_d_tr_trsp04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 337 'd_tr_trsp04','cvltr d_tr_trsp CODU','', (/ ('', i=1, 9) /)) 338 type(ctrl_out),save :: o_d_tr_trsp05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 339 'd_tr_trsp05','cvltr d_tr_trsp SCDU','', (/ ('', i=1, 9) /)) 340 341 342 type(ctrl_out),save :: o_d_tr_sscav01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 343 'd_tr_sscav01','cvltr d_tr_sscav PREC','', (/ ('', i=1, 9) /)) 344 type(ctrl_out),save :: o_d_tr_sscav02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 345 'd_tr_sscav02','cvltr d_tr_sscav FINE','', (/ ('', i=1, 9) /)) 346 type(ctrl_out),save :: o_d_tr_sscav03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 347 'd_tr_sscav03','cvltr d_tr_sscav COSS','', (/ ('', i=1, 9) /)) 348 type(ctrl_out),save :: o_d_tr_sscav04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 349 'd_tr_sscav04','cvltr d_tr_sscav CODU','', (/ ('', i=1, 9) /)) 350 type(ctrl_out),save :: o_d_tr_sscav05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 351 'd_tr_sscav05','cvltr d_tr_sscav SCDU','', (/ ('', i=1, 9) /)) 352 353 type(ctrl_out),save :: o_d_tr_sat01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 354 'd_tr_sat01','cvltr d_tr_sat PREC','', (/ ('', i=1, 9) /)) 355 type(ctrl_out),save :: o_d_tr_sat02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 356 'd_tr_sat02','cvltr d_tr_sat FINE','', (/ ('', i=1, 9) /)) 357 type(ctrl_out),save :: o_d_tr_sat03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 358 'd_tr_sat03','cvltr d_tr_sat COSS','', (/ ('', i=1, 9) /)) 359 type(ctrl_out),save :: o_d_tr_sat04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 360 'd_tr_sat04','cvltr d_tr_sat CODU','', (/ ('', i=1, 9) /)) 361 type(ctrl_out),save :: o_d_tr_sat05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 362 'd_tr_sat05','cvltr d_tr_sat SCDU','', (/ ('', i=1, 9) /)) 363 364 type(ctrl_out),save :: o_d_tr_uscav01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 365 'd_tr_uscav01','cvltr d_tr_uscav PREC','', (/ ('', i=1, 9) /)) 366 type(ctrl_out),save :: o_d_tr_uscav02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 367 'd_tr_uscav02','cvltr d_tr_uscav FINE','', (/ ('', i=1, 9) /)) 368 type(ctrl_out),save :: o_d_tr_uscav03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 369 'd_tr_uscav03','cvltr d_tr_uscav COSS','', (/ ('', i=1, 9) /)) 370 type(ctrl_out),save :: o_d_tr_uscav04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 371 'd_tr_uscav04','cvltr d_tr_uscav CODU','', (/ ('', i=1, 9) /)) 372 type(ctrl_out),save :: o_d_tr_uscav05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 373 'd_tr_uscav05','cvltr d_tr_uscav SCDU','', (/ ('', i=1, 9) /)) 298 type(ctrl_out),save :: o_dhkelsc01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 299 'dhkelsc01','KE dep hum large scale PREC','', (/ ('', i=1, 10) /)) 300 type(ctrl_out),save :: o_dhkelsc02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 301 'dhkelsc02','KE dep hum large scale FINE','', (/ ('', i=1, 10) /)) 302 type(ctrl_out),save :: o_dhkelsc03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 303 'dhkelsc03','KE dep hum large scale COSS','', (/ ('', i=1, 10) /)) 304 type(ctrl_out),save :: o_dhkelsc04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 305 'dhkelsc04','KE dep hum large scale CODU','', (/ ('', i=1, 10) /)) 306 type(ctrl_out),save :: o_dhkelsc05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 307 'dhkelsc05','KE dep hum large scale SCDU','', (/ ('', i=1, 10) /)) 308 309 type(ctrl_out),save :: o_d_tr_cv01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 310 'd_tr_cv01','cvltr d_tr_cv PREC','', (/ ('', i=1, 10) /)) 311 type(ctrl_out),save :: o_d_tr_cv02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 312 'd_tr_cv02','cvltr d_tr_cv FINE','', (/ ('', i=1, 10) /)) 313 type(ctrl_out),save :: o_d_tr_cv03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 314 'd_tr_cv03','cvltr d_tr_cv COSS','', (/ ('', i=1, 10) /)) 315 type(ctrl_out),save :: o_d_tr_cv04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 316 'd_tr_cv04','cvltr d_tr_cv CODU','', (/ ('', i=1, 10) /)) 317 type(ctrl_out),save :: o_d_tr_cv05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 318 'd_tr_cv05','cvltr d_tr_cv SCDU','', (/ ('', i=1, 10) /)) 319 320 type(ctrl_out),save :: o_d_tr_trsp01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 321 'd_tr_trsp01','cvltr d_tr_trsp PREC','', (/ ('', i=1, 10) /)) 322 type(ctrl_out),save :: o_d_tr_trsp02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 323 'd_tr_trsp02','cvltr d_tr_trsp FINE','', (/ ('', i=1, 10) /)) 324 type(ctrl_out),save :: o_d_tr_trsp03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 325 'd_tr_trsp03','cvltr d_tr_trsp COSS','', (/ ('', i=1, 10) /)) 326 type(ctrl_out),save :: o_d_tr_trsp04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 327 'd_tr_trsp04','cvltr d_tr_trsp CODU','', (/ ('', i=1, 10) /)) 328 type(ctrl_out),save :: o_d_tr_trsp05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 329 'd_tr_trsp05','cvltr d_tr_trsp SCDU','', (/ ('', i=1, 10) /)) 330 331 type(ctrl_out),save :: o_d_tr_sscav01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 332 'd_tr_sscav01','cvltr d_tr_sscav PREC','', (/ ('', i=1, 10) /)) 333 type(ctrl_out),save :: o_d_tr_sscav02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 334 'd_tr_sscav02','cvltr d_tr_sscav FINE','', (/ ('', i=1, 10) /)) 335 type(ctrl_out),save :: o_d_tr_sscav03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 336 'd_tr_sscav03','cvltr d_tr_sscav COSS','', (/ ('', i=1, 10) /)) 337 type(ctrl_out),save :: o_d_tr_sscav04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 338 'd_tr_sscav04','cvltr d_tr_sscav CODU','', (/ ('', i=1, 10) /)) 339 type(ctrl_out),save :: o_d_tr_sscav05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 340 'd_tr_sscav05','cvltr d_tr_sscav SCDU','', (/ ('', i=1, 10) /)) 341 342 type(ctrl_out),save :: o_d_tr_sat01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 343 'd_tr_sat01','cvltr d_tr_sat PREC','', (/ ('', i=1, 10) /)) 344 type(ctrl_out),save :: o_d_tr_sat02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 345 'd_tr_sat02','cvltr d_tr_sat FINE','', (/ ('', i=1, 10) /)) 346 type(ctrl_out),save :: o_d_tr_sat03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 347 'd_tr_sat03','cvltr d_tr_sat COSS','', (/ ('', i=1, 10) /)) 348 type(ctrl_out),save :: o_d_tr_sat04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 349 'd_tr_sat04','cvltr d_tr_sat CODU','', (/ ('', i=1, 10) /)) 350 type(ctrl_out),save :: o_d_tr_sat05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 351 'd_tr_sat05','cvltr d_tr_sat SCDU','', (/ ('', i=1, 10) /)) 352 353 type(ctrl_out),save :: o_d_tr_uscav01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 354 'd_tr_uscav01','cvltr d_tr_uscav PREC','', (/ ('', i=1, 10) /)) 355 type(ctrl_out),save :: o_d_tr_uscav02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 356 'd_tr_uscav02','cvltr d_tr_uscav FINE','', (/ ('', i=1, 10) /)) 357 type(ctrl_out),save :: o_d_tr_uscav03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 358 'd_tr_uscav03','cvltr d_tr_uscav COSS','', (/ ('', i=1, 10) /)) 359 type(ctrl_out),save :: o_d_tr_uscav04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 360 'd_tr_uscav04','cvltr d_tr_uscav CODU','', (/ ('', i=1, 10) /)) 361 type(ctrl_out),save :: o_d_tr_uscav05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 362 'd_tr_uscav05','cvltr d_tr_uscav SCDU','', (/ ('', i=1, 10) /)) 374 363 375 364 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 365 type(ctrl_out),save :: o_d_tr_insc01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 366 'd_tr_insc01','large-scale d_tr_insc PREC','', (/ ('', i=1, 10) /)) 367 type(ctrl_out),save :: o_d_tr_insc02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 368 'd_tr_insc02','large-scale d_tr_insc FINE','', (/ ('', i=1, 10) /)) 369 type(ctrl_out),save :: o_d_tr_insc03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 370 'd_tr_insc03','large-scale d_tr_insc COSS','', (/ ('', i=1, 10) /)) 371 type(ctrl_out),save :: o_d_tr_insc04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 372 'd_tr_insc04','large-scale d_tr_insc CODU','', (/ ('', i=1, 10) /)) 373 type(ctrl_out),save :: o_d_tr_insc05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 374 'd_tr_insc05','large-scale d_tr_insc SCDU','', (/ ('', i=1, 10) /)) 376 375 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 377 type(ctrl_out),save :: o_d_tr_ insc01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &378 'd_tr_ insc01','large-scale d_tr_insc PREC','', (/ ('', i=1, 9) /))379 type(ctrl_out),save :: o_d_tr_ insc02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &380 'd_tr_ insc02','large-scale d_tr_insc FINE','', (/ ('', i=1, 9) /))381 type(ctrl_out),save :: o_d_tr_ insc03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &382 'd_tr_ insc03','large-scale d_tr_insc COSS','', (/ ('', i=1, 9) /))383 type(ctrl_out),save :: o_d_tr_ insc04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &384 'd_tr_ insc04','large-scale d_tr_insc CODU','', (/ ('', i=1, 9) /))385 type(ctrl_out),save :: o_d_tr_ insc05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &386 'd_tr_ insc05','large-scale d_tr_insc SCDU','', (/ ('', i=1, 9) /))376 type(ctrl_out),save :: o_d_tr_bcscav01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 377 'd_tr_bcscav01','large-scale d_tr_bcscav PREC','', (/ ('', i=1, 10) /)) 378 type(ctrl_out),save :: o_d_tr_bcscav02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 379 'd_tr_bcscav02','large-scale d_tr_bcscav FINE','', (/ ('', i=1, 10) /)) 380 type(ctrl_out),save :: o_d_tr_bcscav03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 381 'd_tr_bcscav03','large-scale d_tr_bcscav COSS','', (/ ('', i=1, 10) /)) 382 type(ctrl_out),save :: o_d_tr_bcscav04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 383 'd_tr_bcscav04','large-scale d_tr_bcscav CODU','', (/ ('', i=1, 10) /)) 384 type(ctrl_out),save :: o_d_tr_bcscav05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 385 'd_tr_bcscav05','large-scale d_tr_bcscav SCDU','', (/ ('', i=1, 10) /)) 387 386 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 388 type(ctrl_out),save :: o_d_tr_ bcscav01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &389 'd_tr_ bcscav01','large-scale d_tr_bcscav PREC','', (/ ('', i=1, 9) /))390 type(ctrl_out),save :: o_d_tr_ bcscav02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &391 'd_tr_ bcscav02','large-scale d_tr_bcscav FINE','', (/ ('', i=1, 9) /))392 type(ctrl_out),save :: o_d_tr_ bcscav03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &393 'd_tr_ bcscav03','large-scale d_tr_bcscav COSS','', (/ ('', i=1, 9) /))394 type(ctrl_out),save :: o_d_tr_ bcscav04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &395 'd_tr_ bcscav04','large-scale d_tr_bcscav CODU','', (/ ('', i=1, 9) /))396 type(ctrl_out),save :: o_d_tr_ bcscav05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &397 'd_tr_ bcscav05','large-scale d_tr_bcscav SCDU','', (/ ('', i=1, 9) /))387 type(ctrl_out),save :: o_d_tr_evapls01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 388 'd_tr_evapls01','large-scale d_tr_evapls PREC','', (/ ('', i=1, 10) /)) 389 type(ctrl_out),save :: o_d_tr_evapls02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 390 'd_tr_evapls02','large-scale d_tr_evapls FINE','', (/ ('', i=1, 10) /)) 391 type(ctrl_out),save :: o_d_tr_evapls03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 392 'd_tr_evapls03','large-scale d_tr_evapls COSS','', (/ ('', i=1, 10) /)) 393 type(ctrl_out),save :: o_d_tr_evapls04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 394 'd_tr_evapls04','large-scale d_tr_evapls CODU','', (/ ('', i=1, 10) /)) 395 type(ctrl_out),save :: o_d_tr_evapls05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 396 'd_tr_evapls05','large-scale d_tr_evapls SCDU','', (/ ('', i=1, 10) /)) 398 397 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 399 type(ctrl_out),save :: o_d_tr_ evapls01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &400 'd_tr_ evapls01','large-scale d_tr_evapls PREC','', (/ ('', i=1, 9) /))401 type(ctrl_out),save :: o_d_tr_ evapls02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &402 'd_tr_ evapls02','large-scale d_tr_evapls FINE','', (/ ('', i=1, 9) /))403 type(ctrl_out),save :: o_d_tr_ evapls03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &404 'd_tr_ evapls03','large-scale d_tr_evapls COSS','', (/ ('', i=1, 9) /))405 type(ctrl_out),save :: o_d_tr_ evapls04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &406 'd_tr_ evapls04','large-scale d_tr_evapls CODU','', (/ ('', i=1, 9) /))407 type(ctrl_out),save :: o_d_tr_ evapls05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &408 'd_tr_ evapls05','large-scale d_tr_evapls SCDU','', (/ ('', i=1, 9) /))398 type(ctrl_out),save :: o_d_tr_ls01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 399 'd_tr_ls01','large-scale d_tr_ls PREC','', (/ ('', i=1, 10) /)) 400 type(ctrl_out),save :: o_d_tr_ls02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 401 'd_tr_ls02','large-scale d_tr_ls FINE','', (/ ('', i=1, 10) /)) 402 type(ctrl_out),save :: o_d_tr_ls03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 403 'd_tr_ls03','large-scale d_tr_ls COSS','', (/ ('', i=1, 10) /)) 404 type(ctrl_out),save :: o_d_tr_ls04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 405 'd_tr_ls04','large-scale d_tr_ls CODU','', (/ ('', i=1, 10) /)) 406 type(ctrl_out),save :: o_d_tr_ls05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 407 'd_tr_ls05','large-scale d_tr_ls SCDU','', (/ ('', i=1, 10) /)) 409 408 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 410 type(ctrl_out),save :: o_d_tr_ ls01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &411 'd_tr_ ls01','large-scale d_tr_ls PREC','', (/ ('', i=1, 9) /))412 type(ctrl_out),save :: o_d_tr_ ls02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &413 'd_tr_ ls02','large-scale d_tr_ls FINE','', (/ ('', i=1, 9) /))414 type(ctrl_out),save :: o_d_tr_ ls03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &415 'd_tr_ ls03','large-scale d_tr_ls COSS','', (/ ('', i=1, 9) /))416 type(ctrl_out),save :: o_d_tr_ ls04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &417 'd_tr_ ls04','large-scale d_tr_ls CODU','', (/ ('', i=1, 9) /))418 type(ctrl_out),save :: o_d_tr_ ls05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &419 'd_tr_ ls05','large-scale d_tr_ls SCDU','', (/ ('', i=1, 9) /))409 type(ctrl_out),save :: o_d_tr_dyn01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 410 'd_tr_dyn01','cvltr d_tr_dyn PREC','', (/ ('', i=1, 10) /)) 411 type(ctrl_out),save :: o_d_tr_dyn02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 412 'd_tr_dyn02','cvltr d_tr_dyn FINE','', (/ ('', i=1, 10) /)) 413 type(ctrl_out),save :: o_d_tr_dyn03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 414 'd_tr_dyn03','cvltr d_tr_dyn COSS','', (/ ('', i=1, 10) /)) 415 type(ctrl_out),save :: o_d_tr_dyn04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 416 'd_tr_dyn04','cvltr d_tr_dyn CODU','', (/ ('', i=1, 10) /)) 417 type(ctrl_out),save :: o_d_tr_dyn05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 418 'd_tr_dyn05','cvltr d_tr_dyn SCDU','', (/ ('', i=1, 10) /)) 420 419 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 421 422 420 type(ctrl_out),save :: o_d_tr_cl01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 421 'd_tr_cl01','cvltr d_tr_cl PREC','', (/ ('', i=1, 10) /)) 422 type(ctrl_out),save :: o_d_tr_cl02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 423 'd_tr_cl02','cvltr d_tr_cl FINE','', (/ ('', i=1, 10) /)) 424 type(ctrl_out),save :: o_d_tr_cl03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 425 'd_tr_cl03','cvltr d_tr_cl COSS','', (/ ('', i=1, 10) /)) 426 type(ctrl_out),save :: o_d_tr_cl04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 427 'd_tr_cl04','cvltr d_tr_cl CODU','', (/ ('', i=1, 10) /)) 428 type(ctrl_out),save :: o_d_tr_cl05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 429 'd_tr_cl05','cvltr d_tr_cl SCDU','', (/ ('', i=1, 10) /)) 423 430 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 424 type(ctrl_out),save :: o_d_tr_dyn01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 425 'd_tr_dyn01','cvltr d_tr_dyn PREC','', (/ ('', i=1, 9) /)) 426 type(ctrl_out),save :: o_d_tr_dyn02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 427 'd_tr_dyn02','cvltr d_tr_dyn FINE','', (/ ('', i=1, 9) /)) 428 type(ctrl_out),save :: o_d_tr_dyn03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 429 'd_tr_dyn03','cvltr d_tr_dyn COSS','', (/ ('', i=1, 9) /)) 430 type(ctrl_out),save :: o_d_tr_dyn04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 431 'd_tr_dyn04','cvltr d_tr_dyn CODU','', (/ ('', i=1, 9) /)) 432 type(ctrl_out),save :: o_d_tr_dyn05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 433 'd_tr_dyn05','cvltr d_tr_dyn SCDU','', (/ ('', i=1, 9) /)) 434 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 435 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 436 type(ctrl_out),save :: o_d_tr_cl01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 437 'd_tr_cl01','cvltr d_tr_cl PREC','', (/ ('', i=1, 9) /)) 438 type(ctrl_out),save :: o_d_tr_cl02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 439 'd_tr_cl02','cvltr d_tr_cl FINE','', (/ ('', i=1, 9) /)) 440 type(ctrl_out),save :: o_d_tr_cl03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 441 'd_tr_cl03','cvltr d_tr_cl COSS','', (/ ('', i=1, 9) /)) 442 type(ctrl_out),save :: o_d_tr_cl04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 443 'd_tr_cl04','cvltr d_tr_cl CODU','', (/ ('', i=1, 9) /)) 444 type(ctrl_out),save :: o_d_tr_cl05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 445 'd_tr_cl05','cvltr d_tr_cl SCDU','', (/ ('', i=1, 9) /)) 446 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 447 type(ctrl_out),save :: o_d_tr_th01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 448 'd_tr_th01','cvltr d_tr_th PREC','', (/ ('', i=1, 9) /)) 449 type(ctrl_out),save :: o_d_tr_th02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 450 'd_tr_th02','cvltr d_tr_th FINE','', (/ ('', i=1, 9) /)) 451 type(ctrl_out),save :: o_d_tr_th03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 452 'd_tr_th03','cvltr d_tr_th COSS','', (/ ('', i=1, 9) /)) 453 type(ctrl_out),save :: o_d_tr_th04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 454 'd_tr_th04','cvltr d_tr_th CODU','', (/ ('', i=1, 9) /)) 455 type(ctrl_out),save :: o_d_tr_th05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 456 'd_tr_th05','cvltr d_tr_th SCDU','', (/ ('', i=1, 9) /)) 457 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 458 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 431 type(ctrl_out),save :: o_d_tr_th01 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 432 'd_tr_th01','cvltr d_tr_th PREC','', (/ ('', i=1, 10) /)) 433 type(ctrl_out),save :: o_d_tr_th02 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 434 'd_tr_th02','cvltr d_tr_th FINE','', (/ ('', i=1, 10) /)) 435 type(ctrl_out),save :: o_d_tr_th03 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 436 'd_tr_th03','cvltr d_tr_th COSS','', (/ ('', i=1, 10) /)) 437 type(ctrl_out),save :: o_d_tr_th04 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 438 'd_tr_th04','cvltr d_tr_th CODU','', (/ ('', i=1, 10) /)) 439 type(ctrl_out),save :: o_d_tr_th05 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 440 'd_tr_th05','cvltr d_tr_th SCDU','', (/ ('', i=1, 10) /)) 459 441 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 460 type(ctrl_out),save :: o_sed_ss3D = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &461 'sed_ss3D','Tendance Sedmet. COSS','', (/ ('', i=1, 9) /))462 type(ctrl_out),save :: o_sed_dust3D = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &463 'sed_dust3D','Tendance Sedmet. CODU','', (/ ('', i=1, 9) /))464 type(ctrl_out),save :: o_sed_dustsco3D = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &465 'sed_dustsco3D','Tendance Sedmet. SCDU','', (/ ('', i=1, 9) /))442 type(ctrl_out),save :: o_sed_ss3D = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 443 'sed_ss3D','Tendance Sedmet. COSS','', (/ ('', i=1, 10) /)) 444 type(ctrl_out),save :: o_sed_dust3D = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 445 'sed_dust3D','Tendance Sedmet. CODU','', (/ ('', i=1, 10) /)) 446 type(ctrl_out),save :: o_sed_dustsco3D = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 447 'sed_dustsco3D','Tendance Sedmet. SCDU','', (/ ('', i=1, 10) /)) 466 448 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 467 468 type(ctrl_out),save :: o_sed_ss = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 469 'sed_ss','Sedmet. COSS','', (/ ('', i=1, 9) /)) 470 type(ctrl_out),save :: o_sed_dust = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 471 'sed_dust','Sedmet. CODU','', (/ ('', i=1, 9) /)) 472 type(ctrl_out),save :: o_sed_dustsco = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 473 'sed_dustsco','Sedmet. SCDU','', (/ ('', i=1, 9) /)) 474 type(ctrl_out),save :: o_g2p_gas = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 475 'g2p_gas','Gas2particle gas sink','', (/ ('', i=1, 9) /)) 476 type(ctrl_out),save :: o_g2p_aer = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 477 'g2p_aer','Gas2particle tr2 src','', (/ ('', i=1, 9) /)) 478 449 type(ctrl_out),save :: o_sed_ss = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 450 'sed_ss','Sedmet. COSS','', (/ ('', i=1, 10) /)) 451 type(ctrl_out),save :: o_sed_dust = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 452 'sed_dust','Sedmet. CODU','', (/ ('', i=1, 10) /)) 453 type(ctrl_out),save :: o_sed_dustsco = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 454 'sed_dustsco','Sedmet. SCDU','', (/ ('', i=1, 10) /)) 455 type(ctrl_out),save :: o_g2p_gas = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 456 'g2p_gas','Gas2particle gas sink','', (/ ('', i=1, 10) /)) 457 type(ctrl_out),save :: o_g2p_aer = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 458 'g2p_aer','Gas2particle tr2 src','', (/ ('', i=1, 10) /)) 479 459 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 480 460 ! histrac 481 461 482 type(ctrl_out),save :: o_fluxbb = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 483 'fluxbb','Flux BB','mg/m2/s', (/ ('', i=1, 9) /)) 484 485 type(ctrl_out),save :: o_fluxff = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 486 'fluxff','Flux FF','mg/m2/s', (/ ('', i=1, 9) /)) 487 488 type(ctrl_out),save :: o_fluxbcbb = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 489 'fluxbcbb','Flux BC-BB','mg/m2/s', (/ ('', i=1, 9) /)) 490 491 type(ctrl_out),save :: o_fluxbcff = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 492 'fluxbcff','Flux BC-FF','mg/m2/s', (/ ('', i=1, 9) /)) 493 494 type(ctrl_out),save :: o_fluxbcnff = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 495 'fluxbcnff','Flux BC-NFF','mg/m2/s', (/ ('', i=1, 9) /)) 496 497 type(ctrl_out),save :: o_fluxbcba = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 498 'fluxbcba','Flux BC-BA','mg/m2/s', (/ ('', i=1, 9) /)) 499 500 type(ctrl_out),save :: o_fluxbc = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 501 'fluxbc','Flux BC','mg/m2/s', (/ ('', i=1, 9) /)) 502 503 type(ctrl_out),save :: o_fluxombb = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 504 'fluxombb','Flux OM-BB','mg/m2/s', (/ ('', i=1, 9) /)) 505 506 type(ctrl_out),save :: o_fluxomff = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 507 'fluxomff','Flux OM-FF','mg/m2/s', (/ ('', i=1, 9) /)) 508 509 type(ctrl_out),save :: o_fluxomnff = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 510 'fluxomnff','Flux OM-NFF','mg/m2/s', (/ ('', i=1, 9) /)) 511 512 type(ctrl_out),save :: o_fluxomba = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 513 'fluxomba','Flux OM-BA','mg/m2/s', (/ ('', i=1, 9) /)) 514 515 type(ctrl_out),save :: o_fluxomnat = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 516 'fluxomnat','Flux OM-NT','mg/m2/s', (/ ('', i=1, 9) /)) 517 518 type(ctrl_out),save :: o_fluxom = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 519 'fluxom','Flux OM','mg/m2/s', (/ ('', i=1, 9) /)) 520 521 type(ctrl_out),save :: o_fluxh2sff = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 522 'fluxh2sff','Flux H2S FF','mgS/m2/s', (/ ('', i=1, 9) /)) 523 524 type(ctrl_out),save :: o_fluxh2snff = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 525 'fluxh2snff','Flux H2S non-FF','mgS/m2/s', (/ ('', i=1, 9) /)) 526 527 type(ctrl_out),save :: o_fluxso2ff = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 528 'fluxso2ff','Flux SO2 FF','mgS/m2/s', (/ ('', i=1, 9) /)) 529 530 type(ctrl_out),save :: o_fluxso2nff = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 531 'fluxso2nff','Flux SO2 non-FF','mgS/m2/s', (/ ('', i=1, 9) /)) 532 533 type(ctrl_out),save :: o_fluxso2bb = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 534 'fluxso2bb','Flux SO2 BB','mgS/m2/s', (/ ('', i=1, 9) /)) 535 536 type(ctrl_out),save :: o_fluxso2vol = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 537 'fluxso2vol','Flux SO2 Vol','mgS/m2/s', (/ ('', i=1, 9) /)) 538 539 type(ctrl_out),save :: o_fluxso2ba = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 540 'fluxso2ba','Flux SO2 Ba','mgS/m2/s', (/ ('', i=1, 9) /)) 541 542 type(ctrl_out),save :: o_fluxso2 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 543 'fluxso2','Flux SO2','mgS/m2/s', (/ ('', i=1, 9) /)) 544 545 type(ctrl_out),save :: o_fluxso4ff = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 546 'fluxso4ff','Flux SO4 FF','mgS/m2/s', (/ ('', i=1, 9) /)) 547 548 type(ctrl_out),save :: o_fluxso4nff = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 549 'fluxso4nff','Flux SO4 non-FF','mgS/m2/s', (/ ('', i=1, 9) /)) 550 551 type(ctrl_out),save :: o_fluxso4bb = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 552 'fluxso4bb','Flux SO4 BB','mgS/m2/s', (/ ('', i=1, 9) /)) 553 554 type(ctrl_out),save :: o_fluxso4ba = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 555 'fluxso4ba','Flux SO4 Ba','mgS/m2/s', (/ ('', i=1, 9) /)) 556 557 type(ctrl_out),save :: o_fluxso4 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 558 'fluxso4','Flux SO4','mgS/m2/s', (/ ('', i=1, 9) /)) 559 560 type(ctrl_out),save :: o_fluxdms = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 561 'fluxdms','Flux DMS','mgS/m2/s', (/ ('', i=1, 9) /)) 562 563 type(ctrl_out),save :: o_fluxh2sbio = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 564 'fluxh2sbio','Flux H2S Bio','mgS/m2/s', (/ ('', i=1, 9) /)) 565 566 type(ctrl_out),save :: o_fluxdustec = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 567 'fluxdustec','Flux Dust EC','mg/m2/s', (/ ('', i=1, 9) /)) 568 569 type(ctrl_out),save :: o_fluxddfine = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 570 'fluxddfine','DD Fine Mode','mg/m2/s', (/ ('', i=1, 9) /)) 571 572 type(ctrl_out),save :: o_fluxddcoa = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 573 'fluxddcoa','DD Coarse Mode','mg/m2/s', (/ ('', i=1, 9) /)) 574 575 type(ctrl_out),save :: o_fluxddsco = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 576 'fluxddsco','DD SCoarse Mode','mg/m2/s', (/ ('', i=1, 9) /)) 577 578 type(ctrl_out),save :: o_fluxdd = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 579 'fluxdd','Flux DD','mg/m2/s', (/ ('', i=1, 9) /)) 580 581 type(ctrl_out),save :: o_fluxssfine = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 582 'fluxssfine','SS Fine Mode','mg/m2/s', (/ ('', i=1, 9) /)) 583 584 type(ctrl_out),save :: o_fluxsscoa = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 585 'fluxsscoa','SS Coarse Mode','mg/m2/s', (/ ('', i=1, 9) /)) 586 587 type(ctrl_out),save :: o_fluxss = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 588 'fluxss','Flux SS','mg/m2/s', (/ ('', i=1, 9) /)) 589 590 type(ctrl_out),save :: o_flux_sparam_ind = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 591 'flux_sparam_ind','Ind emiss','mg/m2/s', (/ ('', i=1, 9) /)) 592 593 type(ctrl_out),save :: o_flux_sparam_bb = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 594 'flux_sparam_bb','BB emiss','mg/m2/s', (/ ('', i=1, 9) /)) 595 596 type(ctrl_out),save :: o_flux_sparam_ff = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 597 'flux_sparam_ff','FF emiss','mg/m2/s', (/ ('', i=1, 9) /)) 598 599 type(ctrl_out),save :: o_flux_sparam_ddfine = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 600 'flux_sparam_ddfine','DD fine emiss','mg/m2/s', (/ ('', i=1, 9) /)) 601 602 type(ctrl_out),save :: o_flux_sparam_ddcoa = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 603 'flux_sparam_ddcoa','DD coarse emiss','mg/m2/s', (/ ('', i=1, 9) /)) 604 605 type(ctrl_out),save :: o_flux_sparam_ddsco = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 606 'flux_sparam_ddsco','DD Scoarse emiss','mg/m2/s', (/ ('', i=1, 9) /)) 607 608 type(ctrl_out),save :: o_flux_sparam_ssfine = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 609 'flux_sparam_ssfine','SS fine emiss','mg/m2/s', (/ ('', i=1, 9) /)) 610 611 612 type(ctrl_out),save :: o_flux_sparam_sscoa = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 613 'flux_sparam_sscoa','SS coarse emiss','mg/m2/s', (/ ('', i=1, 9) /)) 614 615 type(ctrl_out),save :: o_u10m_ss = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 616 'u10m_ss','Zonal wind at 10 m SS','m/s', (/ ('', i=1, 9) /)) 617 618 type(ctrl_out),save :: o_v10m_ss = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 619 'v10m_ss','Meridional wind at 10 m SS','m/s', (/ ('', i=1, 9) /)) 620 621 622 ! type(ctrl_out),save :: o_ = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), & 623 ! '','','', (/ ('', i=1, 9) /)) 624 625 626 627 !example TYPE(ctrl_out), SAVE :: o_psbg = ctrl_out((/ 11, 11, 11, 11, 11, 11, 5, 5, 5 /), & 462 type(ctrl_out),save :: o_fluxbb = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 463 'fluxbb','Flux BB','mg/m2/s', (/ ('', i=1, 10) /)) 464 465 type(ctrl_out),save :: o_fluxff = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 466 'fluxff','Flux FF','mg/m2/s', (/ ('', i=1, 10) /)) 467 468 type(ctrl_out),save :: o_fluxbcbb = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 469 'fluxbcbb','Flux BC-BB','mg/m2/s', (/ ('', i=1, 10) /)) 470 471 type(ctrl_out),save :: o_fluxbcff = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 472 'fluxbcff','Flux BC-FF','mg/m2/s', (/ ('', i=1, 10) /)) 473 474 type(ctrl_out),save :: o_fluxbcnff = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 475 'fluxbcnff','Flux BC-NFF','mg/m2/s', (/ ('', i=1, 10) /)) 476 477 type(ctrl_out),save :: o_fluxbcba = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 478 'fluxbcba','Flux BC-BA','mg/m2/s', (/ ('', i=1, 10) /)) 479 480 type(ctrl_out),save :: o_fluxbc = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 481 'fluxbc','Flux BC','mg/m2/s', (/ ('', i=1, 10) /)) 482 483 type(ctrl_out),save :: o_fluxombb = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 484 'fluxombb','Flux OM-BB','mg/m2/s', (/ ('', i=1, 10) /)) 485 486 type(ctrl_out),save :: o_fluxomff = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 487 'fluxomff','Flux OM-FF','mg/m2/s', (/ ('', i=1, 10) /)) 488 489 type(ctrl_out),save :: o_fluxomnff = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 490 'fluxomnff','Flux OM-NFF','mg/m2/s', (/ ('', i=1, 10) /)) 491 492 type(ctrl_out),save :: o_fluxomba = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 493 'fluxomba','Flux OM-BA','mg/m2/s', (/ ('', i=1, 10) /)) 494 495 type(ctrl_out),save :: o_fluxomnat = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 496 'fluxomnat','Flux OM-NT','mg/m2/s', (/ ('', i=1, 10) /)) 497 498 type(ctrl_out),save :: o_fluxom = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 499 'fluxom','Flux OM','mg/m2/s', (/ ('', i=1, 10) /)) 500 501 type(ctrl_out),save :: o_fluxh2sff = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 502 'fluxh2sff','Flux H2S FF','mgS/m2/s', (/ ('', i=1, 10) /)) 503 504 type(ctrl_out),save :: o_fluxh2snff = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 505 'fluxh2snff','Flux H2S non-FF','mgS/m2/s', (/ ('', i=1, 10) /)) 506 507 type(ctrl_out),save :: o_fluxso2ff = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 508 'fluxso2ff','Flux SO2 FF','mgS/m2/s', (/ ('', i=1, 10) /)) 509 510 type(ctrl_out),save :: o_fluxso2nff = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 511 'fluxso2nff','Flux SO2 non-FF','mgS/m2/s', (/ ('', i=1, 10) /)) 512 513 type(ctrl_out),save :: o_fluxso2bb = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 514 'fluxso2bb','Flux SO2 BB','mgS/m2/s', (/ ('', i=1, 10) /)) 515 516 type(ctrl_out),save :: o_fluxso2vol = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 517 'fluxso2vol','Flux SO2 Vol','mgS/m2/s', (/ ('', i=1, 10) /)) 518 519 type(ctrl_out),save :: o_fluxso2ba = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 520 'fluxso2ba','Flux SO2 Ba','mgS/m2/s', (/ ('', i=1, 10) /)) 521 522 type(ctrl_out),save :: o_fluxso2 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 523 'fluxso2','Flux SO2','mgS/m2/s', (/ ('', i=1, 10) /)) 524 525 type(ctrl_out),save :: o_fluxso4ff = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 526 'fluxso4ff','Flux SO4 FF','mgS/m2/s', (/ ('', i=1, 10) /)) 527 528 type(ctrl_out),save :: o_fluxso4nff = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 529 'fluxso4nff','Flux SO4 non-FF','mgS/m2/s', (/ ('', i=1, 10) /)) 530 531 type(ctrl_out),save :: o_fluxso4bb = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 532 'fluxso4bb','Flux SO4 BB','mgS/m2/s', (/ ('', i=1, 10) /)) 533 534 type(ctrl_out),save :: o_fluxso4ba = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 535 'fluxso4ba','Flux SO4 Ba','mgS/m2/s', (/ ('', i=1, 10) /)) 536 537 type(ctrl_out),save :: o_fluxso4 = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 538 'fluxso4','Flux SO4','mgS/m2/s', (/ ('', i=1, 10) /)) 539 540 type(ctrl_out),save :: o_fluxdms = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 541 'fluxdms','Flux DMS','mgS/m2/s', (/ ('', i=1, 10) /)) 542 543 type(ctrl_out),save :: o_fluxh2sbio = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 544 'fluxh2sbio','Flux H2S Bio','mgS/m2/s', (/ ('', i=1, 10) /)) 545 546 type(ctrl_out),save :: o_fluxdustec = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 547 'fluxdustec','Flux Dust EC','mg/m2/s', (/ ('', i=1, 10) /)) 548 549 type(ctrl_out),save :: o_fluxddfine = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 550 'fluxddfine','DD Fine Mode','mg/m2/s', (/ ('', i=1, 10) /)) 551 552 type(ctrl_out),save :: o_fluxddcoa = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 553 'fluxddcoa','DD Coarse Mode','mg/m2/s', (/ ('', i=1, 10) /)) 554 555 type(ctrl_out),save :: o_fluxddsco = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 556 'fluxddsco','DD SCoarse Mode','mg/m2/s', (/ ('', i=1, 10) /)) 557 558 type(ctrl_out),save :: o_fluxdd = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 559 'fluxdd','Flux DD','mg/m2/s', (/ ('', i=1, 10) /)) 560 561 type(ctrl_out),save :: o_fluxssfine = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 562 'fluxssfine','SS Fine Mode','mg/m2/s', (/ ('', i=1, 10) /)) 563 564 type(ctrl_out),save :: o_fluxsscoa = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 565 'fluxsscoa','SS Coarse Mode','mg/m2/s', (/ ('', i=1, 10) /)) 566 567 type(ctrl_out),save :: o_fluxss = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 568 'fluxss','Flux SS','mg/m2/s', (/ ('', i=1, 10) /)) 569 570 type(ctrl_out),save :: o_flux_sparam_ind = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 571 'flux_sparam_ind','Ind emiss','mg/m2/s', (/ ('', i=1, 10) /)) 572 573 type(ctrl_out),save :: o_flux_sparam_bb = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 574 'flux_sparam_bb','BB emiss','mg/m2/s', (/ ('', i=1, 10) /)) 575 576 type(ctrl_out),save :: o_flux_sparam_ff = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 577 'flux_sparam_ff','FF emiss','mg/m2/s', (/ ('', i=1, 10) /)) 578 579 type(ctrl_out),save :: o_flux_sparam_ddfine = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 580 'flux_sparam_ddfine','DD fine emiss','mg/m2/s', (/ ('', i=1, 10) /)) 581 582 type(ctrl_out),save :: o_flux_sparam_ddcoa = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 583 'flux_sparam_ddcoa','DD coarse emiss','mg/m2/s', (/ ('', i=1, 10) /)) 584 585 type(ctrl_out),save :: o_flux_sparam_ddsco = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 586 'flux_sparam_ddsco','DD Scoarse emiss','mg/m2/s', (/ ('', i=1, 10) /)) 587 588 type(ctrl_out),save :: o_flux_sparam_ssfine = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 589 'flux_sparam_ssfine','SS fine emiss','mg/m2/s', (/ ('', i=1, 10) /)) 590 591 type(ctrl_out),save :: o_flux_sparam_sscoa = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 592 'flux_sparam_sscoa','SS coarse emiss','mg/m2/s', (/ ('', i=1, 10) /)) 593 594 type(ctrl_out),save :: o_u10m_ss = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 595 'u10m_ss','Zonal wind at 10 m SS','m/s', (/ ('', i=1, 10) /)) 596 597 type(ctrl_out),save :: o_v10m_ss = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 598 'v10m_ss','Meridional wind at 10 m SS','m/s', (/ ('', i=1, 10) /)) 599 600 ! type(ctrl_out),save :: o_ = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), & 601 ! '','','', (/ ('', i=1, 10) /)) 602 603 !example TYPE(ctrl_out), SAVE :: o_psbg = ctrl_out((/ 11, 11, 11, 11, 11, 11, 5, 5, 5, 11 /), & 628 604 ! 'psbg', 'Pressure sfce below ground', '%', (/ "inst(X)", "inst(X)", "inst(X)", & 629 ! "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /)) 630 631 632 605 ! "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /)) -
LMDZ5/branches/testing/libf/phylmd/Dust/spla_output_write.h
r2641 r2787 6 6 ! traceur 7 7 8 CALL histwrite_phy( o_taue550 , diff_aod550_tot 9 CALL histwrite_phy( o_taue670 , diag_aod670_tot 8 CALL histwrite_phy( o_taue550 , diff_aod550_tot ) 9 CALL histwrite_phy( o_taue670 , diag_aod670_tot ) 10 10 CALL histwrite_phy( o_taue865 , diag_aod865_tot ) 11 11 IF(id_fine>0) CALL histwrite_phy( o_taue550_tr2 , diff_aod550_tr2 ) 12 IF(id_fine>0) CALL histwrite_phy( o_taue670_tr2 , diag_aod670_tr2 12 IF(id_fine>0) CALL histwrite_phy( o_taue670_tr2 , diag_aod670_tr2 ) 13 13 IF(id_fine>0) CALL histwrite_phy( o_taue865_tr2 , diag_aod865_tr2 ) 14 IF(id_coss>0) CALL histwrite_phy( o_taue550_ss , diag_aod550_ss 14 IF(id_coss>0) CALL histwrite_phy( o_taue550_ss , diag_aod550_ss ) 15 15 IF(id_coss>0) CALL histwrite_phy( o_taue670_ss , diag_aod670_ss ) 16 IF(id_coss>0) CALL histwrite_phy( o_taue865_ss , diag_aod865_ss 17 IF(id_codu>0) CALL histwrite_phy( o_taue550_dust , diag_aod550_dust 16 IF(id_coss>0) CALL histwrite_phy( o_taue865_ss , diag_aod865_ss ) 17 IF(id_codu>0) CALL histwrite_phy( o_taue550_dust , diag_aod550_dust ) 18 18 IF(id_codu>0) CALL histwrite_phy( o_taue670_dust , diag_aod670_dust ) 19 19 IF(id_codu>0) CALL histwrite_phy( o_taue865_dust , diag_aod865_dust ) … … 22 22 IF(id_scdu>0) CALL histwrite_phy( o_taue865_dustsco , diag_aod865_dustsco ) 23 23 CALL histwrite_phy( o_taue550_aqua , aod550_aqua ) 24 CALL histwrite_phy( o_taue550_terra , aod550_terra 25 CALL histwrite_phy( o_taue670_aqua , aod670_aqua 26 CALL histwrite_phy( o_taue670_terra , aod670_terra 27 CALL histwrite_phy( o_taue865_aqua , aod865_aqua 28 CALL histwrite_phy( o_taue865_terra , aod865_terra 24 CALL histwrite_phy( o_taue550_terra , aod550_terra ) 25 CALL histwrite_phy( o_taue670_aqua , aod670_aqua ) 26 CALL histwrite_phy( o_taue670_terra , aod670_terra ) 27 CALL histwrite_phy( o_taue865_aqua , aod865_aqua ) 28 CALL histwrite_phy( o_taue865_terra , aod865_terra ) 29 29 30 30 IF(id_fine>0) CALL histwrite_phy( o_taue550_fine_aqua ,aod550_tr2_aqua ) 31 31 IF(id_fine>0) CALL histwrite_phy( o_taue670_fine_aqua ,aod670_tr2_aqua ) 32 32 IF(id_fine>0) CALL histwrite_phy( o_taue865_fine_aqua ,aod865_tr2_aqua ) 33 IF(id_coss>0) CALL histwrite_phy( o_taue550_coss_aqua ,aod550_ss_aqua 34 IF(id_coss>0) CALL histwrite_phy( o_taue670_coss_aqua ,aod670_ss_aqua 35 IF(id_coss>0) CALL histwrite_phy( o_taue865_coss_aqua ,aod865_ss_aqua 36 IF(id_codu>0) CALL histwrite_phy( o_taue550_codu_aqua ,aod550_dust_aqua 37 IF(id_codu>0) CALL histwrite_phy( o_taue670_codu_aqua ,aod670_dust_aqua 38 IF(id_codu>0) CALL histwrite_phy( o_taue865_codu_aqua ,aod865_dust_aqua 39 IF(id_scdu>0) CALL histwrite_phy( o_taue670_scdu_aqua ,aod670_dustsco_aqua 40 IF(id_scdu>0) CALL histwrite_phy( o_taue550_scdu_aqua ,aod550_dustsco_aqua 41 IF(id_scdu>0) CALL histwrite_phy( o_taue865_scdu_aqua ,aod865_dustsco_aqua 33 IF(id_coss>0) CALL histwrite_phy( o_taue550_coss_aqua ,aod550_ss_aqua ) 34 IF(id_coss>0) CALL histwrite_phy( o_taue670_coss_aqua ,aod670_ss_aqua ) 35 IF(id_coss>0) CALL histwrite_phy( o_taue865_coss_aqua ,aod865_ss_aqua ) 36 IF(id_codu>0) CALL histwrite_phy( o_taue550_codu_aqua ,aod550_dust_aqua ) 37 IF(id_codu>0) CALL histwrite_phy( o_taue670_codu_aqua ,aod670_dust_aqua ) 38 IF(id_codu>0) CALL histwrite_phy( o_taue865_codu_aqua ,aod865_dust_aqua ) 39 IF(id_scdu>0) CALL histwrite_phy( o_taue670_scdu_aqua ,aod670_dustsco_aqua ) 40 IF(id_scdu>0) CALL histwrite_phy( o_taue550_scdu_aqua ,aod550_dustsco_aqua ) 41 IF(id_scdu>0) CALL histwrite_phy( o_taue865_scdu_aqua ,aod865_dustsco_aqua ) 42 42 43 43 IF(id_fine>0) CALL histwrite_phy( o_taue550_fine_terra ,aod550_tr2_terra ) 44 44 IF(id_fine>0) CALL histwrite_phy( o_taue670_fine_terra ,aod670_tr2_terra ) 45 45 IF(id_fine>0) CALL histwrite_phy( o_taue865_fine_terra ,aod865_tr2_terra ) 46 IF(id_coss>0) CALL histwrite_phy( o_taue550_coss_terra ,aod550_ss_terra 47 IF(id_coss>0) CALL histwrite_phy( o_taue670_coss_terra ,aod670_ss_terra 48 IF(id_coss>0) CALL histwrite_phy( o_taue865_coss_terra ,aod865_ss_terra 49 IF(id_codu>0) CALL histwrite_phy( o_taue550_codu_terra ,aod550_dust_terra 50 IF(id_codu>0) CALL histwrite_phy( o_taue670_codu_terra ,aod670_dust_terra 51 IF(id_codu>0) CALL histwrite_phy( o_taue865_codu_terra ,aod865_dust_terra 52 IF(id_scdu>0) CALL histwrite_phy( o_taue670_scdu_terra ,aod670_dustsco_terra 53 IF(id_scdu>0) CALL histwrite_phy( o_taue550_scdu_terra ,aod550_dustsco_terra 54 IF(id_scdu>0) CALL histwrite_phy( o_taue865_scdu_terra ,aod865_dustsco_terra 46 IF(id_coss>0) CALL histwrite_phy( o_taue550_coss_terra ,aod550_ss_terra ) 47 IF(id_coss>0) CALL histwrite_phy( o_taue670_coss_terra ,aod670_ss_terra ) 48 IF(id_coss>0) CALL histwrite_phy( o_taue865_coss_terra ,aod865_ss_terra ) 49 IF(id_codu>0) CALL histwrite_phy( o_taue550_codu_terra ,aod550_dust_terra ) 50 IF(id_codu>0) CALL histwrite_phy( o_taue670_codu_terra ,aod670_dust_terra ) 51 IF(id_codu>0) CALL histwrite_phy( o_taue865_codu_terra ,aod865_dust_terra ) 52 IF(id_scdu>0) CALL histwrite_phy( o_taue670_scdu_terra ,aod670_dustsco_terra ) 53 IF(id_scdu>0) CALL histwrite_phy( o_taue550_scdu_terra ,aod550_dustsco_terra ) 54 IF(id_scdu>0) CALL histwrite_phy( o_taue865_scdu_terra ,aod865_dustsco_terra ) 55 55 56 56 … … 133 133 IF(id_codu>0) CALL histwrite_phy( o_d_tr_uscav04 , d_tr_uscav04 ) 134 134 IF(id_scdu>0) CALL histwrite_phy( o_d_tr_uscav05 , d_tr_uscav05 ) 135 IF(id_prec>0) CALL histwrite_phy( o_d_tr_insc01 , d_tr_insc01 ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!135 IF(id_prec>0) CALL histwrite_phy( o_d_tr_insc01 , d_tr_insc01 ) 136 136 IF(id_fine>0) CALL histwrite_phy( o_d_tr_insc02 , d_tr_insc02 ) 137 137 IF(id_coss>0) CALL histwrite_phy( o_d_tr_insc03 , d_tr_insc03 ) … … 152 152 IF(id_coss>0) CALL histwrite_phy( o_d_tr_ls03 , d_tr_ls03 ) 153 153 IF(id_codu>0) CALL histwrite_phy( o_d_tr_ls04 , d_tr_ls04 ) 154 IF(id_scdu>0) CALL histwrite_phy( o_d_tr_ls05 , d_tr_ls05 ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 155 154 IF(id_scdu>0) CALL histwrite_phy( o_d_tr_ls05 , d_tr_ls05 ) 156 155 157 156 IF(id_prec>0) CALL histwrite_phy( o_d_tr_dyn01 , d_tr_dyn01 ) … … 159 158 IF(id_coss>0) CALL histwrite_phy( o_d_tr_dyn03 , d_tr_dyn03 ) 160 159 IF(id_codu>0) CALL histwrite_phy( o_d_tr_dyn04 , d_tr_dyn04 ) 161 IF(id_scdu>0) CALL histwrite_phy( o_d_tr_dyn05 , d_tr_dyn05 ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!160 IF(id_scdu>0) CALL histwrite_phy( o_d_tr_dyn05 , d_tr_dyn05 ) 162 161 163 162 IF(id_prec>0) CALL histwrite_phy( o_d_tr_cl01 , d_tr_cl01 ) … … 165 164 IF(id_coss>0) CALL histwrite_phy( o_d_tr_cl03 , d_tr_cl03 ) 166 165 IF(id_codu>0) CALL histwrite_phy( o_d_tr_cl04 , d_tr_cl04 ) 167 IF(id_scdu>0) CALL histwrite_phy( o_d_tr_cl05 , d_tr_cl05 ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!166 IF(id_scdu>0) CALL histwrite_phy( o_d_tr_cl05 , d_tr_cl05 ) 168 167 IF(id_prec>0) CALL histwrite_phy( o_d_tr_th01 , d_tr_th01 ) 169 168 IF(id_fine>0) CALL histwrite_phy( o_d_tr_th02 , d_tr_th02 ) 170 169 IF(id_coss>0) CALL histwrite_phy( o_d_tr_th03 , d_tr_th03 ) 171 170 IF(id_codu>0) CALL histwrite_phy( o_d_tr_th04 , d_tr_th04 ) 172 IF(id_scdu>0) CALL histwrite_phy( o_d_tr_th05 , d_tr_th05 ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 173 171 IF(id_scdu>0) CALL histwrite_phy( o_d_tr_th05 , d_tr_th05 ) 174 172 175 173 IF(id_coss>0) CALL histwrite_phy( o_sed_ss , sed_ss) … … 232 230 CALL histwrite_phy( o_u10m_ss , u10m_ss ) 233 231 CALL histwrite_phy( o_v10m_ss , v10m_ss ) 234 235 ! CALL histwrite_phy( ,)236 !237 238 239 240
Note: See TracChangeset
for help on using the changeset viewer.
