Changeset 1398 for LMDZ4/trunk/libf/phylmd/phys_state_var_mod.F90

Timestamp:

Jun 4, 2010, 6:56:18 PM (14 years ago)

Author:

musat

Message:

Last corrections for CMIP5:

• Add O3 at standard level files histmthNMC.nc
• Add positive attribute "down" for vertical axes for all output files
• Replace "inst" by "ave" for hist*NMC.nc files to have time_counter and bounds for time axis (Marie-Alice's hint)
• Correct units for vertical axes : mb instead of hPa
• Add mass flux at the bottom of clouds
• Comment non initialized variables (s_capCL, s_oliqCL, s_cteiCL, s_trmb1, s_trmb2, s_trmb3) for the output files
• Geopotential field phy850, phi700, phi500, etc are modified to "geopotential height and are called z850, z700, z500, etc
• Meaning of specific humidity outputs - ovapinit and ovap - were interchanged
• Fields albs, albslw become alb1, alb2 in output files
• Correct title for rugs_* fields
• Correct units for pbase and ptop are Pa (not mb)
• Correct ndayrain field

FH/JLD/JYG/MAF/IM

File:

• LMDZ4/trunk/libf/phylmd/phys_state_var_mod.F90 (modified) (7 diffs)

LMDZ4/trunk/libf/phylmd/phys_state_var_mod.F90

@@ -117 +117 @@
       REAL,ALLOCATABLE,SAVE :: T2sumSTD(:,:,:)
 !$OMP THREADPRIVATE(T2sumSTD)
+      REAL,ALLOCATABLE,SAVE :: O3sumSTD(:,:,:), O3daysumSTD(:,:,:)
+!$OMP THREADPRIVATE(O3sumSTD,O3daysumSTD) 
 !IM begin
       REAL,ALLOCATABLE,SAVE :: wlevSTD(:,:), ulevSTD(:,:), vlevSTD(:,:)
@@ -142 +144 @@
       REAL,ALLOCATABLE,SAVE :: T2STD(:,:)
 !$OMP THREADPRIVATE(T2STD)
+      REAL,ALLOCATABLE,SAVE :: O3STD(:,:), O3daySTD(:,:)
+!$OMP THREADPRIVATE(O3STD,O3daySTD)
 !IM end
       INTEGER,ALLOCATABLE,SAVE :: seed_old(:,:)
@@ -150 +154 @@
 !$OMP THREADPRIVATE(alb_neig)
 !cloud base mass flux
-      REAL,ALLOCATABLE,SAVE :: ema_workcbmf(:), ema_cbmf(:)
-!$OMP THREADPRIVATE(ema_workcbmf,ema_cbmf)
+      REAL,ALLOCATABLE,SAVE :: ema_cbmf(:)
+!$OMP THREADPRIVATE(ema_cbmf)
 !cloud base pressure & cloud top pressure
       REAL,ALLOCATABLE,SAVE :: ema_pcb(:), ema_pct(:)
@@ -379 +383 @@
       ALLOCATE(v2sumSTD(klon,nlevSTD,nout))
       ALLOCATE(T2sumSTD(klon,nlevSTD,nout))
+      ALLOCATE(O3sumSTD(klon,nlevSTD,nout))
+      ALLOCATE(O3daysumSTD(klon,nlevSTD,nout))
 !IM beg
       ALLOCATE(wlevSTD(klon,nlevSTD), ulevSTD(klon,nlevSTD), vlevSTD(klon,nlevSTD))
@@ -388 +394 @@
       ALLOCATE(u2STD(klon,nlevSTD),v2STD(klon,nlevSTD))
       ALLOCATE(T2STD(klon,nlevSTD))
+      ALLOCATE(O3STD(klon,nlevSTD))
+      ALLOCATE(O3daySTD(klon,nlevSTD))
 !IM end
       ALLOCATE(seed_old(klon,napisccp))
@@ -393 +401 @@
       ALLOCATE(alb_neig(klon))
 !cloud base mass flux
-      ALLOCATE(ema_workcbmf(klon), ema_cbmf(klon))
+      ALLOCATE(ema_cbmf(klon))
 !cloud base pressure & cloud top pressure
       ALLOCATE(ema_pcb(klon), ema_pct(klon))
@@ -490 +498 @@
       deallocate(v2sumSTD)
       deallocate(T2sumSTD)
+      deallocate(O3sumSTD)
+      deallocate(O3daysumSTD)
 !IM beg
       deallocate(wlevSTD,ulevSTD,vlevSTD,tlevSTD,qlevSTD,rhlevSTD,philevSTD)
-      deallocate(uvSTD,vqSTD,vTSTD,wqSTD,vphiSTD,wTSTD,u2STD,v2STD,T2STD)
+      deallocate(uvSTD,vqSTD,vTSTD,wqSTD,vphiSTD,wTSTD,u2STD,v2STD,T2STD,O3STD,O3daySTD)
 !IM end
       deallocate(seed_old)
       deallocate(zuthe, zvthe)
       deallocate(alb_neig)
-      deallocate(ema_workcbmf, ema_cbmf)
+      deallocate(ema_cbmf)
       deallocate(ema_pcb, ema_pct)
       deallocate(Ma, qcondc)