Changeset 1352 for LMDZ4/trunk/libf/phylmd/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/undefSTD.F (modified) (5 diffs)

LMDZ4/trunk/libf/phylmd/undefSTD.F

@@ -2 +2 @@
 ! $Header$
 !
-      SUBROUTINE undefSTD(nlevSTD,itap,tlevSTD,
-     $           dtime,ecrit_hf,
-     $           oknondef,tnondef,tsumSTD)
+      SUBROUTINE undefSTD(itap,freq_calNMC)
       USE netcdf
       USE dimphy
+      USE phys_state_var_mod ! Variables sauvegardees de la physique
       IMPLICIT none
 c
@@ -14 +13 @@
 c
 c Calcul * du nombre de pas de temps (FLOAT(ecrit_XXX)-tnondef)) 
-c          ou la variable tlevSTD est bien definie (.NE.1.E+20), 
+c          ou la variable tlevSTD est bien definie (.NE.missing_val), 
 c et 
 c        * de la somme de tlevSTD => tsumSTD
@@ -22 +21 @@
 c nout=2 !var. mensuelle "mth" moyennee sur tous les pas de temps
 c        ! de la physique
-c nout=3 !var. mensuelle "NMC" moyennee toutes les 6heures
+c nout=3 !var. mensuelle "NMC" moyennee toutes les ecrit_hf
 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
@@ -34 +33 @@
 c variables Input
 c
-      INTEGER nlevSTD, klevSTD, itap
-      PARAMETER(klevSTD=17)
-      REAL dtime, ecrit_hf
+c     INTEGER nlevSTD, klevSTD, itap
+c     PARAMETER(klevSTD=17)
+      INTEGER itap
+c     REAL dtime
 c
 c variables locales
-      INTEGER i, k, nout
-      PARAMETER(nout=3) !nout=1 : day; =2 : mth; =3 : NMC
+c     INTEGER i, k, nout, n
+c     PARAMETER(nout=3) !nout=1 : day; =2 : mth; =3 : NMC
+      INTEGER i, k, n
+      REAL freq_calNMC(nout)
 c
 c variables Output
-      REAL tlevSTD(klon,klevSTD), tsumSTD(klon,klevSTD,nout)
-      LOGICAL oknondef(klon,klevSTD,nout)
-      REAL tnondef(klon,klevSTD,nout)
+c     REAL tlevSTD(klon,klevSTD), tsumSTD(klon,klevSTD,nout)
+c     LOGICAL oknondef(klon,klevSTD,nout)
+c     REAL tnondef(klon,klevSTD,nout)
 c
       REAL missing_val
@@ -51 +53 @@
       missing_val=nf90_fill_real
 c
-c calcul variables tous les pas de temps de la physique 
+      DO n=1, nout
 c
-      DO k=1, nlevSTD
-       DO i=1, klon
-        IF(tlevSTD(i,k).EQ.missing_val) THEN
-         IF(oknondef(i,k,1)) THEN          
-          tnondef(i,k,1)=tnondef(i,k,1)+1.
-         ENDIF !oknondef(i,k)
 c
-         IF(oknondef(i,k,2)) THEN          
-          tnondef(i,k,2)=tnondef(i,k,2)+1.
-         ENDIF !oknondef(i,k)
+c calcul variables tous les freq_calNMC(n)/dtime pas de temps 
+c de la physique
 c
-        ELSE IF(tlevSTD(i,k).NE.missing_val) THEN
-         tsumSTD(i,k,1)=tsumSTD(i,k,1)+tlevSTD(i,k)
-         tsumSTD(i,k,2)=tsumSTD(i,k,2)+tlevSTD(i,k)
-        ENDIF 
-       ENDDO !i
-      ENDDO !k
+       IF(MOD(itap,NINT(freq_calNMC(n)/dtime)).EQ.0) THEN
 c
-c calcul variables toutes les 6h
+c       print*,'n freq_cal=',n,itap,freq_calNMC(n)/dtime
 c
-      IF(MOD(itap,NINT(ecrit_hf/dtime)).EQ.0) THEN
+        DO k=1, nlevSTD
+         DO i=1, klon
+          IF(tlevSTD(i,k).EQ.missing_val) THEN
+c          IF(oknondef(i,k,n)) THEN          
+            tnondef(i,k,n)=tnondef(i,k,n)+1.
+c          ENDIF !oknondef(i,k)
 c
-       DO k=1, nlevSTD
-        DO i=1, klon
-         IF(tlevSTD(i,k).EQ.missing_val) THEN
-          IF(oknondef(i,k,3)) THEN          
-           tnondef(i,k,3)=tnondef(i,k,3)+1.
-          ENDIF !oknondef(i,k)
+          ELSE IF(tlevSTD(i,k).NE.missing_val) THEN
+           tsumSTD(i,k,n)=tsumSTD(i,k,n)+tlevSTD(i,k)
+           usumSTD(i,k,n)=usumSTD(i,k,n)+ulevSTD(i,k)
+           vsumSTD(i,k,n)=vsumSTD(i,k,n)+vlevSTD(i,k)
+           wsumSTD(i,k,n)=wsumSTD(i,k,n)+wlevSTD(i,k)
+           phisumSTD(i,k,n)=phisumSTD(i,k,n)+philevSTD(i,k)
+           qsumSTD(i,k,n)=qsumSTD(i,k,n)+qlevSTD(i,k)
+           rhsumSTD(i,k,n)=rhsumSTD(i,k,n)+rhlevSTD(i,k)
+           uvsumSTD(i,k,n)=uvsumSTD(i,k,n)+uvSTD(i,k)
+           vqsumSTD(i,k,n)=vqsumSTD(i,k,n)+vqSTD(i,k)
+           vTsumSTD(i,k,n)=vTsumSTD(i,k,n)+vTSTD(i,k)
+           wqsumSTD(i,k,n)=wqsumSTD(i,k,n)+wqSTD(i,k)
+           vphisumSTD(i,k,n)=vphisumSTD(i,k,n)+vphiSTD(i,k)
+           wTsumSTD(i,k,n)=wTsumSTD(i,k,n)+wTSTD(i,k)
+           u2sumSTD(i,k,n)=u2sumSTD(i,k,n)+u2STD(i,k)
+           v2sumSTD(i,k,n)=v2sumSTD(i,k,n)+v2STD(i,k)
+           T2sumSTD(i,k,n)=T2sumSTD(i,k,n)+T2STD(i,k)
 c
-         ELSE IF(tlevSTD(i,k).NE.missing_val) THEN
-         tsumSTD(i,k,3)=tsumSTD(i,k,3)+tlevSTD(i,k)
-         ENDIF 
-        ENDDO !i
-       ENDDO !k
-
-      ENDIF !MOD(itap,NINT(ecrit_hf/dtime)).EQ.0
+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
+          ENDIF 
+         ENDDO !i
+        ENDDO !k
+c
+       ENDIF !MOD(itap,NINT(freq_calNMC(n)/dtime)).EQ.0
+c
+      ENDDO !n
 c
       RETURN