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

Timestamp:

Apr 14, 2010, 11:41:19 AM (14 years ago)

Author:

musat

Message:

Add 3 output files for standard pressure levels AR5 exercice and flags
to manage their computation and output frequencies
histhfNMC.nc with 3 standard pressure levels
histdayNMC.nc with 8 (or may have 17) standard pressure levels
histmthNMC.nc with 17 standard pressure levels
Add 3 flags in the physiq.def file: freq_calNMC(3), freq_outNMC(3) and lev_histdayNMC
freq_calNMC(3) : computation frequency of variables on standard pressure levels

and by default has fallowing values (in fact physics' time step dtime)

freq_calNMC(1)=900.
freq_calNMC(2)=900.
freq_calNMC(3)=900.
freq_outNMC(3) : output frequency of variables on standard pressure levels

with following default values

freq_out(1) = 2592000. (30 days)
freq_out(2) = 86400. (1 day)
freq_out(3) = 21600. (6 hours)
lev_histdayNMC is 8 by default but may be switched to 17 (if we need more levels for a particular run)
IM

File:

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

LMDZ4/trunk/libf/phylmd/phys_state_var_mod.F90

@@ -8 +8 @@
       USE dimphy
       INTEGER, PARAMETER :: nlevSTD=17
+      INTEGER, PARAMETER :: nlevSTD8=8
+      INTEGER, PARAMETER :: nlevSTD3=3
       INTEGER, PARAMETER :: nout=3
       INTEGER, PARAMETER :: napisccp=1
       INTEGER, SAVE :: radpas
+!$OMP THREADPRIVATE(radpas)
       REAL, SAVE :: dtime, solaire_etat0
-!$OMP THREADPRIVATE(radpas)
 !$OMP THREADPRIVATE(dtime, solaire_etat0)
 
@@ -115 +117 @@
       REAL,ALLOCATABLE,SAVE :: T2sumSTD(:,:,:)
 !$OMP THREADPRIVATE(T2sumSTD)
+!IM begin
+      REAL,ALLOCATABLE,SAVE :: wlevSTD(:,:), ulevSTD(:,:), vlevSTD(:,:)
+      REAL,ALLOCATABLE,SAVE :: tlevSTD(:,:), qlevSTD(:,:), rhlevSTD(:,:)
+      REAL,ALLOCATABLE,SAVE :: philevSTD(:,:)
+      REAL,ALLOCATABLE,SAVE :: uvSTD(:,:)
+      REAL,ALLOCATABLE,SAVE :: vqSTD(:,:)
+      REAL,ALLOCATABLE,SAVE :: vTSTD(:,:)
+      REAL,ALLOCATABLE,SAVE :: wqSTD(:,:)
+      REAL,ALLOCATABLE,SAVE :: vphiSTD(:,:)
+      REAL,ALLOCATABLE,SAVE :: wTSTD(:,:)
+      REAL,ALLOCATABLE,SAVE :: u2STD(:,:)
+      REAL,ALLOCATABLE,SAVE :: v2STD(:,:) 
+      REAL,ALLOCATABLE,SAVE :: T2STD(:,:)
+!IM end
       INTEGER,ALLOCATABLE,SAVE :: seed_old(:,:)
 !$OMP THREADPRIVATE(seed_old)
@@ -351 +367 @@
       ALLOCATE(v2sumSTD(klon,nlevSTD,nout))
       ALLOCATE(T2sumSTD(klon,nlevSTD,nout))
+!IM beg
+      ALLOCATE(wlevSTD(klon,nlevSTD), ulevSTD(klon,nlevSTD), vlevSTD(klon,nlevSTD))
+      ALLOCATE(tlevSTD(klon,nlevSTD), qlevSTD(klon,nlevSTD), rhlevSTD(klon,nlevSTD))
+      ALLOCATE(philevSTD(klon,nlevSTD))
+      ALLOCATE(uvSTD(klon,nlevSTD),vqSTD(klon,nlevSTD))
+      ALLOCATE(vTSTD(klon,nlevSTD),wqSTD(klon,nlevSTD))
+      ALLOCATE(vphiSTD(klon,nlevSTD),wTSTD(klon,nlevSTD))
+      ALLOCATE(u2STD(klon,nlevSTD),v2STD(klon,nlevSTD))
+      ALLOCATE(T2STD(klon,nlevSTD))
+!IM end
       ALLOCATE(seed_old(klon,napisccp))
       ALLOCATE(zuthe(klon),zvthe(klon))
@@ -452 +478 @@
       deallocate(v2sumSTD)
       deallocate(T2sumSTD)
+!IM beg
+      deallocate(wlevSTD,ulevSTD,vlevSTD,tlevSTD,qlevSTD,rhlevSTD,philevSTD)
+      deallocate(uvSTD,vqSTD,vTSTD,wqSTD,vphiSTD,wTSTD,u2STD,v2STD,T2STD)
+!IM end
       deallocate(seed_old)
       deallocate(zuthe, zvthe)