- Timestamp:
- Mar 24, 2004, 9:59:48 AM (21 years ago)
- Location:
- LMDZ.3.3/branches/rel-LF/libf/phylmd
- Files:
-
- 2 edited
-
ini_histmthNMC.h (modified) (2 diffs)
-
write_histmthNMC.h (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histmthNMC.h
r500 r508 4 4 IF(region.EQ.'3d') THEN 5 5 zsto = dtime 6 cIM 150304 BEG 6 7 zout = dtime * ecrit_mth 8 c zout = dtime * REAL(NINT(86400./dtime*ecrit_hf)) 9 cIM 150304 END 10 cIM 120304 BEG 11 c zout = dtime * REAL(NINT(86400./dtime*ecrit_hf)) 12 c zsto = zout 13 c zsto = dtime 14 cIM 120304 END 7 15 ELSEIF(region.EQ.'1d') THEN 8 16 zsto = dtime … … 40 48 41 49 c Champs interpolles sur des niveaux de pression du NMC 42 CALL histdef(nid_nmc, "u1000", 43 . "Zonal wind 1000mb", 44 . "m/s", 45 . iim,jjmp1,nhori, 1,1,1, -99, 32, 46 . "ave(X)", zsto,zout) 47 48 CALL histdef(nid_nmc, "v1000", 49 . "Meridional wind 1000mb", "m/s", 50 . iim,jjmp1,nhori, 1,1,1, -99, 32, 51 . "ave(X)", zsto,zout) 52 53 CALL histdef(nid_nmc, "u925", "Zonal wind 925mb", 54 . "m/s", 55 . iim,jjmp1,nhori, 1,1,1, -99, 32, 56 . "ave(X)", zsto,zout) 57 58 CALL histdef(nid_nmc, "v925", "Meridional wind 925mb", 59 . "m/s", 60 . iim,jjmp1,nhori, 1,1,1, -99, 32, 61 . "ave(X)", zsto,zout) 62 63 CALL histdef(nid_nmc, "u850", "Zonal wind 850mb", 64 . "m/s", 65 . iim,jjmp1,nhori, 1,1,1, -99, 32, 66 . "ave(X)", zsto,zout) 67 68 CALL histdef(nid_nmc, "v850", "Meridional wind 850mb", 69 . "m/s", 50 c IMIMIM 110304 BEG 51 c 52 DO k=1, nlevSTD 53 c 54 bb=clevSTD(k) 55 c 56 IF(k.GE.2) THEN 57 aa=clevSTD(k) 58 bb=aa(1:lnblnk1(aa)) 59 ENDIF 60 c 61 CALL histdef(nid_nmc, "t"//bb, 62 . "Temperature"//bb//"mb","K", 70 63 . iim,jjmp1,nhori, 1,1,1, -99, 32, 71 64 . "ave(X)", zsto,zout) 72 65 c 73 CALL histdef(nid_nmc, "u700", "Zonal wind 700mb", 74 . "m/s", 75 . iim,jjmp1,nhori, 1,1,1, -99, 32, 76 . "ave(X)", zsto,zout) 77 78 CALL histdef(nid_nmc, "v700", "Meridional wind 700mb", 79 . "m/s", 66 CALL histdef(nid_nmc, "phi"//bb, 67 . "Geopotential"//bb//"mb", "m2/s2", 80 68 . iim,jjmp1,nhori, 1,1,1, -99, 32, 81 69 . "ave(X)", zsto,zout) 82 70 c 83 CALL histdef(nid_nmc, "u600", "Zonal wind 600mb", 84 . "m/s", 85 . iim,jjmp1,nhori, 1,1,1, -99, 32, 86 . "ave(X)", zsto,zout) 87 88 CALL histdef(nid_nmc, "v600", "Meridional wind 600mb", 89 . "m/s", 71 CALL histdef(nid_nmc, "q"//bb, 72 . "Specific humidity"//bb//"mb","kg/kg", 90 73 . iim,jjmp1,nhori, 1,1,1, -99, 32, 91 74 . "ave(X)", zsto,zout) 92 75 c 93 CALL histdef(nid_nmc, "u500", "Zonal wind 500mb", 94 . "m/s", 95 . iim,jjmp1,nhori, 1,1,1, -99, 32, 96 . "ave(X)", zsto,zout) 97 98 CALL histdef(nid_nmc, "v500", "Meridional wind 500mb", 99 . "m/s", 76 CALL histdef(nid_nmc, "rh"//bb, 77 . "Relative humidity"//bb//"mb", "%", 100 78 . iim,jjmp1,nhori, 1,1,1, -99, 32, 101 79 . "ave(X)", zsto,zout) 102 80 c 103 CALL histdef(nid_nmc, "u400", "Zonal wind 400mb", 104 . "m/s", 105 . iim,jjmp1,nhori, 1,1,1, -99, 32, 106 . "ave(X)", zsto,zout) 107 108 CALL histdef(nid_nmc, "v400", "Meridional wind 400mb", 109 . "m/s", 81 CALL histdef(nid_nmc, "u"//bb, 82 . "Zonal wind"//bb//"mb","K", 110 83 . iim,jjmp1,nhori, 1,1,1, -99, 32, 111 84 . "ave(X)", zsto,zout) 112 85 c 113 CALL histdef(nid_nmc, "u300", "Zonal wind 300mb", 114 . "m/s", 115 . iim,jjmp1,nhori, 1,1,1, -99, 32, 116 . "ave(X)", zsto,zout) 117 118 CALL histdef(nid_nmc, "v300", "Meridional wind 300mb", 119 . "m/s", 86 CALL histdef(nid_nmc, "v"//bb, 87 . "Meridional wind"//bb//"mb","K", 120 88 . iim,jjmp1,nhori, 1,1,1, -99, 32, 121 89 . "ave(X)", zsto,zout) 122 90 c 123 CALL histdef(nid_nmc, "u250", "Zonal wind 250mb", 124 . "m/s", 125 . iim,jjmp1,nhori, 1,1,1, -99, 32, 126 . "ave(X)", zsto,zout) 91 ENDDO !nlevSTD 92 c IMIMIM 110304 END 127 93 128 CALL histdef(nid_nmc, "v250", "Meridional wind 250mb",129 . "m/s",130 . iim,jjmp1,nhori, 1,1,1, -99, 32,131 . "ave(X)", zsto,zout)132 c133 CALL histdef(nid_nmc, "u200", "Zonal wind 200mb",134 . "m/s",135 . iim,jjmp1,nhori, 1,1,1, -99, 32,136 . "ave(X)", zsto,zout)137 138 CALL histdef(nid_nmc, "v200", "Meridional wind 200mb",139 . "m/s",140 . iim,jjmp1,nhori, 1,1,1, -99, 32,141 . "ave(X)", zsto,zout)142 c143 CALL histdef(nid_nmc, "u150", "Zonal wind 150mb",144 . "m/s",145 . iim,jjmp1,nhori, 1,1,1, -99, 32,146 . "ave(X)", zsto,zout)147 148 CALL histdef(nid_nmc, "v150", "Meridional wind 150mb",149 . "m/s",150 . iim,jjmp1,nhori, 1,1,1, -99, 32,151 . "ave(X)", zsto,zout)152 c153 CALL histdef(nid_nmc, "u100", "Zonal wind 100mb",154 . "m/s",155 . iim,jjmp1,nhori, 1,1,1, -99, 32,156 . "ave(X)", zsto,zout)157 158 CALL histdef(nid_nmc, "v100", "Meridional wind 100mb",159 . "m/s",160 . iim,jjmp1,nhori, 1,1,1, -99, 32,161 . "ave(X)", zsto,zout)162 c163 CALL histdef(nid_nmc, "u70", "Zonal wind 70mb",164 . "m/s",165 . iim,jjmp1,nhori, 1,1,1, -99, 32,166 . "ave(X)", zsto,zout)167 168 CALL histdef(nid_nmc, "v70", "Meridional wind 70mb",169 . "m/s",170 . iim,jjmp1,nhori, 1,1,1, -99, 32,171 . "ave(X)", zsto,zout)172 c173 CALL histdef(nid_nmc, "u50", "Zonal wind 50mb",174 . "m/s",175 . iim,jjmp1,nhori, 1,1,1, -99, 32,176 . "ave(X)", zsto,zout)177 178 CALL histdef(nid_nmc, "v50", "Meridional wind 50mb",179 . "m/s",180 . iim,jjmp1,nhori, 1,1,1, -99, 32,181 . "ave(X)", zsto,zout)182 c183 CALL histdef(nid_nmc, "u30", "Zonal wind 30mb",184 . "m/s",185 . iim,jjmp1,nhori, 1,1,1, -99, 32,186 . "ave(X)", zsto,zout)187 188 CALL histdef(nid_nmc, "v30", "Meridional wind 30mb",189 . "m/s",190 . iim,jjmp1,nhori, 1,1,1, -99, 32,191 . "ave(X)", zsto,zout)192 c193 CALL histdef(nid_nmc, "u20", "Zonal wind 20mb",194 . "m/s",195 . iim,jjmp1,nhori, 1,1,1, -99, 32,196 . "ave(X)", zsto,zout)197 198 CALL histdef(nid_nmc, "v20", "Meridional wind 20mb",199 . "m/s",200 . iim,jjmp1,nhori, 1,1,1, -99, 32,201 . "ave(X)", zsto,zout)202 c203 CALL histdef(nid_nmc, "u10", "Zonal wind 10mb",204 . "m/s",205 . iim,jjmp1,nhori, 1,1,1, -99, 32,206 . "ave(X)", zsto,zout)207 208 CALL histdef(nid_nmc, "v10", "Meridional wind 10mb",209 . "m/s",210 . iim,jjmp1,nhori, 1,1,1, -99, 32,211 . "ave(X)", zsto,zout)212 c213 94 CALL histend(nid_nmc) 214 95 c -
LMDZ.3.3/branches/rel-LF/libf/phylmd/write_histmthNMC.h
r500 r508 15 15 c 16 16 c Champs interpolles sur des niveaux de pression du NMC 17 c110304 BEG 18 DO k=1, nlevSTD 17 19 c 18 CALL gr_fi_ecrit(1, klon,iim,jjmp1, u1000,zx_tmp_2d) 19 CALL histwrite(nid_nmc, 20 ."u1000",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 20 bb=clevSTD(k) 21 21 c 22 CALL gr_fi_ecrit(1, klon,iim,jjmp1, v1000,zx_tmp_2d) 23 CALL histwrite(nid_nmc, 24 ."v1000",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 22 IF(k.GE.2) THEN 23 aa=clevSTD(k) 24 bb=aa(1:lnblnk1(aa)) 25 ENDIF 25 26 c 26 CALL gr_fi_ecrit(1, klon,iim,jjmp1, u925,zx_tmp_2d)27 CALL histwrite(nid_nmc, 28 ."u925",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)27 CALL gr_fi_ecrit(1, klon,iim,jjmp1, tlevSTD(:,k),zx_tmp_2d) 28 CALL histwrite(nid_nmc,"t"//bb,itau_w,zx_tmp_2d, 29 $ iim*jjmp1,ndex2d) 29 30 c 30 CALL gr_fi_ecrit(1, klon,iim,jjmp1, v925,zx_tmp_2d)31 CALL histwrite(nid_nmc, 32 ."v925",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)31 CALL gr_fi_ecrit(1, klon,iim,jjmp1, philevSTD(:,k),zx_tmp_2d) 32 CALL histwrite(nid_nmc,"phi"//bb,itau_w,zx_tmp_2d, 33 $ iim*jjmp1,ndex2d) 33 34 c 34 CALL gr_fi_ecrit(1, klon,iim,jjmp1, u850,zx_tmp_2d)35 CALL histwrite(nid_nmc, 36 ."u850",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)35 CALL gr_fi_ecrit(1, klon,iim,jjmp1, qlevSTD(:,k),zx_tmp_2d) 36 CALL histwrite(nid_nmc,"q"//bb,itau_w,zx_tmp_2d, 37 $ iim*jjmp1,ndex2d) 37 38 c 38 CALL gr_fi_ecrit(1, klon,iim,jjmp1, v850,zx_tmp_2d)39 CALL histwrite(nid_nmc," v850",itau_w,40 .zx_tmp_2d,iim*jjmp1,ndex2d)39 CALL gr_fi_ecrit(1, klon,iim,jjmp1, rhlevSTD(:,k),zx_tmp_2d) 40 CALL histwrite(nid_nmc,"rh"//bb,itau_w,zx_tmp_2d, 41 $ iim*jjmp1,ndex2d) 41 42 c 42 CALL gr_fi_ecrit(1, klon,iim,jjmp1, u 700,zx_tmp_2d)43 CALL histwrite(nid_nmc,"u 700",itau_w,44 .zx_tmp_2d,iim*jjmp1,ndex2d)43 CALL gr_fi_ecrit(1, klon,iim,jjmp1, ulevSTD(:,k),zx_tmp_2d) 44 CALL histwrite(nid_nmc,"u"//bb,itau_w,zx_tmp_2d, 45 $ iim*jjmp1,ndex2d) 45 46 c 46 CALL gr_fi_ecrit(1, klon,iim,jjmp1, v 700,zx_tmp_2d)47 CALL histwrite(nid_nmc,"v 700",itau_w,48 .zx_tmp_2d,iim*jjmp1,ndex2d)47 CALL gr_fi_ecrit(1, klon,iim,jjmp1, vlevSTD(:,k),zx_tmp_2d) 48 CALL histwrite(nid_nmc,"v"//bb,itau_w,zx_tmp_2d, 49 $ iim*jjmp1,ndex2d) 49 50 c 50 CALL gr_fi_ecrit(1, klon,iim,jjmp1, u600,zx_tmp_2d)51 CALL histwrite(nid_nmc,"u600",itau_w, 52 .zx_tmp_2d,iim*jjmp1,ndex2d) 51 ENDDO !k=1, nlevSTD 52 ccc 53 c110304 END 53 54 c 54 CALL gr_fi_ecrit(1, klon,iim,jjmp1, v600,zx_tmp_2d)55 CALL histwrite(nid_nmc,"v600",itau_w,56 .zx_tmp_2d,iim*jjmp1,ndex2d)57 c58 CALL gr_fi_ecrit(1, klon,iim,jjmp1, u500,zx_tmp_2d)59 CALL histwrite(nid_nmc,"u500",itau_w,60 .zx_tmp_2d,iim*jjmp1,ndex2d)61 c62 CALL gr_fi_ecrit(1, klon,iim,jjmp1, v500,zx_tmp_2d)63 CALL histwrite(nid_nmc,"v500",itau_w,64 .zx_tmp_2d,iim*jjmp1,ndex2d)65 c66 CALL gr_fi_ecrit(1, klon,iim,jjmp1, u400,zx_tmp_2d)67 CALL histwrite(nid_nmc,"u400",itau_w,68 .zx_tmp_2d,iim*jjmp1,ndex2d)69 c70 CALL gr_fi_ecrit(1, klon,iim,jjmp1, v400,zx_tmp_2d)71 CALL histwrite(nid_nmc,"v400",itau_w,72 .zx_tmp_2d,iim*jjmp1,ndex2d)73 c74 CALL gr_fi_ecrit(1, klon,iim,jjmp1, u300,zx_tmp_2d)75 CALL histwrite(nid_nmc,"u300",itau_w,76 .zx_tmp_2d,iim*jjmp1,ndex2d)77 c78 CALL gr_fi_ecrit(1, klon,iim,jjmp1, v300,zx_tmp_2d)79 CALL histwrite(nid_nmc,"v300",itau_w,80 .zx_tmp_2d,iim*jjmp1,ndex2d)81 c82 CALL gr_fi_ecrit(1, klon,iim,jjmp1, u250,zx_tmp_2d)83 CALL histwrite(nid_nmc,"u250",itau_w,84 .zx_tmp_2d,iim*jjmp1,ndex2d)85 c86 CALL gr_fi_ecrit(1, klon,iim,jjmp1, v250,zx_tmp_2d)87 CALL histwrite(nid_nmc,"v250",itau_w,88 .zx_tmp_2d,iim*jjmp1,ndex2d)89 c90 CALL gr_fi_ecrit(1, klon,iim,jjmp1, u200,zx_tmp_2d)91 CALL histwrite(nid_nmc,"u200",itau_w,92 .zx_tmp_2d,iim*jjmp1,ndex2d)93 c94 CALL gr_fi_ecrit(1, klon,iim,jjmp1, v200,zx_tmp_2d)95 CALL histwrite(nid_nmc,"v200",itau_w,96 .zx_tmp_2d,iim*jjmp1,ndex2d)97 c98 CALL gr_fi_ecrit(1, klon,iim,jjmp1, u150,zx_tmp_2d)99 CALL histwrite(nid_nmc,"u150",itau_w,100 .zx_tmp_2d,iim*jjmp1,ndex2d)101 c102 CALL gr_fi_ecrit(1, klon,iim,jjmp1, v150,zx_tmp_2d)103 CALL histwrite(nid_nmc,"v150",itau_w,104 .zx_tmp_2d,iim*jjmp1,ndex2d)105 c106 CALL gr_fi_ecrit(1, klon,iim,jjmp1, u100,zx_tmp_2d)107 CALL histwrite(nid_nmc,"u100",itau_w,108 .zx_tmp_2d,iim*jjmp1,ndex2d)109 c110 CALL gr_fi_ecrit(1, klon,iim,jjmp1, v100,zx_tmp_2d)111 CALL histwrite(nid_nmc,"v100",itau_w,112 .zx_tmp_2d,iim*jjmp1,ndex2d)113 c114 CALL gr_fi_ecrit(1, klon,iim,jjmp1, u70,zx_tmp_2d)115 CALL histwrite(nid_nmc,"u70",itau_w,116 .zx_tmp_2d,iim*jjmp1,ndex2d)117 c118 CALL gr_fi_ecrit(1, klon,iim,jjmp1, v70,zx_tmp_2d)119 CALL histwrite(nid_nmc,"v70",itau_w,120 .zx_tmp_2d,iim*jjmp1,ndex2d)121 c122 CALL gr_fi_ecrit(1, klon,iim,jjmp1, u50,zx_tmp_2d)123 CALL histwrite(nid_nmc,"u50",itau_w,124 .zx_tmp_2d,iim*jjmp1,ndex2d)125 c126 CALL gr_fi_ecrit(1, klon,iim,jjmp1, v50,zx_tmp_2d)127 CALL histwrite(nid_nmc,"v50",itau_w,128 .zx_tmp_2d,iim*jjmp1,ndex2d)129 c130 CALL gr_fi_ecrit(1, klon,iim,jjmp1, u30,zx_tmp_2d)131 CALL histwrite(nid_nmc,"u30",itau_w,132 .zx_tmp_2d,iim*jjmp1,ndex2d)133 c134 CALL gr_fi_ecrit(1, klon,iim,jjmp1, v30,zx_tmp_2d)135 CALL histwrite(nid_nmc,"v30",itau_w,136 .zx_tmp_2d,iim*jjmp1,ndex2d)137 c138 CALL gr_fi_ecrit(1, klon,iim,jjmp1, u20,zx_tmp_2d)139 CALL histwrite(nid_nmc,"u20",itau_w,140 .zx_tmp_2d,iim*jjmp1,ndex2d)141 c142 CALL gr_fi_ecrit(1, klon,iim,jjmp1, v20,zx_tmp_2d)143 CALL histwrite(nid_nmc,"v20",itau_w,144 .zx_tmp_2d,iim*jjmp1,ndex2d)145 c146 CALL gr_fi_ecrit(1, klon,iim,jjmp1, u10,zx_tmp_2d)147 CALL histwrite(nid_nmc,"u10",itau_w,148 .zx_tmp_2d,iim*jjmp1,ndex2d)149 c150 CALL gr_fi_ecrit(1, klon,iim,jjmp1, v10,zx_tmp_2d)151 CALL histwrite(nid_nmc,"v10",itau_w,152 .zx_tmp_2d,iim*jjmp1,ndex2d)153 154 55 if (ok_sync) then 155 56 call histsync(nid_nmc)
Note: See TracChangeset
for help on using the changeset viewer.
