Changeset 1352 for LMDZ4/trunk/libf/phylmd/ini_undefSTD.F

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/ini_undefSTD.F (modified) (3 diffs)

LMDZ4/trunk/libf/phylmd/ini_undefSTD.F

@@ -3 +3 @@
 !
 
-      SUBROUTINE ini_undefSTD(nlevSTD,itap,
-     $           dtime,ecrit_day,ecrit_mth,
-     $           tnondef,tsumSTD)
+      SUBROUTINE ini_undefSTD(itap,
+     $           freq_outNMC)
       USE dimphy
+      USE phys_state_var_mod ! Variables sauvegardees de la physique
       IMPLICIT none
 c
@@ -27 +27 @@
 c nout=3 !var. mensuelle "NMC" moyennee toutes les 6heures
 c
-c
-c NB: mettre "inst(X)" dans le write_histXXX.h !
+c NB: mettre "inst(X)" dans le write_hist*NMC.h !
 c====================================================================
 c
@@ -36 +35 @@
 cym #include "dimphy.h"
 c variables Input/Output
-      INTEGER nlevSTD, klevSTD, itap
-      PARAMETER(klevSTD=17)
-      REAL dtime
-      REAL ecrit_day,ecrit_mth
+c     INTEGER nlevSTD, klevSTD, itap
+      INTEGER itap
+c     PARAMETER(klevSTD=17)
+c     REAL dtime
 c
 c variables locales
-      INTEGER i, k, nout
-      PARAMETER(nout=3) !nout=1 day/nout=2 mth/nout=3 NMC
+c     INTEGER i, k, nout, n
+      INTEGER i, k, n
+c     PARAMETER(nout=3) !nout=1 day/nout=2 mth/nout=3 NMC
+      REAL freq_outNMC(nout)
 c
 c variables Output
-      REAL tnondef(klon,klevSTD,nout)
-      REAL tsumSTD(klon,klevSTD,nout)
+c     REAL tnondef(klon,klevSTD,nout)
+c     REAL tsumSTD(klon,klevSTD,nout)
 c
-c initialisation variables journalieres en debut de journee
+      DO n=1, nout
 c
-      IF(MOD(itap,NINT(ecrit_day/dtime)).EQ.1.) THEN
-       DO k=1, nlevSTD
-        DO i=1, klon
-         tnondef(i,k,1)=0.
-         tsumSTD(i,k,1)=0.
-        ENDDO !i
-       ENDDO !k
-      ENDIF
+c initialisation variables en debut de la journee ou du mois
 c
-c initialisation variables mensuelles (calculees a chaque pas de temps) 
-c en debut de mois : nout=2
+       IF(MOD(itap,NINT(freq_outNMC(n)/dtime)).EQ.1.) THEN
 c
-      IF(MOD(itap,NINT(ecrit_mth/dtime)).EQ.1.) THEN
+c       print*,'n freq_ini=',n,itap,freq_outNMC(n)/dtime
 c
-       DO k=1, nlevSTD
-        DO i=1, klon
-         tnondef(i,k,2)=0.
-         tsumSTD(i,k,2)=0.
-        ENDDO !i
-       ENDDO !k
+        DO k=1, nlevSTD
+         DO i=1, klon
+          tnondef(i,k,n)=0.
+          tsumSTD(i,k,n)=0.
+          usumSTD(i,k,n)=0.
+          vsumSTD(i,k,n)=0.
+          wsumSTD(i,k,n)=0.
+          phisumSTD(i,k,n)=0.
+          qsumSTD(i,k,n)=0.
+          rhsumSTD(i,k,n)=0.
+          uvsumSTD(i,k,n)=0.
+          vqsumSTD(i,k,n)=0.
+          vTsumSTD(i,k,n)=0.
+          wqsumSTD(i,k,n)=0.
+          vphisumSTD(i,k,n)=0.
+          wTsumSTD(i,k,n)=0.
+          u2sumSTD(i,k,n)=0.
+          v2sumSTD(i,k,n)=0.
+          T2sumSTD(i,k,n)=0.
+         ENDDO !i
+        ENDDO !k
 c
-c initialisation variables mensuelles - runs type Amip - (calculees toutes les 6h)
-c en debut de mois : nout = 3
+c          if(n.EQ.1.AND.i.EQ.4513.AND.k.EQ.17) THEN
+           if(n.EQ.1.AND.i.EQ.1128.AND.k.EQ.17) THEN
+            print*,'itap rlon rlat tlevSTD=',itap,rlon(i),rlat(i),
+     $tlevSTD(i,k)
+           endif
 c
-       DO k=1, nlevSTD
-        DO i=1, klon
-         tnondef(i,k,3)=0.
-         tsumSTD(i,k,3)=0.
-        ENDDO !i
-       ENDDO !k
+       ENDIF !MOD(itap,NINT(freq_outNMC(n)/dtime)).EQ.1.
 c
-      ENDIF
+      ENDDO !n
 c
       RETURN