Changeset 644 for LMDZ4/trunk/libf/phylmd/ini_histmthNMC.h
- Timestamp:
- May 25, 2005, 3:10:09 PM (19 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/trunk/libf/phylmd/ini_histmthNMC.h
r524 r644 4 4 IF (ok_mensuel) THEN 5 5 c 6 zsto = dtime 7 zout = dtime * ecrit_mth 6 zstophy = dtime 7 zstohf = dtime * FLOAT(ecrit_hf) 8 zstomth = dtime * FLOAT(ecrit_mth) 9 zout = dtime * FLOAT(ecrit_mth) 8 10 c 9 11 idayref = day_ref … … 22 24 . 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 23 25 . nhori, nid_nmc) 24 write(*,*)'Mensuel NMC ', itau_phy, zjulian 26 c 25 27 CALL histvert(nid_nmc, "presnivs", "Vertical levels", "mb", 26 . klev, presnivs, nvert) 27 c call histvert(nid_nmc, 'sig_s', 'Niveaux sigma','-', 28 c . klev, znivsig, nvert) 28 . nlevSTD, rlevSTD/100., nvert) 29 ccc 30 ccc Champs 3D interpolles sur des niveaux de pression du NMC 31 ccc 32 IF(type_run.EQ."ENSP".OR.type_run.EQ."CLIM") THEN 29 33 c 30 c Champs 2D: 31 32 c Champs interpolles sur des niveaux de pression du NMC 33 c IMIMIM 110304 BEG 34 CALL histdef(nid_nmc, "temp", 35 . "Temperature","K", 36 . iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32, 37 . "inst(X)", zout,zout) 34 38 c 35 DO k=1, nlevSTD 39 CALL histdef(nid_nmc, "phi", 40 . "Geopotential", "m", 41 . iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32, 42 . "inst(X)", zout,zout) 36 43 c 37 bb=clevSTD(k) 44 CALL histdef(nid_nmc, "q", 45 . "Specific humidity","kg/kg", 46 . iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32, 47 . "inst(X)", zout,zout) 38 48 c 39 IF(k.GE.2) THEN40 aa=clevSTD(k)41 bb=aa(1:lnblnk1(aa))42 ENDIF49 CALL histdef(nid_nmc, "rh", 50 . "Relative humidity", "%", 51 . iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32, 52 . "inst(X)", zout,zout) 43 53 c 44 CALL histdef(nid_nmc, " t"//bb,45 . " Temperature"//bb//"mb","K",46 . iim,jjmp1,nhori, 1,1,1, -99, 32,47 . " ave(X)", zsto,zout)54 CALL histdef(nid_nmc, "u", 55 . "Zonal wind","m/s", 56 . iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32, 57 . "inst(X)", zout,zout) 48 58 c 49 CALL histdef(nid_nmc, "phi"//bb,50 . "Geopotential"//bb//"mb", "m2/s2",51 . iim,jjmp1,nhori, 1,1,1, -99, 32,52 . " ave(X)", zsto,zout)59 CALL histdef(nid_nmc, "v", 60 . "Meridional wind","m/s", 61 . iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32, 62 . "inst(X)", zout,zout) 53 63 c 54 CALL histdef(nid_nmc, "q"//bb, 55 . "Specific humidity"//bb//"mb","kg/kg", 56 . iim,jjmp1,nhori, 1,1,1, -99, 32, 57 . "ave(X)", zsto,zout) 64 ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN 58 65 c 59 CALL histdef(nid_nmc, "rh"//bb, 60 . "Relative humidity"//bb//"mb", "%", 61 . iim,jjmp1,nhori, 1,1,1, -99, 32, 62 . "ave(X)", zsto,zout) 66 c ATTENTION : pour AMIP2 on interpole t,u,v,wphi,q,rh 67 c sur les niveaux du NMC et on somme & moyenne 68 c toutes les 6 heures par des routines undefSTD et 69 c moy_undefSTD pour eliminer les valeurs "undef" 70 c de la moyenne mensuelle 71 c ======> le "inst(X)" ci-dessous est par consequence factice ! 63 72 c 64 CALL histdef(nid_nmc, " u"//bb,65 . " Zonal wind"//bb//"mb","K",66 . iim,jjmp1,nhori, 1,1,1, -99, 32,67 . " ave(X)", zsto,zout)73 CALL histdef(nid_nmc, "temp", 74 . "Temperature","K", 75 . iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32, 76 . "inst(X)", zout,zout) 68 77 c 69 CALL histdef(nid_nmc, "v"//bb,70 . "Meridional wind"//bb//"mb","K",71 . iim,jjmp1,nhori, 1,1,1, -99, 32,72 . " ave(X)", zsto,zout)78 CALL histdef(nid_nmc, "phi", 79 . "Geopotential ", "m", 80 . iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32, 81 . "inst(X)", zout,zout) 73 82 c 74 ENDDO !nlevSTD 75 c IMIMIM 110304 END 83 CALL histdef(nid_nmc, "q", 84 . "Specific humidity","kg/kg", 85 . iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32, 86 . "inst(X)", zout,zout) 87 c 88 CALL histdef(nid_nmc, "rh", 89 . "Relative humidity", "%", 90 . iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32, 91 . "inst(X)", zout,zout) 92 c 93 CALL histdef(nid_nmc, "u", 94 . "Zonal wind","m/s", 95 . iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32, 96 . "inst(X)", zout,zout) 97 c 98 CALL histdef(nid_nmc, "v", 99 . "Meridional wind","m/s", 100 . iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32, 101 . "inst(X)", zout,zout) 102 c 103 CALL histdef(nid_nmc, "w", 104 . "Vertical motion","Pa/s", 105 . iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32, 106 . "inst(X)", zout,zout) 107 c 108 c ATTENTION : pour AMIP2 on interpole t,u,v,wphi,q,rh 109 c sur les niveaux du NMC et on somme & moyenne 110 c toutes les 6 heures par des routines undefSTD et 111 c moy_undefSTD pour eliminer les valeurs "undef" 112 c de la moyenne mensuelle 113 c ======> le "inst(X)" ci-dessus est par consequence factice ! 114 c 115 c 116 CALL histdef(nid_nmc, "psbg", 117 . "Pressure sfce below ground","%", 118 . iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32, 119 . "inst(X)", zout,zout) 120 c 121 CALL histdef(nid_nmc, "uv", 122 . "uv ", 123 . "m2/s2",iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32, 124 . "inst(X)", zout,zout) 125 c 126 CALL histdef(nid_nmc, "vq", 127 . "vq ", 128 . "m/s * (kg/kg)",iim,jjmp1,nhori, 129 . nlevSTD,1,nlevSTD, nvert, 32, 130 . "inst(X)", zout,zout) 131 c 132 CALL histdef(nid_nmc, "vT", 133 . "vT ", 134 . "mK/s",iim,jjmp1,nhori, 135 . nlevSTD,1,nlevSTD, nvert, 32, 136 . "inst(X)", zout,zout) 137 c 138 CALL histdef(nid_nmc, "wq", 139 . "wq ", 140 . "(Pa/s)*(kg/kg)",iim,jjmp1,nhori, 141 . nlevSTD,1,nlevSTD, nvert, 32, 142 . "inst(X)", zout,zout) 143 c 144 CALL histdef(nid_nmc, "vphi", 145 . "vphi ", 146 . "m2/s",iim,jjmp1,nhori, 147 . nlevSTD,1,nlevSTD, nvert, 32, 148 . "inst(X)", zout,zout) 149 c 150 CALL histdef(nid_nmc, "wT", 151 . "wT ", 152 . "K*Pa/s",iim,jjmp1,nhori, 153 . nlevSTD,1,nlevSTD, nvert, 32, 154 . "inst(X)", zout,zout) 155 c 156 CALL histdef(nid_nmc, "uxu", 157 . "u2 ", 158 . "m2/s2",iim,jjmp1,nhori, 159 . nlevSTD,1,nlevSTD, nvert, 32, 160 . "inst(X)", zout,zout) 161 c 162 CALL histdef(nid_nmc, "vxv", 163 . "v2 ", 164 . "m2/s2",iim,jjmp1,nhori, 165 . nlevSTD,1,nlevSTD, nvert, 32, 166 . "inst(X)", zout,zout) 167 c 168 CALL histdef(nid_nmc, "TxT", 169 . "T2 ", 170 . "K2",iim,jjmp1,nhori, 171 . nlevSTD,1,nlevSTD, nvert, 32, 172 . "inst(X)", zout,zout) 173 c 174 ENDIF !(type_run.EQ."AMIP") 76 175 77 176 CALL histend(nid_nmc) 78 177 c 79 c ndex2d = 080 c81 178 ENDIF ! fin de test sur ok_mensuel 82
Note: See TracChangeset
for help on using the changeset viewer.