- Timestamp:
- May 18, 2015, 2:27:29 PM (9 years ago)
- Location:
- LMDZ5/trunk
- Files:
-
- 5 edited
-
DefLists/field_def_lmdz.xml (modified) (2 diffs)
-
DefLists/file_def_histdayNMC_lmdz.xml (modified) (2 diffs)
-
DefLists/file_def_histhfNMC_lmdz.xml (modified) (2 diffs)
-
DefLists/file_def_histmthNMC_lmdz.xml (modified) (2 diffs)
-
libf/phylmd/phys_output_write_mod.F90 (modified) (6 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/DefLists/field_def_lmdz.xml
r2279 r2282 547 547 548 548 <field_group id="fields_NMC" domain_ref="dom_glo" axis_ref="plev"> 549 <field id="tnondef" long_name="Undefined value of T" unit="K" />549 <!-- <field id="tnondef" long_name="Undefined value of T" unit="K" /> --> 550 550 <field id="ta" long_name="Air temperature" unit="K" /> 551 551 <field id="zg" long_name="Geopotential height" unit="m" /> … … 555 555 <field id="va" long_name="Northward wind" unit="m s-1" /> 556 556 <field id="wap" long_name="Lagrangian tendency of air pressure" unit="Pa s-1" /> 557 <field id="psbg" long_name="Pressure sfce below ground" unit="%" />557 <!-- <field id="psbg" long_name="Pressure sfce below ground" unit="%" /> --> 558 558 <field id="tro3" long_name="Ozone mole fraction" unit="1e-9" /> 559 559 <field id="tro3_daylight" long_name="Daylight ozone mole fraction" unit="1e-9" /> -
LMDZ5/trunk/DefLists/file_def_histdayNMC_lmdz.xml
r2122 r2282 5 5 <!-- VARS 3D --> 6 6 <field_group operation="instant" freq_op="1ts" > 7 <field field_ref="tnondef" level="5" />7 <!-- <field field_ref="tnondef" level="5" /> --> 8 8 <field field_ref="ta" level="1" /> 9 9 <field field_ref="zg" level="5" /> … … 13 13 <field field_ref="va" level="5" /> 14 14 <field field_ref="wap" level="5" /> 15 <field field_ref="psbg" level="5" />15 <!-- <field field_ref="psbg" level="5" /> --> 16 16 <field field_ref="tro3" level="5" /> 17 17 <field field_ref="tro3_daylight" level="5" /> -
LMDZ5/trunk/DefLists/file_def_histhfNMC_lmdz.xml
r2122 r2282 5 5 <!-- VARS 3D --> 6 6 <field_group operation="instant" freq_op="1ts" > 7 <field field_ref="tnondef" level="5" />7 <!-- <field field_ref="tnondef" level="5" /> --> 8 8 <field field_ref="ta" level="1" /> 9 9 <field field_ref="zg" level="5" /> … … 13 13 <field field_ref="va" level="5" /> 14 14 <field field_ref="wap" level="5" /> 15 <field field_ref="psbg" level="5" />15 <!-- <field field_ref="psbg" level="5" /> --> 16 16 <field field_ref="tro3" level="5" /> 17 17 <field field_ref="tro3_daylight" level="5" /> -
LMDZ5/trunk/DefLists/file_def_histmthNMC_lmdz.xml
r2122 r2282 5 5 <!-- VARS 3D --> 6 6 <field_group operation="instant" freq_op="1ts" > 7 <field field_ref="tnondef" level="5" />7 <!-- <field field_ref="tnondef" level="5" /> --> 8 8 <field field_ref="ta" level="1" /> 9 9 <field field_ref="zg" level="5" /> … … 13 13 <field field_ref="va" level="5" /> 14 14 <field field_ref="wap" level="5" /> 15 <field field_ref="psbg" level="5" />15 <!-- <field field_ref="psbg" level="5" /> --> 16 16 <field field_ref="tro3" level="5" /> 17 17 <field field_ref="tro3_daylight" level="5" /> -
LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90
r2271 r2282 176 176 T2sumSTD, nlevSTD, du_gwd_rando, dv_gwd_rando, & 177 177 ulevSTD, vlevSTD, wlevSTD, philevSTD, qlevSTD, tlevSTD, & 178 missing_val_nf90 178 rhlevSTD, O3STD, O3daySTD, uvSTD, vqSTD, vTSTD, wqSTD, & 179 vphiSTD, wTSTD, u2STD, v2STD, T2STD, missing_val_nf90 179 180 180 181 USE phys_local_var_mod, only: zxfluxlat, slp, zxtsol, zt2m, & … … 663 664 ! ENDIF 664 665 666 #ifdef CPP_IOIPSL 667 #ifndef CPP_XIOS 665 668 IF (.NOT.ok_all_xml) THEN 666 669 ! ATTENTION, LES ANCIENS HISTWRITE ONT ETES CONSERVES EN ATTENDANT MIEUX: 667 670 ! Champs interpolles sur des niveaux de pression 668 #ifndef CPP_XIOS669 671 missing_val=missing_val_nf90 670 #endif671 672 DO iff=1, nfiles 672 673 ll=0 … … 690 691 ENDDO 691 692 ENDDO 692 ELSE 693 ENDIF 694 #endif 695 #endif 696 #ifdef CPP_XIOS 697 IF(ok_all_xml) THEN 693 698 !XIOS CALL xios_get_field_attr("u850",default_value=missing_val) 694 699 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 695 DO iff=1, nfiles696 700 ll=0 697 701 DO k=1, nlevSTD … … 702 706 bb2.EQ."50".OR.bb2.EQ."10") THEN 703 707 ll=ll+1 704 CALL histwrite_phy(o_uSTDlevs(ll),ulevSTD(:,k) , iff)705 CALL histwrite_phy(o_vSTDlevs(ll),vlevSTD(:,k) , iff)706 CALL histwrite_phy(o_wSTDlevs(ll),wlevSTD(:,k) , iff)707 CALL histwrite_phy(o_zSTDlevs(ll),philevSTD(:,k) , iff)708 CALL histwrite_phy(o_qSTDlevs(ll),qlevSTD(:,k) , iff)709 CALL histwrite_phy(o_tSTDlevs(ll),tlevSTD(:,k) , iff)708 CALL histwrite_phy(o_uSTDlevs(ll),ulevSTD(:,k)) 709 CALL histwrite_phy(o_vSTDlevs(ll),vlevSTD(:,k)) 710 CALL histwrite_phy(o_wSTDlevs(ll),wlevSTD(:,k)) 711 CALL histwrite_phy(o_zSTDlevs(ll),philevSTD(:,k)) 712 CALL histwrite_phy(o_qSTDlevs(ll),qlevSTD(:,k)) 713 CALL histwrite_phy(o_tSTDlevs(ll),tlevSTD(:,k)) 710 714 ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR. 711 715 ENDDO 712 ENDDO713 716 ENDIF 717 #endif 714 718 IF (vars_defined) THEN 715 719 DO i=1, klon … … 1304 1308 ENDIF 1305 1309 !!!!!!!!!!!! Sorties niveaux de pression NMC !!!!!!!!!!!!!!!!!!!! 1310 #ifdef CPP_IOIPSL 1311 #ifndef CPP_XIOS 1312 IF (.NOT.ok_all_xml) THEN 1313 ! ATTENTION, LES ANCIENS HISTWRITE ONT ETES CONSERVES EN ATTENDANT MIEUX: 1314 ! Champs interpolles sur des niveaux de pression 1315 missing_val=missing_val_nf90 1306 1316 DO iff=7, nfiles 1307 1317 … … 1367 1377 CALL histwrite_phy(o_TxT,T2sumSTD(:,:,iff-6),iff) 1368 1378 ENDDO !nfiles 1379 ENDIF 1380 #endif 1381 #endif 1382 #ifdef CPP_XIOS 1383 IF(ok_all_xml) THEN 1384 ! DO iff=7, nfiles 1385 1386 ! CALL histwrite_phy(o_tnondef,tnondef(:,:,3)) 1387 CALL histwrite_phy(o_ta,tlevSTD(:,:)) 1388 CALL histwrite_phy(o_zg,philevSTD(:,:)) 1389 CALL histwrite_phy(o_hus,qlevSTD(:,:)) 1390 CALL histwrite_phy(o_hur,rhlevSTD(:,:)) 1391 CALL histwrite_phy(o_ua,ulevSTD(:,:)) 1392 CALL histwrite_phy(o_va,vlevSTD(:,:)) 1393 CALL histwrite_phy(o_wap,wlevSTD(:,:)) 1394 ! IF(vars_defined) THEN 1395 ! DO k=1, nlevSTD 1396 ! DO i=1, klon 1397 ! IF(tnondef(i,k,3).NE.missing_val) THEN 1398 ! IF(freq_outNMC(iff-6).LT.0) THEN 1399 ! freq_moyNMC(iff-6)=(mth_len*un_jour)/freq_calNMC(iff-6) 1400 ! ELSE 1401 ! freq_moyNMC(iff-6)=freq_outNMC(iff-6)/freq_calNMC(iff-6) 1402 ! ENDIF 1403 ! zx_tmp_fi3d_STD(i,k) = (100.*tnondef(i,k,3))/freq_moyNMC(iff-6) 1404 ! ELSE 1405 ! zx_tmp_fi3d_STD(i,k) = missing_val 1406 ! ENDIF 1407 ! ENDDO 1408 ! ENDDO 1409 ! ENDIF 1410 ! CALL histwrite_phy(o_psbg,zx_tmp_fi3d_STD) 1411 IF(vars_defined) THEN 1412 DO k=1, nlevSTD 1413 DO i=1, klon 1414 IF(O3STD(i,k).NE.missing_val) THEN 1415 zx_tmp_fi3d_STD(i,k) = O3STD(i,k) * 1.e+9 1416 ELSE 1417 zx_tmp_fi3d_STD(i,k) = missing_val 1418 ENDIF 1419 ENDDO 1420 ENDDO !k=1, nlevSTD 1421 ENDIF 1422 CALL histwrite_phy(o_tro3,zx_tmp_fi3d_STD) 1423 if (read_climoz == 2) THEN 1424 IF(vars_defined) THEN 1425 DO k=1, nlevSTD 1426 DO i=1, klon 1427 IF(O3daySTD(i,k).NE.missing_val) THEN 1428 zx_tmp_fi3d_STD(i,k) = O3daySTD(i,k) * 1.e+9 1429 ELSE 1430 zx_tmp_fi3d_STD(i,k) = missing_val 1431 ENDIF 1432 ENDDO 1433 ENDDO !k=1, nlevSTD 1434 ENDIF 1435 CALL histwrite_phy(o_tro3_daylight,zx_tmp_fi3d_STD) 1436 endif 1437 CALL histwrite_phy(o_uxv,uvSTD(:,:)) 1438 CALL histwrite_phy(o_vxq,vqSTD(:,:)) 1439 CALL histwrite_phy(o_vxT,vTSTD(:,:)) 1440 CALL histwrite_phy(o_wxq,wqSTD(:,:)) 1441 CALL histwrite_phy(o_vxphi,vphiSTD(:,:)) 1442 CALL histwrite_phy(o_wxT,wTSTD(:,:)) 1443 CALL histwrite_phy(o_uxu,u2STD(:,:)) 1444 CALL histwrite_phy(o_vxv,v2STD(:,:)) 1445 CALL histwrite_phy(o_TxT,T2STD(:,:)) 1446 ! ENDDO !nfiles 1447 ENDIF 1448 #endif 1369 1449 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 1370 1450 IF (nqtot.GE.nqo+1) THEN
Note: See TracChangeset
for help on using the changeset viewer.
