Changeset 1400

Timestamp:

Jun 14, 2010, 5:03:34 PM (15 years ago)

Author:

musat

Message:

Change units for vertical axes from "mb" to (SI) "Pa"
Change titles for variables at standard pressure levels from "mb" to "hPa"
IM

Location:

LMDZ4/trunk/libf/phylmd

Files:

• ini_histdayNMC.h (modified) (1 diff)
• ini_histhfNMC.h (modified) (1 diff)
• ini_histmthNMC.h (modified) (1 diff)
• phys_output_mod.F90 (modified) (2 diffs)

LMDZ4/trunk/libf/phylmd/ini_histdayNMC.h

@@ -30 +30 @@
 c
         IF(lev_histdayNMC.EQ.nlevSTD) THEN 
-         CALL histvert(nid_daynmc, "plev", "pressure", "mb",
-     .                 nlevSTD, rlevSTD/100., nvert,"down")
+         CALL histvert(nid_daynmc, "plev", "pressure", "Pa",
+     .                 nlevSTD, rlevSTD, nvert,"down")
         ELSE IF(lev_histdayNMC.EQ.nlevSTD8) THEN 
-         CALL histvert(nid_daynmc, "plev", "pressure", "mb",
-     .                 nlevSTD8, rlevSTD8/100., nvert,"down")
+         CALL histvert(nid_daynmc, "plev", "pressure", "Pa",
+     .                 nlevSTD8, rlevSTD8, nvert,"down")
         ENDIF
 c

LMDZ4/trunk/libf/phylmd/ini_histhfNMC.h

@@ -31 +31 @@
      .                 nhori, nid_hfnmc)
 c
-         CALL histvert(nid_hfnmc, "plev", "pressure", "mb",
-     .                 nlevSTD3, rlevSTD3/100., nvert,"down")
+         CALL histvert(nid_hfnmc, "plev", "pressure", "Pa",
+     .                 nlevSTD3, rlevSTD3, nvert,"down")
 c
 cIM Astuce MAF: remplacer inst par ave pour les variables NMC pour avoir

LMDZ4/trunk/libf/phylmd/ini_histmthNMC.h

@@ -27 +27 @@
      .                 nhori, nid_mthnmc)
 c
-         CALL histvert(nid_mthnmc, "plev", "pressure", "mb",
-     .                 nlevSTD, rlevSTD/100., nvert,"down")
+         CALL histvert(nid_mthnmc, "plev", "pressure", "Pa",
+     .                 nlevSTD, rlevSTD, nvert,"down")
 c
 cIM Astuce MAF: remplacer inst par ave pour les variables NMC pour avoir

LMDZ4/trunk/libf/phylmd/phys_output_mod.F90

@@ -745 +745 @@
        endif
  
-      CALL histvert(nid_files(iff), "presnivs", "Vertical levels", "mb", &
+      CALL histvert(nid_files(iff), "presnivs", "Vertical levels", "Pa", &
            levmax(iff) - levmin(iff) + 1, &
-           presnivs(levmin(iff):levmax(iff))/100., nvertm(iff),"down")
+           presnivs(levmin(iff):levmax(iff)), nvertm(iff),"down")
 
 !!!!!!!!!!!!! Traitement des champs 3D pour histhf !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -1041 +1041 @@
       ilev=ilev+1
       print*,'ilev k bb2 flag name ',ilev,k, bb2,o_uSTDlevs(ilev)%flag,o_uSTDlevs(ilev)%name
- CALL histdef2d(iff,o_uSTDlevs(ilev)%flag,o_uSTDlevs(ilev)%name,"Zonal wind "//bb2//"mb", "m/s")
- CALL histdef2d(iff,o_vSTDlevs(ilev)%flag,o_vSTDlevs(ilev)%name,"Meridional wind "//bb2//"mb", "m/s")
- CALL histdef2d(iff,o_wSTDlevs(ilev)%flag,o_wSTDlevs(ilev)%name,"Vertical wind "//bb2//"mb", "Pa/s")
- CALL histdef2d(iff,o_zSTDlevs(ilev)%flag,o_zSTDlevs(ilev)%name,"Geopotential height "//bb2//"mb", "m")
- CALL histdef2d(iff,o_qSTDlevs(ilev)%flag,o_qSTDlevs(ilev)%name,"Specific humidity "//bb2//"mb", "kg/kg" )
- CALL histdef2d(iff,o_tSTDlevs(ilev)%flag,o_tSTDlevs(ilev)%name,"Temperature "//bb2//"mb", "K")
+ CALL histdef2d(iff,o_uSTDlevs(ilev)%flag,o_uSTDlevs(ilev)%name,"Zonal wind "//bb2//"hPa", "m/s")
+ CALL histdef2d(iff,o_vSTDlevs(ilev)%flag,o_vSTDlevs(ilev)%name,"Meridional wind "//bb2//"hPa", "m/s")
+ CALL histdef2d(iff,o_wSTDlevs(ilev)%flag,o_wSTDlevs(ilev)%name,"Vertical wind "//bb2//"hPa", "Pa/s")
+ CALL histdef2d(iff,o_zSTDlevs(ilev)%flag,o_zSTDlevs(ilev)%name,"Geopotential height "//bb2//"hPa", "m")
+ CALL histdef2d(iff,o_qSTDlevs(ilev)%flag,o_qSTDlevs(ilev)%name,"Specific humidity "//bb2//"hPa", "kg/kg" )
+ CALL histdef2d(iff,o_tSTDlevs(ilev)%flag,o_tSTDlevs(ilev)%name,"Temperature "//bb2//"hPa", "K")
      ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR."500".OR.bb2.EQ."200".OR.bb2.EQ."50".OR.bb2.EQ."10")
        ENDDO