Changeset 1350 for LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/ini_undefSTD.F
- Timestamp:
- Apr 14, 2010, 10:55:22 AM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/ini_undefSTD.F
r776 r1350 3 3 ! 4 4 5 SUBROUTINE ini_undefSTD(nlevSTD,itap, 6 $ dtime,ecrit_day,ecrit_mth, 7 $ tnondef,tsumSTD) 5 SUBROUTINE ini_undefSTD(itap, 6 $ freq_outNMC) 8 7 USE dimphy 8 USE phys_state_var_mod ! Variables sauvegardees de la physique 9 9 IMPLICIT none 10 10 c … … 27 27 c nout=3 !var. mensuelle "NMC" moyennee toutes les 6heures 28 28 c 29 c 30 c NB: mettre "inst(X)" dans le write_histXXX.h ! 29 c NB: mettre "inst(X)" dans le write_hist*NMC.h ! 31 30 c==================================================================== 32 31 c … … 36 35 cym #include "dimphy.h" 37 36 c variables Input/Output 38 39 PARAMETER(klevSTD=17)40 REAL dtime 41 REAL ecrit_day,ecrit_mth 37 c INTEGER nlevSTD, klevSTD, itap 38 INTEGER itap 39 c PARAMETER(klevSTD=17) 40 c REAL dtime 42 41 c 43 42 c variables locales 44 INTEGER i, k, nout 45 PARAMETER(nout=3) !nout=1 day/nout=2 mth/nout=3 NMC 43 c INTEGER i, k, nout, n 44 INTEGER i, k, n 45 c PARAMETER(nout=3) !nout=1 day/nout=2 mth/nout=3 NMC 46 REAL freq_outNMC(nout) 46 47 c 47 48 c variables Output 48 49 49 c REAL tnondef(klon,klevSTD,nout) 50 c REAL tsumSTD(klon,klevSTD,nout) 50 51 c 51 c initialisation variables journalieres en debut de journee 52 DO n=1, nout 52 53 c 53 IF(MOD(itap,NINT(ecrit_day/dtime)).EQ.1.) THEN 54 DO k=1, nlevSTD 55 DO i=1, klon 56 tnondef(i,k,1)=0. 57 tsumSTD(i,k,1)=0. 58 ENDDO !i 59 ENDDO !k 60 ENDIF 54 c initialisation variables en debut de la journee ou du mois 61 55 c 62 c initialisation variables mensuelles (calculees a chaque pas de temps) 63 c en debut de mois : nout=2 56 IF(MOD(itap,NINT(freq_outNMC(n)/dtime)).EQ.1.) THEN 64 57 c 65 IF(MOD(itap,NINT(ecrit_mth/dtime)).EQ.1.) THEN 58 c print*,'n freq_ini=',n,itap,freq_outNMC(n)/dtime 66 59 c 67 DO k=1, nlevSTD 68 DO i=1, klon 69 tnondef(i,k,2)=0. 70 tsumSTD(i,k,2)=0. 71 ENDDO !i 72 ENDDO !k 60 DO k=1, nlevSTD 61 DO i=1, klon 62 tnondef(i,k,n)=0. 63 tsumSTD(i,k,n)=0. 64 usumSTD(i,k,n)=0. 65 vsumSTD(i,k,n)=0. 66 wsumSTD(i,k,n)=0. 67 phisumSTD(i,k,n)=0. 68 qsumSTD(i,k,n)=0. 69 rhsumSTD(i,k,n)=0. 70 uvsumSTD(i,k,n)=0. 71 vqsumSTD(i,k,n)=0. 72 vTsumSTD(i,k,n)=0. 73 wqsumSTD(i,k,n)=0. 74 vphisumSTD(i,k,n)=0. 75 wTsumSTD(i,k,n)=0. 76 u2sumSTD(i,k,n)=0. 77 v2sumSTD(i,k,n)=0. 78 T2sumSTD(i,k,n)=0. 79 ENDDO !i 80 ENDDO !k 73 81 c 74 c initialisation variables mensuelles - runs type Amip - (calculees toutes les 6h) 75 c en debut de mois : nout = 3 82 c if(n.EQ.1.AND.i.EQ.4513.AND.k.EQ.17) THEN 83 if(n.EQ.1.AND.i.EQ.1128.AND.k.EQ.17) THEN 84 print*,'itap rlon rlat tlevSTD=',itap,rlon(i),rlat(i), 85 $tlevSTD(i,k) 86 endif 76 87 c 77 DO k=1, nlevSTD 78 DO i=1, klon 79 tnondef(i,k,3)=0. 80 tsumSTD(i,k,3)=0. 81 ENDDO !i 82 ENDDO !k 88 ENDIF !MOD(itap,NINT(freq_outNMC(n)/dtime)).EQ.1. 83 89 c 84 END IF90 ENDDO !n 85 91 c 86 92 RETURN
Note: See TracChangeset
for help on using the changeset viewer.