Changeset 1539 for LMDZ5/trunk/libf/phylmd/phys_output_write.h
- Timestamp:
- Jun 9, 2011, 12:13:33 AM (13 years ago)
- File:
-
- 1 edited
-
LMDZ5/trunk/libf/phylmd/phys_output_write.h (modified) (80 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/phylmd/phys_output_write.h
r1518 r1539 9 9 !!! Champs 1D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 10 10 IF (o_phis%flag(iff)<=lev_files(iff)) THEN 11 CALL histwrite_phy(nid_files(iff), 12 $ o_phis%name,itau_w,pphis) 13 ENDIF 14 11 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 12 $ o_phis%name,itau_w,pphis) 13 ENDIF 14 15 IF (.NOT.clef_stations(iff)) THEN 15 16 IF (o_aire%flag(iff)<=lev_files(iff)) THEN 16 CALL histwrite_phy(nid_files(iff),o_aire%name,itau_w,airephy) 17 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 18 $o_aire%name,itau_w,airephy) 17 19 ENDIF 18 20 … … 21 23 zx_tmp_fi2d(i)=pctsrf(i,is_ter)+pctsrf(i,is_lic) 22 24 ENDDO 23 CALL histwrite_phy(nid_files(iff), 25 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 24 26 $ o_contfracATM%name,itau_w,zx_tmp_fi2d) 25 27 ENDIF 28 ENDIF 26 29 27 30 IF (o_contfracOR%flag(iff)<=lev_files(iff)) THEN 28 CALL histwrite_phy(nid_files(iff), o_contfracOR%name,itau_w,29 $ 31 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 32 $ o_contfracOR%name,itau_w,pctsrf(:,is_ter)) 30 33 ENDIF 31 34 32 35 IF (o_aireTER%flag(iff)<=lev_files(iff)) THEN 33 CALL histwrite_phy(nid_files(iff), 36 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 34 37 $ o_aireTER%name,itau_w,paire_ter) 35 38 ENDIF … … 38 41 39 42 IF (o_flat%flag(iff)<=lev_files(iff)) THEN 40 CALL histwrite_phy(nid_files(iff),o_flat%name,itau_w,zxfluxlat) 43 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 44 $ o_flat%name,itau_w,zxfluxlat) 41 45 ENDIF 42 46 43 47 IF (o_slp%flag(iff)<=lev_files(iff)) THEN 44 CALL histwrite_phy(nid_files(iff),o_slp%name,itau_w,slp) 48 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 49 $ o_slp%name,itau_w,slp) 45 50 ENDIF 46 51 47 52 IF (o_tsol%flag(iff)<=lev_files(iff)) THEN 48 CALL histwrite_phy(nid_files(iff),o_tsol%name,itau_w,zxtsol) 53 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 54 $ o_tsol%name,itau_w,zxtsol) 49 55 ENDIF 50 56 51 57 IF (o_t2m%flag(iff)<=lev_files(iff)) THEN 52 CALL histwrite_phy(nid_files(iff),o_t2m%name,itau_w,zt2m) 53 ENDIF 54 58 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 59 $ o_t2m%name,itau_w,zt2m) 60 ENDIF 61 62 IF (.NOT.clef_stations(iff)) THEN 55 63 IF (o_t2m_min%flag(iff)<=lev_files(iff)) THEN 56 CALL histwrite_phy(nid_files(iff),o_t2m_min%name,itau_w,zt2m) 64 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 65 $o_t2m_min%name,itau_w,zt2m) 57 66 ENDIF 58 67 59 68 IF (o_t2m_max%flag(iff)<=lev_files(iff)) THEN 60 CALL histwrite_phy(nid_files(iff),o_t2m_max%name,itau_w,zt2m) 69 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 70 $o_t2m_max%name,itau_w,zt2m) 71 ENDIF 61 72 ENDIF 62 73 … … 65 76 zx_tmp_fi2d(i)=SQRT(zu10m(i)*zu10m(i)+zv10m(i)*zv10m(i)) 66 77 ENDDO 67 CALL histwrite_phy(nid_files(iff), 78 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 68 79 s o_wind10m%name,itau_w,zx_tmp_fi2d) 69 80 ENDIF 70 81 82 IF (.NOT.clef_stations(iff)) THEN 71 83 IF (o_wind10max%flag(iff)<=lev_files(iff)) THEN 72 84 DO i=1, klon 73 85 zx_tmp_fi2d(i)=SQRT(zu10m(i)*zu10m(i)+zv10m(i)*zv10m(i)) 74 86 ENDDO 75 CALL histwrite_phy(nid_files(iff),o_wind10max%name, 76 $ itau_w,zx_tmp_fi2d) 77 ENDIF 87 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 88 $o_wind10max%name,itau_w,zx_tmp_fi2d) 89 ENDIF 90 ENDIF 78 91 79 92 IF (o_sicf%flag(iff)<=lev_files(iff)) THEN … … 81 94 zx_tmp_fi2d(i) = pctsrf(i,is_sic) 82 95 ENDDO 83 CALL histwrite_phy(nid_files(iff), 96 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 84 97 $ o_sicf%name,itau_w,zx_tmp_fi2d) 85 98 ENDIF 86 99 87 100 IF (o_q2m%flag(iff)<=lev_files(iff)) THEN 88 CALL histwrite_phy(nid_files(iff),o_q2m%name,itau_w,zq2m) 101 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 102 $o_q2m%name,itau_w,zq2m) 89 103 ENDIF 90 104 91 105 IF (o_u10m%flag(iff)<=lev_files(iff)) THEN 92 CALL histwrite_phy(nid_files(iff),o_u10m%name,itau_w,zu10m) 106 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 107 $o_u10m%name,itau_w,zu10m) 93 108 ENDIF 94 109 95 110 IF (o_v10m%flag(iff)<=lev_files(iff)) THEN 96 CALL histwrite_phy(nid_files(iff),o_v10m%name,itau_w,zv10m) 111 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 112 $o_v10m%name,itau_w,zv10m) 97 113 ENDIF 98 114 … … 101 117 zx_tmp_fi2d(i) = paprs(i,1) 102 118 ENDDO 103 CALL histwrite_phy(nid_files(iff), 119 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 104 120 s o_psol%name,itau_w,zx_tmp_fi2d) 105 121 ENDIF 106 122 107 123 IF (o_mass%flag(iff)<=lev_files(iff)) THEN 108 CALL histwrite_phy(nid_files(iff),o_mass%name,itau_w,zmasse) 124 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 125 s o_mass%name,itau_w,zmasse) 109 126 ENDIF 110 127 111 128 112 129 IF (o_qsurf%flag(iff)<=lev_files(iff)) THEN 113 CALL histwrite_phy(nid_files(iff),o_qsurf%name,itau_w,zxqsurf) 130 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 131 $o_qsurf%name,itau_w,zxqsurf) 114 132 ENDIF 115 133 116 134 if (.not. ok_veget) then 117 135 IF (o_qsol%flag(iff)<=lev_files(iff)) THEN 118 CALL histwrite_phy(nid_files(iff),o_qsol%name,itau_w,qsol) 136 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 137 $o_qsol%name,itau_w,qsol) 119 138 ENDIF 120 139 endif … … 124 143 zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i) 125 144 ENDDO 126 CALL histwrite_phy(nid_files(iff), o_precip%name,127 sitau_w,zx_tmp_fi2d)145 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 146 $o_precip%name,itau_w,zx_tmp_fi2d) 128 147 ENDIF 129 148 130 149 IF (o_ndayrain%flag(iff)<=lev_files(iff)) THEN 131 CALL histwrite_phy(nid_files(iff), o_ndayrain%name,132 sitau_w,nday_rain)150 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 151 $o_ndayrain%name,itau_w,nday_rain) 133 152 ENDIF 134 153 … … 137 156 zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i) 138 157 ENDDO 139 CALL histwrite_phy(nid_files(iff),o_plul%name,itau_w,zx_tmp_fi2d) 158 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 159 $o_plul%name,itau_w,zx_tmp_fi2d) 140 160 ENDIF 141 161 … … 144 164 zx_tmp_fi2d(i) = rain_con(i) + snow_con(i) 145 165 ENDDO 146 CALL histwrite_phy(nid_files(iff),o_pluc%name,itau_w,zx_tmp_fi2d) 166 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 167 $o_pluc%name,itau_w,zx_tmp_fi2d) 147 168 ENDIF 148 169 149 170 IF (o_snow%flag(iff)<=lev_files(iff)) THEN 150 CALL histwrite_phy(nid_files(iff),o_snow%name,itau_w,snow_fall) 171 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 172 $o_snow%name,itau_w,snow_fall) 151 173 ENDIF 152 174 153 175 IF (o_msnow%flag(iff)<=lev_files(iff)) THEN 154 CALL histwrite_phy(nid_files(iff),o_msnow%name,itau_w,snow_o) 176 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 177 $o_msnow%name,itau_w,snow_o) 155 178 ENDIF 156 179 157 180 IF (o_fsnow%flag(iff)<=lev_files(iff)) THEN 158 CALL histwrite_phy(nid_files(iff),o_fsnow%name,itau_w,zfra_o) 181 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 182 $o_fsnow%name,itau_w,zfra_o) 159 183 ENDIF 160 184 161 185 IF (o_evap%flag(iff)<=lev_files(iff)) THEN 162 CALL histwrite_phy(nid_files(iff),o_evap%name,itau_w,evap) 186 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 187 $o_evap%name,itau_w,evap) 163 188 ENDIF 164 189 165 190 IF (o_tops%flag(iff)<=lev_files(iff)) THEN 166 CALL histwrite_phy(nid_files(iff),o_tops%name,itau_w,topsw) 191 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 192 $o_tops%name,itau_w,topsw) 167 193 ENDIF 168 194 169 195 IF (o_tops0%flag(iff)<=lev_files(iff)) THEN 170 CALL histwrite_phy(nid_files(iff),o_tops0%name,itau_w,topsw0) 196 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 197 $o_tops0%name,itau_w,topsw0) 171 198 ENDIF 172 199 173 200 IF (o_topl%flag(iff)<=lev_files(iff)) THEN 174 CALL histwrite_phy(nid_files(iff),o_topl%name,itau_w,toplw) 201 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 202 $o_topl%name,itau_w,toplw) 175 203 ENDIF 176 204 177 205 IF (o_topl0%flag(iff)<=lev_files(iff)) THEN 178 CALL histwrite_phy(nid_files(iff),o_topl0%name,itau_w,toplw0) 206 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 207 $o_topl0%name,itau_w,toplw0) 179 208 ENDIF 180 209 181 210 IF (o_SWupTOA%flag(iff)<=lev_files(iff)) THEN 182 211 zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, klevp1 ) 183 CALL histwrite_phy(nid_files(iff), o_SWupTOA%name,184 sitau_w,zx_tmp_fi2d)212 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 213 $o_SWupTOA%name,itau_w,zx_tmp_fi2d) 185 214 ENDIF 186 215 187 216 IF (o_SWupTOAclr%flag(iff)<=lev_files(iff)) THEN 188 217 zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, klevp1 ) 189 CALL histwrite_phy(nid_files(iff), 218 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 190 219 $ o_SWupTOAclr%name,itau_w,zx_tmp_fi2d) 191 220 ENDIF … … 193 222 IF (o_SWdnTOA%flag(iff)<=lev_files(iff)) THEN 194 223 zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, klevp1 ) 195 CALL histwrite_phy(nid_files(iff), 224 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 196 225 s o_SWdnTOA%name,itau_w,zx_tmp_fi2d) 197 226 ENDIF … … 199 228 IF (o_SWdnTOAclr%flag(iff)<=lev_files(iff)) THEN 200 229 zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, klevp1 ) 201 CALL histwrite_phy(nid_files(iff), 230 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 202 231 $ o_SWdnTOAclr%name,itau_w,zx_tmp_fi2d) 203 232 ENDIF … … 205 234 IF (o_nettop%flag(iff)<=lev_files(iff)) THEN 206 235 zx_tmp_fi2d(:) = topsw(:)-toplw(:) 207 CALL histwrite_phy(nid_files(iff), 236 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 208 237 $ o_nettop%name,itau_w,zx_tmp_fi2d) 209 238 ENDIF 210 239 211 240 IF (o_SWup200%flag(iff)<=lev_files(iff)) THEN 212 CALL histwrite_phy(nid_files(iff),o_SWup200%name,itau_w,SWup200) 241 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 242 $o_SWup200%name,itau_w,SWup200) 213 243 ENDIF 214 244 215 245 IF (o_SWup200clr%flag(iff)<=lev_files(iff)) THEN 216 CALL histwrite_phy(nid_files(iff), 246 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 217 247 s o_SWup200clr%name,itau_w,SWup200clr) 218 248 ENDIF 219 249 220 250 IF (o_SWdn200%flag(iff)<=lev_files(iff)) THEN 221 CALL histwrite_phy(nid_files(iff),o_SWdn200%name,itau_w,SWdn200) 251 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 252 $o_SWdn200%name,itau_w,SWdn200) 222 253 ENDIF 223 254 224 255 IF (o_SWdn200clr%flag(iff)<=lev_files(iff)) THEN 225 CALL histwrite_phy(nid_files(iff), 256 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 226 257 s o_SWdn200clr%name,itau_w,SWdn200clr) 227 258 ENDIF 228 259 229 260 IF (o_LWup200%flag(iff)<=lev_files(iff)) THEN 230 CALL histwrite_phy(nid_files(iff),o_LWup200%name,itau_w,LWup200) 261 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 262 $o_LWup200%name,itau_w,LWup200) 231 263 ENDIF 232 264 233 265 IF (o_LWup200clr%flag(iff)<=lev_files(iff)) THEN 234 CALL histwrite_phy(nid_files(iff), 266 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 235 267 s o_LWup200clr%name,itau_w,LWup200clr) 236 268 ENDIF 237 269 238 270 IF (o_LWdn200%flag(iff)<=lev_files(iff)) THEN 239 CALL histwrite_phy(nid_files(iff), 271 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 240 272 s o_LWdn200%name,itau_w,LWdn200) 241 273 ENDIF 242 274 243 275 IF (o_LWdn200clr%flag(iff)<=lev_files(iff)) THEN 244 CALL histwrite_phy(nid_files(iff), 276 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 245 277 s o_LWdn200clr%name,itau_w,LWdn200clr) 246 278 ENDIF 247 279 248 280 IF (o_sols%flag(iff)<=lev_files(iff)) THEN 249 CALL histwrite_phy(nid_files(iff),o_sols%name,itau_w,solsw) 281 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 282 $o_sols%name,itau_w,solsw) 250 283 ENDIF 251 284 252 285 IF (o_sols0%flag(iff)<=lev_files(iff)) THEN 253 CALL histwrite_phy(nid_files(iff),o_sols0%name,itau_w,solsw0) 286 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 287 $o_sols0%name,itau_w,solsw0) 254 288 ENDIF 255 289 256 290 IF (o_soll%flag(iff)<=lev_files(iff)) THEN 257 CALL histwrite_phy(nid_files(iff),o_soll%name,itau_w,sollw) 291 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 292 $o_soll%name,itau_w,sollw) 258 293 ENDIF 259 294 260 295 IF (o_radsol%flag(iff)<=lev_files(iff)) THEN 261 CALL histwrite_phy(nid_files(iff),o_radsol%name,itau_w,radsol) 296 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 297 $o_radsol%name,itau_w,radsol) 262 298 ENDIF 263 299 264 300 IF (o_soll0%flag(iff)<=lev_files(iff)) THEN 265 CALL histwrite_phy(nid_files(iff),o_soll0%name,itau_w,sollw0) 301 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 302 $o_soll0%name,itau_w,sollw0) 266 303 ENDIF 267 304 268 305 IF (o_SWupSFC%flag(iff)<=lev_files(iff)) THEN 269 306 zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, 1 ) 270 CALL histwrite_phy(nid_files(iff), 307 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 271 308 s o_SWupSFC%name,itau_w,zx_tmp_fi2d) 272 309 ENDIF … … 274 311 IF (o_SWupSFCclr%flag(iff)<=lev_files(iff)) THEN 275 312 zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, 1 ) 276 CALL histwrite_phy(nid_files(iff), 313 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 277 314 $ o_SWupSFCclr%name,itau_w,zx_tmp_fi2d) 278 315 ENDIF … … 280 317 IF (o_SWdnSFC%flag(iff)<=lev_files(iff)) THEN 281 318 zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, 1 ) 282 CALL histwrite_phy(nid_files(iff), 319 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 283 320 $ o_SWdnSFC%name,itau_w,zx_tmp_fi2d) 284 321 ENDIF … … 286 323 IF (o_SWdnSFCclr%flag(iff)<=lev_files(iff)) THEN 287 324 zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, 1 ) 288 CALL histwrite_phy(nid_files(iff), 325 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 289 326 $ o_SWdnSFCclr%name,itau_w,zx_tmp_fi2d) 290 327 ENDIF … … 292 329 IF (o_LWupSFC%flag(iff)<=lev_files(iff)) THEN 293 330 zx_tmp_fi2d(1:klon)=sollwdown(1:klon)-sollw(1:klon) 294 CALL histwrite_phy(nid_files(iff), 331 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 295 332 $ o_LWupSFC%name,itau_w,zx_tmp_fi2d) 296 333 ENDIF 297 334 298 335 IF (o_LWdnSFC%flag(iff)<=lev_files(iff)) THEN 299 CALL histwrite_phy(nid_files(iff), 336 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 300 337 $ o_LWdnSFC%name,itau_w,sollwdown) 301 338 ENDIF … … 304 341 IF (o_LWupSFCclr%flag(iff)<=lev_files(iff)) THEN 305 342 zx_tmp_fi2d(1:klon)=sollwdownclr(1:klon)-sollw0(1:klon) 306 CALL histwrite_phy(nid_files(iff), 343 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 307 344 $ o_LWupSFCclr%name,itau_w,zx_tmp_fi2d) 308 345 ENDIF 309 346 310 347 IF (o_LWdnSFCclr%flag(iff)<=lev_files(iff)) THEN 311 CALL histwrite_phy(nid_files(iff), 348 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 312 349 $ o_LWdnSFCclr%name,itau_w,sollwdownclr) 313 350 ENDIF 314 351 315 352 IF (o_bils%flag(iff)<=lev_files(iff)) THEN 316 CALL histwrite_phy(nid_files(iff),o_bils%name,itau_w,bils) 353 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 354 $o_bils%name,itau_w,bils) 317 355 ENDIF 318 356 319 357 IF (o_sens%flag(iff)<=lev_files(iff)) THEN 320 358 zx_tmp_fi2d(1:klon)=-1*sens(1:klon) 321 CALL histwrite_phy(nid_files(iff),o_sens%name,itau_w,zx_tmp_fi2d) 359 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 360 $o_sens%name,itau_w,zx_tmp_fi2d) 322 361 ENDIF 323 362 324 363 IF (o_fder%flag(iff)<=lev_files(iff)) THEN 325 CALL histwrite_phy(nid_files(iff),o_fder%name,itau_w,fder) 364 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 365 $o_fder%name,itau_w,fder) 326 366 ENDIF 327 367 328 368 IF (o_ffonte%flag(iff)<=lev_files(iff)) THEN 329 CALL histwrite_phy(nid_files(iff),o_ffonte%name,itau_w,zxffonte) 369 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 370 $o_ffonte%name,itau_w,zxffonte) 330 371 ENDIF 331 372 332 373 IF (o_fqcalving%flag(iff)<=lev_files(iff)) THEN 333 CALL histwrite_phy(nid_files(iff), 374 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 334 375 $ o_fqcalving%name,itau_w,zxfqcalving) 335 376 ENDIF 336 377 337 378 IF (o_fqfonte%flag(iff)<=lev_files(iff)) THEN 338 CALL histwrite_phy(nid_files(iff), 379 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 339 380 $ o_fqfonte%name,itau_w,zxfqfonte) 340 381 ENDIF … … 345 386 zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+pctsrf(:,nsrf)*fluxu(:,1,nsrf) 346 387 enddo 347 CALL histwrite_phy(nid_files(iff), 388 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 348 389 $ o_taux%name,itau_w,zx_tmp_fi2d) 349 390 ENDIF … … 354 395 zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+pctsrf(:,nsrf)*fluxv(:,1,nsrf) 355 396 enddo 356 CALL histwrite_phy(nid_files(iff), 397 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 357 398 $ o_tauy%name,itau_w,zx_tmp_fi2d) 358 399 ENDIF … … 363 404 zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)*100. 364 405 CALL histwrite_phy(nid_files(iff), 406 $ clef_stations(iff), 365 407 $ o_pourc_srf(nsrf)%name,itau_w, 366 408 $ zx_tmp_fi2d) … … 369 411 IF (o_fract_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN 370 412 zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf) 371 CALL histwrite_phy(nid_files(iff), 413 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 372 414 $ o_fract_srf(nsrf)%name,itau_w, 373 415 $ zx_tmp_fi2d) … … 377 419 IF (o_taux_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN 378 420 zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf) 379 CALL histwrite_phy(nid_files(iff), 421 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 380 422 $ o_taux_srf(nsrf)%name,itau_w, 381 423 $ zx_tmp_fi2d) … … 384 426 IF (o_tauy_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN 385 427 zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf) 386 CALL histwrite_phy(nid_files(iff), 428 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 387 429 $ o_tauy_srf(nsrf)%name,itau_w, 388 430 $ zx_tmp_fi2d) … … 391 433 IF (o_tsol_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN 392 434 zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf) 393 CALL histwrite_phy(nid_files(iff), 435 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 394 436 $ o_tsol_srf(nsrf)%name,itau_w, 395 437 $ zx_tmp_fi2d) … … 398 440 IF (o_u10m_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN 399 441 zx_tmp_fi2d(1 : klon) = u10m(1 : klon, nsrf) 400 CALL histwrite_phy(nid_files(iff),o_u10m_srf(nsrf)%name, 442 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 443 $o_u10m_srf(nsrf)%name, 401 444 $ itau_w,zx_tmp_fi2d) 402 445 ENDIF … … 404 447 IF (o_v10m_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN 405 448 zx_tmp_fi2d(1 : klon) = v10m(1 : klon, nsrf) 406 CALL histwrite_phy(nid_files(iff),o_v10m_srf(nsrf)%name, 449 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 450 $o_v10m_srf(nsrf)%name, 407 451 $ itau_w,zx_tmp_fi2d) 408 452 ENDIF … … 410 454 IF (o_t2m_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN 411 455 zx_tmp_fi2d(1 : klon) = t2m(1 : klon, nsrf) 412 CALL histwrite_phy(nid_files(iff),o_t2m_srf(nsrf)%name, 456 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 457 $o_t2m_srf(nsrf)%name, 413 458 $ itau_w,zx_tmp_fi2d) 414 459 ENDIF … … 416 461 IF (o_evap_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN 417 462 zx_tmp_fi2d(1 : klon) = fevap(1 : klon, nsrf) 418 CALL histwrite_phy(nid_files(iff),o_evap_srf(nsrf)%name, 463 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 464 $o_evap_srf(nsrf)%name, 419 465 $ itau_w,zx_tmp_fi2d) 420 466 ENDIF … … 422 468 IF (o_sens_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN 423 469 zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf) 424 CALL histwrite_phy(nid_files(iff), 470 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 425 471 $ o_sens_srf(nsrf)%name,itau_w, 426 472 $ zx_tmp_fi2d) … … 429 475 IF (o_lat_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN 430 476 zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf) 431 CALL histwrite_phy(nid_files(iff), 477 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 432 478 $ o_lat_srf(nsrf)%name,itau_w, 433 479 $ zx_tmp_fi2d) … … 436 482 IF (o_flw_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN 437 483 zx_tmp_fi2d(1 : klon) = fsollw( 1 : klon, nsrf) 438 CALL histwrite_phy(nid_files(iff), 484 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 439 485 $ o_flw_srf(nsrf)%name,itau_w, 440 486 $ zx_tmp_fi2d) … … 443 489 IF (o_fsw_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN 444 490 zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, nsrf) 445 CALL histwrite_phy(nid_files(iff), 491 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 446 492 $ o_fsw_srf(nsrf)%name,itau_w, 447 493 $ zx_tmp_fi2d) … … 450 496 IF (o_wbils_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN 451 497 zx_tmp_fi2d(1 : klon) = wfbils( 1 : klon, nsrf) 452 CALL histwrite_phy(nid_files(iff), 498 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 453 499 $ o_wbils_srf(nsrf)%name,itau_w, 454 500 $ zx_tmp_fi2d) … … 457 503 IF (o_wbilo_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN 458 504 zx_tmp_fi2d(1 : klon) = wfbilo( 1 : klon, nsrf) 459 CALL histwrite_phy(nid_files(iff), 505 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 460 506 $ o_wbilo_srf(nsrf)%name,itau_w, 461 507 $ zx_tmp_fi2d) … … 464 510 if (iflag_pbl>1 .and. lev_histday.gt.10 ) then 465 511 IF (o_tke_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN 466 CALL histwrite_phy(nid_files(iff), 512 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 467 513 $ o_tke_srf(nsrf)%name,itau_w, 468 514 $ pbl_tke(:,1:klev,nsrf)) 469 515 ENDIF 470 516 517 IF (.NOT.clef_stations(iff)) THEN 471 518 IF (o_tke_max_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN 472 CALL histwrite_phy(nid_files(iff), 519 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 473 520 $ o_tke_max_srf(nsrf)%name,itau_w, 474 521 $ pbl_tke(:,1:klev,nsrf)) 475 522 ENDIF 523 ENDIF 476 524 endif 477 525 ENDDO 478 526 479 527 IF (o_cdrm%flag(iff)<=lev_files(iff)) THEN 480 CALL histwrite_phy(nid_files(iff),o_cdrm%name,itau_w,cdragm) 528 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 529 $o_cdrm%name,itau_w,cdragm) 481 530 ENDIF 482 531 483 532 IF (o_cdrh%flag(iff)<=lev_files(iff)) THEN 484 CALL histwrite_phy(nid_files(iff),o_cdrh%name,itau_w,cdragh) 533 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 534 $o_cdrh%name,itau_w,cdragh) 485 535 ENDIF 486 536 487 537 IF (o_cldl%flag(iff)<=lev_files(iff)) THEN 488 CALL histwrite_phy(nid_files(iff),o_cldl%name,itau_w,cldl) 538 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 539 $o_cldl%name,itau_w,cldl) 489 540 ENDIF 490 541 491 542 IF (o_cldm%flag(iff)<=lev_files(iff)) THEN 492 CALL histwrite_phy(nid_files(iff),o_cldm%name,itau_w,cldm) 543 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 544 $o_cldm%name,itau_w,cldm) 493 545 ENDIF 494 546 495 547 IF (o_cldh%flag(iff)<=lev_files(iff)) THEN 496 CALL histwrite_phy(nid_files(iff),o_cldh%name,itau_w,cldh) 548 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 549 $o_cldh%name,itau_w,cldh) 497 550 ENDIF 498 551 499 552 IF (o_cldt%flag(iff)<=lev_files(iff)) THEN 500 CALL histwrite_phy(nid_files(iff),o_cldt%name, 553 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 554 $o_cldt%name, 501 555 & itau_w,cldt) 502 556 ENDIF 503 557 504 558 IF (o_cldq%flag(iff)<=lev_files(iff)) THEN 505 CALL histwrite_phy(nid_files(iff),o_cldq%name,itau_w,cldq) 559 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 560 $o_cldq%name,itau_w,cldq) 506 561 ENDIF 507 562 508 563 IF (o_lwp%flag(iff)<=lev_files(iff)) THEN 509 564 zx_tmp_fi2d(1:klon) = flwp(1:klon) 510 CALL histwrite_phy(nid_files(iff), 565 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 511 566 s o_lwp%name,itau_w,zx_tmp_fi2d) 512 567 ENDIF … … 514 569 IF (o_iwp%flag(iff)<=lev_files(iff)) THEN 515 570 zx_tmp_fi2d(1:klon) = fiwp(1:klon) 516 CALL histwrite_phy(nid_files(iff), 571 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 517 572 s o_iwp%name,itau_w,zx_tmp_fi2d) 518 573 ENDIF 519 574 520 575 IF (o_ue%flag(iff)<=lev_files(iff)) THEN 521 CALL histwrite_phy(nid_files(iff),o_ue%name,itau_w,ue) 576 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 577 $o_ue%name,itau_w,ue) 522 578 ENDIF 523 579 524 580 IF (o_ve%flag(iff)<=lev_files(iff)) THEN 525 CALL histwrite_phy(nid_files(iff),o_ve%name,itau_w,ve) 581 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 582 $o_ve%name,itau_w,ve) 526 583 ENDIF 527 584 528 585 IF (o_uq%flag(iff)<=lev_files(iff)) THEN 529 CALL histwrite_phy(nid_files(iff),o_uq%name,itau_w,uq) 586 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 587 $o_uq%name,itau_w,uq) 530 588 ENDIF 531 589 532 590 IF (o_vq%flag(iff)<=lev_files(iff)) THEN 533 CALL histwrite_phy(nid_files(iff),o_vq%name,itau_w,vq) 591 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 592 $o_vq%name,itau_w,vq) 534 593 ENDIF 535 594 536 595 IF(iflag_con.GE.3) THEN ! sb 537 596 IF (o_cape%flag(iff)<=lev_files(iff)) THEN 538 CALL histwrite_phy(nid_files(iff),o_cape%name,itau_w,cape) 597 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 598 $o_cape%name,itau_w,cape) 539 599 ENDIF 540 600 541 601 IF (o_pbase%flag(iff)<=lev_files(iff)) THEN 542 CALL histwrite_phy(nid_files(iff),o_pbase%name,itau_w,ema_pcb) 602 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 603 $o_pbase%name,itau_w,ema_pcb) 543 604 ENDIF 544 605 545 606 IF (o_ptop%flag(iff)<=lev_files(iff)) THEN 546 CALL histwrite_phy(nid_files(iff),o_ptop%name,itau_w,ema_pct) 607 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 608 $o_ptop%name,itau_w,ema_pct) 547 609 ENDIF 548 610 549 611 IF (o_fbase%flag(iff)<=lev_files(iff)) THEN 550 CALL histwrite_phy(nid_files(iff),o_fbase%name,itau_w,ema_cbmf) 612 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 613 $o_fbase%name,itau_w,ema_cbmf) 551 614 ENDIF 552 615 553 616 IF (o_plcl%flag(iff)<=lev_files(iff)) THEN 554 CALL histwrite_phy(nid_files(iff),o_plcl%name,itau_w,plcl) 617 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 618 $o_plcl%name,itau_w,plcl) 555 619 ENDIF 556 620 557 621 IF (o_plfc%flag(iff)<=lev_files(iff)) THEN 558 CALL histwrite_phy(nid_files(iff),o_plfc%name,itau_w,plfc) 622 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 623 $o_plfc%name,itau_w,plfc) 559 624 ENDIF 560 625 561 626 IF (o_wbeff%flag(iff)<=lev_files(iff)) THEN 562 CALL histwrite_phy(nid_files(iff),o_wbeff%name,itau_w,wbeff) 627 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 628 $o_wbeff%name,itau_w,wbeff) 563 629 ENDIF 564 630 565 631 566 632 IF (o_prw%flag(iff)<=lev_files(iff)) THEN 567 CALL histwrite_phy(nid_files(iff),o_prw%name,itau_w,prw) 568 ENDIF 569 633 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 634 $o_prw%name,itau_w,prw) 635 ENDIF 636 637 IF (.NOT.clef_stations(iff)) THEN 570 638 IF (o_cape_max%flag(iff)<=lev_files(iff)) THEN 571 CALL histwrite_phy(nid_files(iff),o_cape_max%name,itau_w,cape) 639 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 640 $o_cape_max%name,itau_w,cape) 641 ENDIF 572 642 ENDIF 573 643 574 644 IF (o_upwd%flag(iff)<=lev_files(iff)) THEN 575 CALL histwrite_phy(nid_files(iff),o_upwd%name,itau_w,upwd) 645 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 646 $o_upwd%name,itau_w,upwd) 576 647 ENDIF 577 648 578 649 IF (o_Ma%flag(iff)<=lev_files(iff)) THEN 579 CALL histwrite_phy(nid_files(iff),o_Ma%name,itau_w,Ma) 650 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 651 $o_Ma%name,itau_w,Ma) 580 652 ENDIF 581 653 582 654 IF (o_dnwd%flag(iff)<=lev_files(iff)) THEN 583 CALL histwrite_phy(nid_files(iff),o_dnwd%name,itau_w,dnwd) 655 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 656 $o_dnwd%name,itau_w,dnwd) 584 657 ENDIF 585 658 586 659 IF (o_dnwd0%flag(iff)<=lev_files(iff)) THEN 587 CALL histwrite_phy(nid_files(iff),o_dnwd0%name,itau_w,dnwd0) 660 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 661 $o_dnwd0%name,itau_w,dnwd0) 588 662 ENDIF 589 663 590 664 IF (o_ftime_con%flag(iff)<=lev_files(iff)) THEN 591 665 zx_tmp_fi2d=float(itau_con)/float(itap) 592 CALL histwrite_phy(nid_files(iff),o_ftime_con%name, 666 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 667 $o_ftime_con%name, 593 668 s itau_w,zx_tmp_fi2d) 594 669 ENDIF … … 600 675 zx_tmp_fi3d=dnwd+dnwd0+upwd 601 676 endif 602 CALL histwrite_phy(nid_files(iff),o_mc%name,itau_w,zx_tmp_fi3d) 677 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 678 $o_mc%name,itau_w,zx_tmp_fi3d) 603 679 ENDIF 604 680 … … 606 682 607 683 IF (o_s_pblh%flag(iff)<=lev_files(iff)) THEN 608 CALL histwrite_phy(nid_files(iff),o_s_pblh%name,itau_w,s_pblh) 684 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 685 $o_s_pblh%name,itau_w,s_pblh) 609 686 ENDIF 610 687 611 688 IF (o_s_pblt%flag(iff)<=lev_files(iff)) THEN 612 CALL histwrite_phy(nid_files(iff),o_s_pblt%name,itau_w,s_pblt) 689 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 690 $o_s_pblt%name,itau_w,s_pblt) 613 691 ENDIF 614 692 615 693 IF (o_s_lcl%flag(iff)<=lev_files(iff)) THEN 616 CALL histwrite_phy(nid_files(iff),o_s_lcl%name,itau_w,s_lcl) 694 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 695 $o_s_lcl%name,itau_w,s_lcl) 617 696 ENDIF 618 697 619 698 IF (o_s_therm%flag(iff)<=lev_files(iff)) THEN 620 CALL histwrite_phy(nid_files(iff),o_s_therm%name,itau_w,s_therm) 699 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 700 $o_s_therm%name,itau_w,s_therm) 621 701 ENDIF 622 702 623 703 !IM : Les champs suivants (s_capCL, s_oliqCL, s_cteiCL, s_trmb1, s_trmb2, s_trmb3) ne sont pas definis dans HBTM.F 624 704 ! IF (o_s_capCL%flag(iff)<=lev_files(iff)) THEN 625 ! CALL histwrite_phy(nid_files(iff),o_s_capCL%name,itau_w,s_capCL) 705 ! CALL histwrite_phy(nid_files(iff),clef_stations(iff), 706 ! $o_s_capCL%name,itau_w,s_capCL) 626 707 ! ENDIF 627 708 628 709 ! IF (o_s_oliqCL%flag(iff)<=lev_files(iff)) THEN 629 ! CALL histwrite_phy(nid_files(iff),o_s_oliqCL%name,itau_w,s_oliqCL) 710 ! CALL histwrite_phy(nid_files(iff),clef_stations(iff), 711 ! $o_s_oliqCL%name,itau_w,s_oliqCL) 630 712 ! ENDIF 631 713 632 714 ! IF (o_s_cteiCL%flag(iff)<=lev_files(iff)) THEN 633 ! CALL histwrite_phy(nid_files(iff),o_s_cteiCL%name,itau_w,s_cteiCL) 715 ! CALL histwrite_phy(nid_files(iff),clef_stations(iff), 716 ! $o_s_cteiCL%name,itau_w,s_cteiCL) 634 717 ! ENDIF 635 718 636 719 ! IF (o_s_trmb1%flag(iff)<=lev_files(iff)) THEN 637 ! CALL histwrite_phy(nid_files(iff),o_s_trmb1%name,itau_w,s_trmb1) 720 ! CALL histwrite_phy(nid_files(iff),clef_stations(iff), 721 ! $o_s_trmb1%name,itau_w,s_trmb1) 638 722 ! ENDIF 639 723 640 724 ! IF (o_s_trmb2%flag(iff)<=lev_files(iff)) THEN 641 ! CALL histwrite_phy(nid_files(iff),o_s_trmb2%name,itau_w,s_trmb2) 725 ! CALL histwrite_phy(nid_files(iff),clef_stations(iff), 726 ! $o_s_trmb2%name,itau_w,s_trmb2) 642 727 ! ENDIF 643 728 644 729 ! IF (o_s_trmb3%flag(iff)<=lev_files(iff)) THEN 645 ! CALL histwrite_phy(nid_files(iff),o_s_trmb3%name,itau_w,s_trmb3) 730 ! CALL histwrite_phy(nid_files(iff),clef_stations(iff), 731 ! $o_s_trmb3%name,itau_w,s_trmb3) 646 732 ! ENDIF 647 733 … … 650 736 ll=0 651 737 DO k=1, nlevSTD 652 ! IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)653 ! IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)654 738 bb2=clevSTD(k) 655 739 IF(bb2.EQ."850".OR.bb2.EQ."700".OR. 656 740 $ bb2.EQ."500".OR.bb2.EQ."200".OR. 741 $ bb2.EQ."100".OR. 657 742 $ bb2.EQ."50".OR.bb2.EQ."10") THEN 658 743 … … 660 745 ll=ll+1 661 746 IF (o_uSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN 662 CALL histwrite_phy(nid_files(iff),o_uSTDlevs(ll)%name, 747 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 748 $o_uSTDlevs(ll)%name, 663 749 & itau_w,uwriteSTD(:,k,iff)) 664 750 ENDIF 665 751 666 752 IF (o_vSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN 667 CALL histwrite_phy(nid_files(iff),o_vSTDlevs(ll)%name, 753 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 754 $o_vSTDlevs(ll)%name, 668 755 & itau_w,vwriteSTD(:,k,iff)) 669 756 ENDIF 670 757 671 758 IF (o_wSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN 672 CALL histwrite_phy(nid_files(iff),o_wSTDlevs(ll)%name, 759 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 760 $o_wSTDlevs(ll)%name, 673 761 & itau_w,wwriteSTD(:,k,iff)) 674 762 ENDIF 675 763 676 764 IF (o_zSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN 677 CALL histwrite_phy(nid_files(iff),o_zSTDlevs(ll)%name, 765 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 766 $o_zSTDlevs(ll)%name, 678 767 & itau_w,phiwriteSTD(:,k,iff)) 679 768 ENDIF 680 769 681 770 IF (o_qSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN 682 CALL histwrite_phy(nid_files(iff),o_qSTDlevs(ll)%name, 771 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 772 $o_qSTDlevs(ll)%name, 683 773 & itau_w, qwriteSTD(:,k,iff)) 684 774 ENDIF 685 775 686 776 IF (o_tSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN 687 CALL histwrite_phy(nid_files(iff),o_tSTDlevs(ll)%name, 777 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 778 $o_tSTDlevs(ll)%name, 688 779 & itau_w, twriteSTD(:,k,iff)) 689 780 ENDIF … … 703 794 ENDIF 704 795 ENDDO 705 CALL histwrite_phy(nid_files(iff), 796 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 706 797 s o_t_oce_sic%name,itau_w,zx_tmp_fi2d) 707 798 ENDIF … … 711 802 IF (iflag_coupl>=1) THEN 712 803 IF (o_ale_bl%flag(iff)<=lev_files(iff)) THEN 713 CALL histwrite_phy(nid_files(iff),o_ale_bl%name,itau_w,ale_bl) 804 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 805 $o_ale_bl%name,itau_w,ale_bl) 714 806 ENDIF 715 807 IF (o_alp_bl%flag(iff)<=lev_files(iff)) THEN 716 CALL histwrite_phy(nid_files(iff),o_alp_bl%name,itau_w,alp_bl) 808 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 809 $o_alp_bl%name,itau_w,alp_bl) 717 810 ENDIF 718 811 ENDIF !iflag_coupl>=1 … … 723 816 IF (iflag_wake>=1) THEN 724 817 IF (o_ale_wk%flag(iff)<=lev_files(iff)) THEN 725 CALL histwrite_phy(nid_files(iff),o_ale_wk%name,itau_w,ale_wake) 818 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 819 $o_ale_wk%name,itau_w,ale_wake) 726 820 ENDIF 727 821 IF (o_alp_wk%flag(iff)<=lev_files(iff)) THEN 728 CALL histwrite_phy(nid_files(iff),o_alp_wk%name,itau_w,alp_wake) 822 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 823 $o_alp_wk%name,itau_w,alp_wake) 729 824 ENDIF 730 825 731 826 IF (o_ale%flag(iff)<=lev_files(iff)) THEN 732 CALL histwrite_phy(nid_files(iff),o_ale%name,itau_w,ale) 827 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 828 $o_ale%name,itau_w,ale) 733 829 ENDIF 734 830 IF (o_alp%flag(iff)<=lev_files(iff)) THEN 735 CALL histwrite_phy(nid_files(iff),o_alp%name,itau_w,alp) 831 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 832 $o_alp%name,itau_w,alp) 736 833 ENDIF 737 834 IF (o_cin%flag(iff)<=lev_files(iff)) THEN 738 CALL histwrite_phy(nid_files(iff),o_cin%name,itau_w,cin) 835 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 836 $o_cin%name,itau_w,cin) 739 837 ENDIF 740 838 IF (o_wape%flag(iff)<=lev_files(iff)) THEN 741 CALL histwrite_phy(nid_files(iff),o_WAPE%name,itau_w,wake_pe) 839 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 840 $o_WAPE%name,itau_w,wake_pe) 742 841 ENDIF 743 842 IF (o_wake_h%flag(iff)<=lev_files(iff)) THEN 744 CALL histwrite_phy(nid_files(iff),o_wake_h%name,itau_w,wake_h) 843 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 844 $o_wake_h%name,itau_w,wake_h) 745 845 ENDIF 746 846 747 847 IF (o_wake_s%flag(iff)<=lev_files(iff)) THEN 748 CALL histwrite_phy(nid_files(iff),o_wake_s%name,itau_w,wake_s) 848 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 849 $o_wake_s%name,itau_w,wake_s) 749 850 ENDIF 750 851 751 852 IF (o_wake_deltat%flag(iff)<=lev_files(iff)) THEN 752 CALL histwrite_phy(nid_files(iff),o_wake_deltat%name, 853 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 854 $o_wake_deltat%name, 753 855 $ itau_w,wake_deltat) 754 856 ENDIF 755 857 756 858 IF (o_wake_deltaq%flag(iff)<=lev_files(iff)) THEN 757 CALL histwrite_phy(nid_files(iff),o_wake_deltaq%name, 859 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 860 $o_wake_deltaq%name, 758 861 $ itau_w,wake_deltaq) 759 862 ENDIF 760 863 761 864 IF (o_wake_omg%flag(iff)<=lev_files(iff)) THEN 762 CALL histwrite_phy(nid_files(iff), 865 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 763 866 s o_wake_omg%name,itau_w,wake_omg) 764 867 ENDIF … … 768 871 & /pdtphys 769 872 CALL histwrite_phy(nid_files(iff), 873 $clef_stations(iff), 770 874 & o_dtwak%name,itau_w,zx_tmp_fi3d) 771 875 ENDIF … … 773 877 IF (o_dqwak%flag(iff)<=lev_files(iff)) THEN 774 878 zx_tmp_fi3d(1:klon,1:klev)=d_q_wake(1:klon,1:klev)/pdtphys 775 CALL histwrite_phy(nid_files(iff), 879 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 776 880 & o_dqwak%name,itau_w,zx_tmp_fi3d) 777 881 ENDIF … … 779 883 780 884 IF (o_Vprecip%flag(iff)<=lev_files(iff)) THEN 781 CALL histwrite_phy(nid_files(iff),o_Vprecip%name,itau_w,Vprecip) 885 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 886 $o_Vprecip%name,itau_w,Vprecip) 782 887 ENDIF 783 888 784 889 IF (o_ftd%flag(iff)<=lev_files(iff)) THEN 785 CALL histwrite_phy(nid_files(iff),o_ftd%name,itau_w,ftd) 890 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 891 $o_ftd%name,itau_w,ftd) 786 892 ENDIF 787 893 788 894 IF (o_fqd%flag(iff)<=lev_files(iff)) THEN 789 CALL histwrite_phy(nid_files(iff),o_fqd%name,itau_w,fqd) 895 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 896 $o_fqd%name,itau_w,fqd) 790 897 ENDIF 791 898 ENDIF !(iflag_con.EQ.3) … … 794 901 IF ( o_slab_bils%flag(iff)<=lev_files(iff)) 795 902 $ CALL histwrite_phy( 796 $ nid_files(iff),o_slab_bils%name,itau_w,slab_wfbils) 903 $ nid_files(iff),clef_stations(iff), 904 $o_slab_bils%name,itau_w,slab_wfbils) 797 905 798 906 ENDIF !type_ocean == force/slab 799 907 800 908 IF (o_weakinv%flag(iff)<=lev_files(iff)) THEN 801 CALL histwrite_phy(nid_files(iff), 909 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 802 910 s o_weakinv%name,itau_w,weak_inversion) 803 911 ENDIF 804 912 805 913 IF (o_dthmin%flag(iff)<=lev_files(iff)) THEN 806 CALL histwrite_phy(nid_files(iff),o_dthmin%name,itau_w,dthmin) 914 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 915 $o_dthmin%name,itau_w,dthmin) 807 916 ENDIF 808 917 809 918 IF (o_cldtau%flag(iff)<=lev_files(iff)) THEN 810 CALL histwrite_phy(nid_files(iff),o_cldtau%name,itau_w,cldtau) 919 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 920 $o_cldtau%name,itau_w,cldtau) 811 921 ENDIF 812 922 813 923 IF (o_cldemi%flag(iff)<=lev_files(iff)) THEN 814 CALL histwrite_phy(nid_files(iff),o_cldemi%name,itau_w,cldemi) 924 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 925 $o_cldemi%name,itau_w,cldemi) 815 926 ENDIF 816 927 817 928 IF (o_pr_con_l%flag(iff)<=lev_files(iff)) THEN 818 CALL histwrite_phy(nid_files(iff), 929 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 819 930 s o_pr_con_l%name,itau_w,pmflxr(:,1:klev)) 820 931 ENDIF 821 932 822 933 IF (o_pr_con_i%flag(iff)<=lev_files(iff)) THEN 823 CALL histwrite_phy(nid_files(iff), 934 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 824 935 s o_pr_con_i%name,itau_w,pmflxs(:,1:klev)) 825 936 ENDIF 826 937 827 938 IF (o_pr_lsc_l%flag(iff)<=lev_files(iff)) THEN 828 CALL histwrite_phy(nid_files(iff), 939 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 829 940 s o_pr_lsc_l%name,itau_w,prfl(:,1:klev)) 830 941 ENDIF 831 942 832 943 IF (o_pr_lsc_i%flag(iff)<=lev_files(iff)) THEN 833 CALL histwrite_phy(nid_files(iff), 944 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 834 945 s o_pr_lsc_i%name,itau_w,psfl(:,1:klev)) 835 946 ENDIF 836 947 837 948 IF (o_re%flag(iff)<=lev_files(iff)) THEN 838 CALL histwrite_phy(nid_files(iff),o_re%name,itau_w,re) 949 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 950 $o_re%name,itau_w,re) 839 951 ENDIF 840 952 841 953 IF (o_fl%flag(iff)<=lev_files(iff)) THEN 842 CALL histwrite_phy(nid_files(iff),o_fl%name,itau_w,fl) 954 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 955 $o_fl%name,itau_w,fl) 843 956 ENDIF 844 957 … … 849 962 zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.) 850 963 ENDDO 851 CALL histwrite_phy(nid_files(iff),o_rh2m%name,itau_w,zx_tmp_fi2d) 852 ENDIF 853 964 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 965 $o_rh2m%name,itau_w,zx_tmp_fi2d) 966 ENDIF 967 968 IF (.NOT.clef_stations(iff)) THEN 854 969 IF (o_rh2m_min%flag(iff)<=lev_files(iff)) THEN 855 970 DO i=1, klon 856 971 zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.) 857 972 ENDDO 858 CALL histwrite_phy(nid_files(iff),o_rh2m_min%name, 973 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 974 $o_rh2m_min%name, 859 975 s itau_w,zx_tmp_fi2d) 860 976 ENDIF … … 864 980 zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.) 865 981 ENDDO 866 CALL histwrite_phy(nid_files(iff),o_rh2m_max%name, 982 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 983 $o_rh2m_max%name, 867 984 s itau_w,zx_tmp_fi2d) 868 985 ENDIF 986 ENDIF 869 987 870 988 871 989 IF (o_qsat2m%flag(iff)<=lev_files(iff)) THEN 872 CALL histwrite_phy(nid_files(iff),o_qsat2m%name,itau_w,qsat2m) 990 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 991 $o_qsat2m%name,itau_w,qsat2m) 873 992 ENDIF 874 993 875 994 IF (o_tpot%flag(iff)<=lev_files(iff)) THEN 876 CALL histwrite_phy(nid_files(iff),o_tpot%name,itau_w,tpot) 995 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 996 $o_tpot%name,itau_w,tpot) 877 997 ENDIF 878 998 879 999 IF (o_tpote%flag(iff)<=lev_files(iff)) THEN 880 CALL histwrite_phy(nid_files(iff),o_tpote%name,itau_w,tpote) 1000 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1001 $o_tpote%name,itau_w,tpote) 881 1002 ENDIF 882 1003 883 1004 IF (o_SWnetOR%flag(iff)<=lev_files(iff)) THEN 884 1005 zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, is_ter) 885 CALL histwrite_phy(nid_files(iff), 1006 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 886 1007 s o_SWnetOR%name,itau_w, zx_tmp_fi2d) 887 1008 ENDIF … … 889 1010 IF (o_SWdownOR%flag(iff)<=lev_files(iff)) THEN 890 1011 zx_tmp_fi2d(1:klon) = solsw(1:klon)/(1.-albsol1(1:klon)) 891 CALL histwrite_phy(nid_files(iff), 1012 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 892 1013 s o_SWdownOR%name,itau_w, zx_tmp_fi2d) 893 1014 ENDIF 894 1015 895 1016 IF (o_LWdownOR%flag(iff)<=lev_files(iff)) THEN 896 CALL histwrite_phy(nid_files(iff), 1017 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 897 1018 s o_LWdownOR%name,itau_w,sollwdown) 898 1019 ENDIF 899 1020 900 1021 IF (o_snowl%flag(iff)<=lev_files(iff)) THEN 901 CALL histwrite_phy(nid_files(iff),o_snowl%name,itau_w,snow_lsc) 1022 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1023 $o_snowl%name,itau_w,snow_lsc) 902 1024 ENDIF 903 1025 904 1026 IF (o_solldown%flag(iff)<=lev_files(iff)) THEN 905 CALL histwrite_phy(nid_files(iff), 1027 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 906 1028 s o_solldown%name,itau_w,sollwdown) 907 1029 ENDIF 908 1030 909 1031 IF (o_dtsvdfo%flag(iff)<=lev_files(iff)) THEN 910 CALL histwrite_phy(nid_files(iff), 1032 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 911 1033 s o_dtsvdfo%name,itau_w,d_ts(:,is_oce)) 912 1034 ENDIF 913 1035 914 1036 IF (o_dtsvdft%flag(iff)<=lev_files(iff)) THEN 915 CALL histwrite_phy(nid_files(iff), 1037 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 916 1038 s o_dtsvdft%name,itau_w,d_ts(:,is_ter)) 917 1039 ENDIF 918 1040 919 1041 IF (o_dtsvdfg%flag(iff)<=lev_files(iff)) THEN 920 CALL histwrite_phy(nid_files(iff), 1042 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 921 1043 $ o_dtsvdfg%name,itau_w, d_ts(:,is_lic)) 922 1044 ENDIF 923 1045 924 1046 IF (o_dtsvdfi%flag(iff)<=lev_files(iff)) THEN 925 CALL histwrite_phy(nid_files(iff), 1047 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 926 1048 s o_dtsvdfi%name,itau_w,d_ts(:,is_sic)) 927 1049 ENDIF 928 1050 929 1051 IF (o_rugs%flag(iff)<=lev_files(iff)) THEN 930 CALL histwrite_phy(nid_files(iff),o_rugs%name,itau_w,zxrugs) 1052 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1053 $o_rugs%name,itau_w,zxrugs) 931 1054 ENDIF 932 1055 … … 936 1059 937 1060 IF (o_od550aer%flag(iff)<=lev_files(iff)) THEN 938 CALL histwrite_phy(nid_files(iff),o_od550aer%name,itau_w, 1061 CALL histwrite_phy(nid_files(iff), 1062 $clef_stations(iff), 1063 $o_od550aer%name,itau_w, 939 1064 $ od550aer) 940 1065 ENDIF 941 1066 IF (o_od865aer%flag(iff)<=lev_files(iff)) THEN 942 CALL histwrite_phy(nid_files(iff),o_od865aer%name,itau_w, 1067 CALL histwrite_phy(nid_files(iff), 1068 $clef_stations(iff), 1069 $o_od865aer%name,itau_w, 943 1070 $ od865aer) 944 1071 ENDIF 945 1072 IF (o_absvisaer%flag(iff)<=lev_files(iff)) THEN 946 CALL histwrite_phy(nid_files(iff),o_absvisaer%name,itau_w, 1073 CALL histwrite_phy(nid_files(iff), 1074 $clef_stations(iff), 1075 $o_absvisaer%name,itau_w, 947 1076 $ absvisaer) 948 1077 ENDIF 949 1078 IF (o_od550lt1aer%flag(iff)<=lev_files(iff)) THEN 950 CALL histwrite_phy(nid_files(iff),o_od550lt1aer%name,itau_w, 1079 CALL histwrite_phy(nid_files(iff), 1080 $clef_stations(iff), 1081 $o_od550lt1aer%name,itau_w, 951 1082 $ od550lt1aer) 952 1083 ENDIF 953 1084 954 1085 IF (o_sconcso4%flag(iff)<=lev_files(iff)) THEN 955 CALL histwrite_phy(nid_files(iff),o_sconcso4%name,itau_w, 1086 CALL histwrite_phy(nid_files(iff), 1087 $clef_stations(iff), 1088 $o_sconcso4%name,itau_w, 956 1089 $ sconcso4) 957 1090 ENDIF 958 1091 IF (o_sconcoa%flag(iff)<=lev_files(iff)) THEN 959 CALL histwrite_phy(nid_files(iff),o_sconcoa%name,itau_w, 1092 CALL histwrite_phy(nid_files(iff), 1093 $clef_stations(iff), 1094 $o_sconcoa%name,itau_w, 960 1095 $ sconcoa) 961 1096 ENDIF 962 1097 IF (o_sconcbc%flag(iff)<=lev_files(iff)) THEN 963 CALL histwrite_phy(nid_files(iff),o_sconcbc%name,itau_w, 1098 CALL histwrite_phy(nid_files(iff), 1099 $clef_stations(iff), 1100 $o_sconcbc%name,itau_w, 964 1101 $ sconcbc) 965 1102 ENDIF 966 1103 IF (o_sconcss%flag(iff)<=lev_files(iff)) THEN 967 CALL histwrite_phy(nid_files(iff),o_sconcss%name,itau_w, 1104 CALL histwrite_phy(nid_files(iff), 1105 $clef_stations(iff), 1106 $o_sconcss%name,itau_w, 968 1107 $ sconcss) 969 1108 ENDIF 970 1109 IF (o_sconcdust%flag(iff)<=lev_files(iff)) THEN 971 CALL histwrite_phy(nid_files(iff),o_sconcdust%name,itau_w, 1110 CALL histwrite_phy(nid_files(iff), 1111 $clef_stations(iff), 1112 $o_sconcdust%name,itau_w, 972 1113 $ sconcdust) 973 1114 ENDIF 974 1115 975 1116 IF (o_concso4%flag(iff)<=lev_files(iff)) THEN 976 CALL histwrite_phy(nid_files(iff),o_concso4%name,itau_w, 1117 CALL histwrite_phy(nid_files(iff), 1118 $clef_stations(iff), 1119 $o_concso4%name,itau_w, 977 1120 $ concso4) 978 1121 ENDIF 979 1122 IF (o_concoa%flag(iff)<=lev_files(iff)) THEN 980 CALL histwrite_phy(nid_files(iff),o_concoa%name,itau_w, 1123 CALL histwrite_phy(nid_files(iff), 1124 $clef_stations(iff), 1125 $o_concoa%name,itau_w, 981 1126 $ concoa) 982 1127 ENDIF 983 1128 IF (o_concbc%flag(iff)<=lev_files(iff)) THEN 984 CALL histwrite_phy(nid_files(iff),o_concbc%name,itau_w, 1129 CALL histwrite_phy(nid_files(iff), 1130 $clef_stations(iff), 1131 $o_concbc%name,itau_w, 985 1132 $ concbc) 986 1133 ENDIF 987 1134 IF (o_concss%flag(iff)<=lev_files(iff)) THEN 988 CALL histwrite_phy(nid_files(iff),o_concss%name,itau_w, 1135 CALL histwrite_phy(nid_files(iff), 1136 $clef_stations(iff), 1137 $o_concss%name,itau_w, 989 1138 $ concss) 990 1139 ENDIF 991 1140 IF (o_concdust%flag(iff)<=lev_files(iff)) THEN 992 CALL histwrite_phy(nid_files(iff),o_concdust%name,itau_w, 1141 CALL histwrite_phy(nid_files(iff), 1142 $clef_stations(iff), 1143 $o_concdust%name,itau_w, 993 1144 $ concdust) 994 1145 ENDIF 995 1146 996 1147 IF (o_loadso4%flag(iff)<=lev_files(iff)) THEN 997 CALL histwrite_phy(nid_files(iff),o_loadso4%name,itau_w, 1148 CALL histwrite_phy(nid_files(iff), 1149 $clef_stations(iff), 1150 $o_loadso4%name,itau_w, 998 1151 $ loadso4) 999 1152 ENDIF 1000 1153 IF (o_loadoa%flag(iff)<=lev_files(iff)) THEN 1001 CALL histwrite_phy(nid_files(iff),o_loadoa%name,itau_w, 1154 CALL histwrite_phy(nid_files(iff), 1155 $clef_stations(iff), 1156 $o_loadoa%name,itau_w, 1002 1157 $ loadoa) 1003 1158 ENDIF 1004 1159 IF (o_loadbc%flag(iff)<=lev_files(iff)) THEN 1005 CALL histwrite_phy(nid_files(iff),o_loadbc%name,itau_w, 1160 CALL histwrite_phy(nid_files(iff), 1161 $clef_stations(iff), 1162 $o_loadbc%name,itau_w, 1006 1163 $ loadbc) 1007 1164 ENDIF 1008 1165 IF (o_loadss%flag(iff)<=lev_files(iff)) THEN 1009 CALL histwrite_phy(nid_files(iff),o_loadss%name,itau_w, 1166 CALL histwrite_phy(nid_files(iff), 1167 $clef_stations(iff), 1168 $o_loadss%name,itau_w, 1010 1169 $ loadss) 1011 1170 ENDIF 1012 1171 IF (o_loaddust%flag(iff)<=lev_files(iff)) THEN 1013 CALL histwrite_phy(nid_files(iff),o_loaddust%name,itau_w, 1172 CALL histwrite_phy(nid_files(iff), 1173 $clef_stations(iff), 1174 $o_loaddust%name,itau_w, 1014 1175 $ loaddust) 1015 1176 ENDIF … … 1018 1179 IF (o_tausumaero(naero)%flag(iff)<=lev_files(iff)) THEN 1019 1180 CALL histwrite_phy(nid_files(iff), 1181 $clef_stations(iff), 1020 1182 $ o_tausumaero(naero)%name,itau_w, 1021 1183 $ tausum_aero(:,2,naero) ) … … 1027 1189 IF (ok_ade) THEN 1028 1190 IF (o_topswad%flag(iff)<=lev_files(iff)) THEN 1029 CALL histwrite_phy(nid_files(iff),o_topswad%name,itau_w, 1191 CALL histwrite_phy(nid_files(iff), 1192 $clef_stations(iff), 1193 $o_topswad%name,itau_w, 1030 1194 $ topswad_aero) 1031 1195 ENDIF 1032 1196 IF (o_solswad%flag(iff)<=lev_files(iff)) THEN 1033 CALL histwrite_phy(nid_files(iff),o_solswad%name,itau_w, 1197 CALL histwrite_phy(nid_files(iff), 1198 $clef_stations(iff), 1199 $o_solswad%name,itau_w, 1034 1200 $ solswad_aero) 1035 1201 ENDIF … … 1038 1204 if (new_aod) then 1039 1205 IF (o_swtoaas_nat%flag(iff)<=lev_files(iff)) THEN 1040 CALL histwrite_phy(nid_files(iff),o_swtoaas_nat%name,itau_w, 1206 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1207 $o_swtoaas_nat%name,itau_w, 1041 1208 $ topsw_aero(:,1)) 1042 1209 ENDIF 1043 1210 1044 1211 IF (o_swsrfas_nat%flag(iff)<=lev_files(iff)) THEN 1045 CALL histwrite_phy(nid_files(iff),o_swsrfas_nat%name,itau_w, 1212 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1213 $o_swsrfas_nat%name,itau_w, 1046 1214 $ solsw_aero(:,1)) 1047 1215 ENDIF 1048 1216 1049 1217 IF (o_swtoacs_nat%flag(iff)<=lev_files(iff)) THEN 1050 CALL histwrite_phy(nid_files(iff),o_swtoacs_nat%name,itau_w, 1218 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1219 $o_swtoacs_nat%name,itau_w, 1051 1220 $ topsw0_aero(:,1)) 1052 1221 ENDIF 1053 1222 1054 1223 IF (o_swsrfcs_nat%flag(iff)<=lev_files(iff)) THEN 1055 CALL histwrite_phy(nid_files(iff),o_swsrfcs_nat%name,itau_w, 1224 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1225 $o_swsrfcs_nat%name,itau_w, 1056 1226 $ solsw0_aero(:,1)) 1057 1227 ENDIF … … 1059 1229 !ant 1060 1230 IF (o_swtoaas_ant%flag(iff)<=lev_files(iff)) THEN 1061 CALL histwrite_phy(nid_files(iff),o_swtoaas_ant%name,itau_w, 1231 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1232 $o_swtoaas_ant%name,itau_w, 1062 1233 $ topsw_aero(:,2)) 1063 1234 ENDIF 1064 1235 1065 1236 IF (o_swsrfas_ant%flag(iff)<=lev_files(iff)) THEN 1066 CALL histwrite_phy(nid_files(iff),o_swsrfas_ant%name,itau_w, 1237 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1238 $o_swsrfas_ant%name,itau_w, 1067 1239 $ solsw_aero(:,2)) 1068 1240 ENDIF 1069 1241 1070 1242 IF (o_swtoacs_ant%flag(iff)<=lev_files(iff)) THEN 1071 CALL histwrite_phy(nid_files(iff),o_swtoacs_ant%name,itau_w, 1243 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1244 $o_swtoacs_ant%name,itau_w, 1072 1245 $ topsw0_aero(:,2)) 1073 1246 ENDIF 1074 1247 1075 1248 IF (o_swsrfcs_ant%flag(iff)<=lev_files(iff)) THEN 1076 CALL histwrite_phy(nid_files(iff),o_swsrfcs_ant%name,itau_w, 1249 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1250 $o_swsrfcs_ant%name,itau_w, 1077 1251 $ solsw0_aero(:,2)) 1078 1252 ENDIF … … 1082 1256 if (.not. aerosol_couple) then 1083 1257 IF (o_swtoacf_nat%flag(iff)<=lev_files(iff)) THEN 1084 CALL histwrite_phy(nid_files(iff),o_swtoacf_nat%name,itau_w, 1258 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1259 $o_swtoacf_nat%name,itau_w, 1085 1260 $ topswcf_aero(:,1)) 1086 1261 ENDIF 1087 1262 1088 1263 IF (o_swsrfcf_nat%flag(iff)<=lev_files(iff)) THEN 1089 CALL histwrite_phy(nid_files(iff),o_swsrfcf_nat%name,itau_w, 1264 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1265 $o_swsrfcf_nat%name,itau_w, 1090 1266 $ solswcf_aero(:,1)) 1091 1267 ENDIF 1092 1268 1093 1269 IF (o_swtoacf_ant%flag(iff)<=lev_files(iff)) THEN 1094 CALL histwrite_phy(nid_files(iff),o_swtoacf_ant%name,itau_w, 1270 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1271 $o_swtoacf_ant%name,itau_w, 1095 1272 $ topswcf_aero(:,2)) 1096 1273 ENDIF 1097 1274 1098 1275 IF (o_swsrfcf_ant%flag(iff)<=lev_files(iff)) THEN 1099 CALL histwrite_phy(nid_files(iff),o_swsrfcf_ant%name,itau_w, 1276 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1277 $o_swsrfcf_ant%name,itau_w, 1100 1278 $ solswcf_aero(:,2)) 1101 1279 ENDIF 1102 1280 1103 1281 IF (o_swtoacf_zero%flag(iff)<=lev_files(iff)) THEN 1104 CALL histwrite_phy(nid_files(iff),o_swtoacf_zero%name,itau_w, 1282 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1283 $o_swtoacf_zero%name,itau_w, 1105 1284 $ topswcf_aero(:,3)) 1106 1285 ENDIF 1107 1286 1108 1287 IF (o_swsrfcf_zero%flag(iff)<=lev_files(iff)) THEN 1109 CALL histwrite_phy(nid_files(iff),o_swsrfcf_zero%name,itau_w, 1288 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1289 $o_swsrfcf_zero%name,itau_w, 1110 1290 $ solswcf_aero(:,3)) 1111 1291 ENDIF … … 1119 1299 IF (ok_aie) THEN 1120 1300 IF (o_topswai%flag(iff)<=lev_files(iff)) THEN 1121 CALL histwrite_phy(nid_files(iff),o_topswai%name,itau_w, 1301 CALL histwrite_phy(nid_files(iff), 1302 $clef_stations(iff), 1303 $o_topswai%name,itau_w, 1122 1304 $ topswai_aero) 1123 1305 ENDIF 1124 1306 IF (o_solswai%flag(iff)<=lev_files(iff)) THEN 1125 CALL histwrite_phy(nid_files(iff),o_solswai%name,itau_w, 1307 CALL histwrite_phy(nid_files(iff), 1308 $clef_stations(iff), 1309 $o_solswai%name,itau_w, 1126 1310 $ solswai_aero) 1127 1311 ENDIF 1128 1312 IF (o_scdnc%flag(iff)<=lev_files(iff)) THEN 1129 CALL histwrite_phy(nid_files(iff),o_scdnc%name,itau_w, 1313 CALL histwrite_phy(nid_files(iff), 1314 $clef_stations(iff), 1315 $o_scdnc%name,itau_w, 1130 1316 $ scdnc) 1131 1317 ENDIF 1132 1318 IF (o_cldncl%flag(iff)<=lev_files(iff)) THEN 1133 CALL histwrite_phy(nid_files(iff),o_cldncl%name,itau_w, 1319 CALL histwrite_phy(nid_files(iff), 1320 $clef_stations(iff), 1321 $o_cldncl%name,itau_w, 1134 1322 $ cldncl) 1135 1323 ENDIF 1136 IF (o_reffclws%flag(iff)<=lev_files(iff)) THEN 1137 CALL histwrite_phy(nid_files(iff),o_reffclws%name,itau_w, 1324 IF (o_reffclws%flag(iff)<=lev_files(iff)) THEN 1325 CALL histwrite_phy(nid_files(iff), 1326 $clef_stations(iff), 1327 $o_reffclws%name,itau_w, 1138 1328 $ reffclws) 1139 ENDIF 1140 IF (o_reffclwc%flag(iff)<=lev_files(iff)) THEN 1141 CALL histwrite_phy(nid_files(iff),o_reffclwc%name,itau_w, 1329 ENDIF 1330 IF (o_reffclwc%flag(iff)<=lev_files(iff)) THEN 1331 CALL histwrite_phy(nid_files(iff), 1332 $clef_stations(iff), 1333 $o_reffclwc%name,itau_w, 1142 1334 $ reffclwc) 1143 1335 ENDIF 1144 1336 IF (o_cldnvi%flag(iff)<=lev_files(iff)) THEN 1145 CALL histwrite_phy(nid_files(iff),o_cldnvi%name,itau_w, 1337 CALL histwrite_phy(nid_files(iff), 1338 $clef_stations(iff), 1339 $o_cldnvi%name,itau_w, 1146 1340 $ cldnvi) 1147 1341 ENDIF 1148 1342 IF (o_lcc%flag(iff)<=lev_files(iff)) THEN 1149 CALL histwrite_phy(nid_files(iff),o_lcc%name,itau_w, 1343 CALL histwrite_phy(nid_files(iff), 1344 $clef_stations(iff), 1345 $o_lcc%name,itau_w, 1150 1346 $ lcc) 1151 1347 ENDIF 1152 1348 IF (o_lcc3d%flag(iff)<=lev_files(iff)) THEN 1153 CALL histwrite_phy(nid_files(iff),o_lcc3d%name,itau_w, 1349 CALL histwrite_phy(nid_files(iff), 1350 $clef_stations(iff), 1351 $o_lcc3d%name,itau_w, 1154 1352 $ lcc3d) 1155 1353 ENDIF 1156 1354 IF (o_lcc3dcon%flag(iff)<=lev_files(iff)) THEN 1157 CALL histwrite_phy(nid_files(iff),o_lcc3dcon%name,itau_w, 1355 CALL histwrite_phy(nid_files(iff), 1356 $clef_stations(iff), 1357 $o_lcc3dcon%name,itau_w, 1158 1358 $ lcc3dcon) 1159 1359 ENDIF 1160 1360 IF (o_lcc3dstra%flag(iff)<=lev_files(iff)) THEN 1161 CALL histwrite_phy(nid_files(iff),o_lcc3dstra%name,itau_w, 1361 CALL histwrite_phy(nid_files(iff), 1362 $clef_stations(iff), 1363 $o_lcc3dstra%name,itau_w, 1162 1364 $ lcc3dstra) 1163 1365 ENDIF 1164 1366 IF (o_reffclwtop%flag(iff)<=lev_files(iff)) THEN 1165 CALL histwrite_phy(nid_files(iff),o_reffclwtop%name,itau_w, 1367 CALL histwrite_phy(nid_files(iff), 1368 $clef_stations(iff), 1369 $o_reffclwtop%name,itau_w, 1166 1370 $ reffclwtop) 1167 1371 ENDIF … … 1171 1375 IF (ok_ade .OR. ok_aie) then 1172 1376 IF (o_ec550aer%flag(iff)<=lev_files(iff)) THEN 1173 CALL histwrite_phy(nid_files(iff),o_ec550aer%name,itau_w, 1377 CALL histwrite_phy(nid_files(iff), 1378 $clef_stations(iff), 1379 $o_ec550aer%name,itau_w, 1174 1380 & ec550aer) 1175 1381 ENDIF … … 1177 1383 1178 1384 IF (o_lwcon%flag(iff)<=lev_files(iff)) THEN 1179 CALL histwrite_phy(nid_files(iff),o_lwcon%name,itau_w,flwc) 1385 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1386 $o_lwcon%name,itau_w,flwc) 1180 1387 ENDIF 1181 1388 1182 1389 IF (o_iwcon%flag(iff)<=lev_files(iff)) THEN 1183 CALL histwrite_phy(nid_files(iff),o_iwcon%name,itau_w,fiwc) 1390 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1391 $o_iwcon%name,itau_w,fiwc) 1184 1392 ENDIF 1185 1393 1186 1394 IF (o_temp%flag(iff)<=lev_files(iff)) THEN 1187 CALL histwrite_phy(nid_files(iff),o_temp%name,itau_w,t_seri) 1395 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1396 $o_temp%name,itau_w,t_seri) 1188 1397 ENDIF 1189 1398 1190 1399 IF (o_theta%flag(iff)<=lev_files(iff)) THEN 1191 CALL histwrite_phy(nid_files(iff),o_theta%name,itau_w,theta) 1400 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1401 $o_theta%name,itau_w,theta) 1192 1402 ENDIF 1193 1403 1194 1404 IF (o_ovapinit%flag(iff)<=lev_files(iff)) THEN 1195 CALL histwrite_phy(nid_files(iff),o_ovapinit%name,itau_w, 1405 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1406 $o_ovapinit%name,itau_w, 1196 1407 $ qx(:,:,ivap)) 1197 1408 ENDIF 1198 1409 1199 1410 IF (o_ovap%flag(iff)<=lev_files(iff)) THEN 1200 CALL histwrite_phy(nid_files(iff), 1411 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1201 1412 $ o_ovap%name,itau_w,q_seri) 1202 1413 ENDIF 1203 1414 1204 1415 IF (o_geop%flag(iff)<=lev_files(iff)) THEN 1205 CALL histwrite_phy(nid_files(iff),o_geop%name,itau_w,zphi) 1416 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1417 $o_geop%name,itau_w,zphi) 1206 1418 ENDIF 1207 1419 1208 1420 IF (o_vitu%flag(iff)<=lev_files(iff)) THEN 1209 CALL histwrite_phy(nid_files(iff),o_vitu%name,itau_w,u_seri) 1421 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1422 $o_vitu%name,itau_w,u_seri) 1210 1423 ENDIF 1211 1424 1212 1425 IF (o_vitv%flag(iff)<=lev_files(iff)) THEN 1213 CALL histwrite_phy(nid_files(iff),o_vitv%name,itau_w,v_seri) 1426 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1427 $o_vitv%name,itau_w,v_seri) 1214 1428 ENDIF 1215 1429 1216 1430 IF (o_vitw%flag(iff)<=lev_files(iff)) THEN 1217 CALL histwrite_phy(nid_files(iff),o_vitw%name,itau_w,omega) 1431 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1432 $o_vitw%name,itau_w,omega) 1218 1433 ENDIF 1219 1434 1220 1435 IF (o_pres%flag(iff)<=lev_files(iff)) THEN 1221 CALL histwrite_phy(nid_files(iff),o_pres%name,itau_w,pplay) 1436 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1437 $o_pres%name,itau_w,pplay) 1222 1438 ENDIF 1223 1439 1224 1440 IF (o_paprs%flag(iff)<=lev_files(iff)) THEN 1225 CALL histwrite_phy(nid_files(iff),o_paprs%name, 1441 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1442 $o_paprs%name, 1226 1443 s itau_w,paprs(:,1:klev)) 1227 1444 ENDIF 1228 1445 1446 IF (o_zfull%flag(iff)<=lev_files(iff)) THEN 1447 DO i=1, klon 1448 zx_tmp_fi3d1(i,1)= pphis(i)/RG 1449 !020611 zx_tmp_fi3d(i,1)= pphis(i)/RG 1450 ENDDO 1451 DO k=1, klev 1452 !020611 DO k=1, klev-1 1453 DO i=1, klon 1454 !020611 zx_tmp_fi3d(i,k+1)= zx_tmp_fi3d(i,k) - (t_seri(i,k) *RD * 1455 zx_tmp_fi3d1(i,k+1)= zx_tmp_fi3d1(i,k) - (t_seri(i,k) *RD * 1456 $ (paprs(i,k+1) - paprs(i,k))) / ( pplay(i,k) * RG ) 1457 ENDDO 1458 ENDDO 1459 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1460 $o_zfull%name,itau_w,zx_tmp_fi3d1(:,2:klevp1)) 1461 !020611 $o_zfull%name,itau_w,zx_tmp_fi3d) 1462 ENDIF 1463 1464 IF (o_zhalf%flag(iff)<=lev_files(iff)) THEN 1465 DO i=1, klon 1466 zx_tmp_fi3d(i,1)= pphis(i)/RG - ( 1467 $ (t_seri(i,1)+zxtsol(i))/2. *RD * 1468 $ (pplay(i,1) - paprs(i,1)))/( (paprs(i,1)+pplay(i,1))/2. * RG) 1469 ENDDO 1470 DO k=1, klev-1 1471 DO i=1, klon 1472 zx_tmp_fi3d(i,k+1)= zx_tmp_fi3d(i,k) - ( 1473 $ (t_seri(i,k)+t_seri(i,k+1))/2. *RD * 1474 $ (pplay(i,k+1) - pplay(i,k))) / ( paprs(i,k) * RG ) 1475 ENDDO 1476 ENDDO 1477 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1478 $o_zhalf%name,itau_w,zx_tmp_fi3d) 1479 ENDIF 1480 1229 1481 IF (o_rneb%flag(iff)<=lev_files(iff)) THEN 1230 CALL histwrite_phy(nid_files(iff),o_rneb%name,itau_w,cldfra) 1482 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1483 $o_rneb%name,itau_w,cldfra) 1231 1484 ENDIF 1232 1485 1233 1486 IF (o_rnebcon%flag(iff)<=lev_files(iff)) THEN 1234 CALL histwrite_phy(nid_files(iff),o_rnebcon%name,itau_w,rnebcon) 1487 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1488 $o_rnebcon%name,itau_w,rnebcon) 1235 1489 ENDIF 1236 1490 1237 1491 IF (o_rhum%flag(iff)<=lev_files(iff)) THEN 1238 CALL histwrite_phy(nid_files(iff),o_rhum%name,itau_w,zx_rh) 1492 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1493 $o_rhum%name,itau_w,zx_rh) 1239 1494 ENDIF 1240 1495 1241 1496 IF (o_ozone%flag(iff)<=lev_files(iff)) THEN 1242 CALL histwrite_phy(nid_files(iff), o_ozone%name, itau_w, 1497 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1498 $o_ozone%name, itau_w, 1243 1499 $ wo(:, :, 1) * dobson_u * 1e3 / zmasse / rmo3 * rmd) 1244 1500 ENDIF … … 1246 1502 IF (o_ozone_light%flag(iff)<=lev_files(iff) .and. 1247 1503 $ read_climoz == 2) THEN 1248 CALL histwrite_phy(nid_files(iff), o_ozone_light%name, itau_w, 1504 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1505 $o_ozone_light%name, itau_w, 1249 1506 $ wo(:, :, 2) * dobson_u * 1e3 / zmasse / rmo3 * rmd) 1250 1507 ENDIF 1251 1508 1252 1509 IF (o_dtphy%flag(iff)<=lev_files(iff)) THEN 1253 CALL histwrite_phy(nid_files(iff),o_dtphy%name,itau_w,d_t) 1510 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1511 $o_dtphy%name,itau_w,d_t) 1254 1512 ENDIF 1255 1513 1256 1514 IF (o_dqphy%flag(iff)<=lev_files(iff)) THEN 1257 CALL histwrite_phy(nid_files(iff), 1515 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1258 1516 s o_dqphy%name,itau_w, d_qx(:,:,ivap)) 1259 1517 ENDIF … … 1262 1520 IF (o_albe_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN 1263 1521 zx_tmp_fi2d(1 : klon) = falb1( 1 : klon, nsrf) 1264 CALL histwrite_phy(nid_files(iff), 1522 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1265 1523 s o_albe_srf(nsrf)%name,itau_w, 1266 1524 $ zx_tmp_fi2d) … … 1269 1527 IF (o_rugs_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN 1270 1528 zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf) 1271 CALL histwrite_phy(nid_files(iff), 1529 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1272 1530 s o_rugs_srf(nsrf)%name,itau_w, 1273 1531 $ zx_tmp_fi2d) … … 1276 1534 IF (o_ages_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN 1277 1535 zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf) 1278 CALL histwrite_phy(nid_files(iff), 1536 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1279 1537 s o_ages_srf(nsrf)%name,itau_w 1280 1538 $ ,zx_tmp_fi2d) … … 1283 1541 1284 1542 IF (o_alb1%flag(iff)<=lev_files(iff)) THEN 1285 CALL histwrite_phy(nid_files(iff),o_alb1%name,itau_w,albsol1) 1543 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1544 $o_alb1%name,itau_w,albsol1) 1286 1545 ENDIF 1287 1546 1288 1547 IF (o_alb2%flag(iff)<=lev_files(iff)) THEN 1289 CALL histwrite_phy(nid_files(iff),o_alb2%name,itau_w,albsol2) 1548 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1549 $o_alb2%name,itau_w,albsol2) 1290 1550 ENDIF 1291 1551 … … 1300 1560 enddo 1301 1561 IF (o_tke%flag(iff)<=lev_files(iff)) THEN 1302 CALL histwrite_phy(nid_files(iff),o_tke%name,itau_w,zx_tmp_fi3d) 1303 ENDIF 1304 1562 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1563 $o_tke%name,itau_w,zx_tmp_fi3d) 1564 ENDIF 1565 1566 IF (.NOT.clef_stations(iff)) THEN 1305 1567 IF (o_tke_max%flag(iff)<=lev_files(iff)) THEN 1306 CALL histwrite_phy(nid_files(iff), 1568 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1307 1569 s o_tke_max%name,itau_w,zx_tmp_fi3d) 1308 1570 ENDIF 1571 ENDIF 1309 1572 endif 1310 1573 1311 1574 IF (o_kz%flag(iff)<=lev_files(iff)) THEN 1312 CALL histwrite_phy(nid_files(iff),o_kz%name,itau_w,coefh) 1313 ENDIF 1314 1575 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1576 $o_kz%name,itau_w,coefh) 1577 ENDIF 1578 1579 IF (.NOT.clef_stations(iff)) THEN 1315 1580 IF (o_kz_max%flag(iff)<=lev_files(iff)) THEN 1316 CALL histwrite_phy(nid_files(iff),o_kz_max%name,itau_w,coefh) 1317 ENDIF 1581 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1582 $o_kz_max%name,itau_w,coefh) 1583 ENDIF 1584 ENDIF 1318 1585 1319 1586 IF (o_clwcon%flag(iff)<=lev_files(iff)) THEN 1320 CALL histwrite_phy(nid_files(iff),o_clwcon%name,itau_w,clwcon0) 1587 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1588 $o_clwcon%name,itau_w,clwcon0) 1321 1589 ENDIF 1322 1590 1323 1591 IF (o_dtdyn%flag(iff)<=lev_files(iff)) THEN 1324 CALL histwrite_phy(nid_files(iff),o_dtdyn%name,itau_w,d_t_dyn) 1592 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1593 $o_dtdyn%name,itau_w,d_t_dyn) 1325 1594 ENDIF 1326 1595 1327 1596 IF (o_dqdyn%flag(iff)<=lev_files(iff)) THEN 1328 CALL histwrite_phy(nid_files(iff),o_dqdyn%name,itau_w,d_q_dyn) 1597 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1598 $o_dqdyn%name,itau_w,d_q_dyn) 1329 1599 ENDIF 1330 1600 1331 1601 IF (o_dudyn%flag(iff)<=lev_files(iff)) THEN 1332 CALL histwrite_phy(nid_files(iff),o_dudyn%name,itau_w,d_u_dyn) 1602 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1603 $o_dudyn%name,itau_w,d_u_dyn) 1333 1604 ENDIF 1334 1605 1335 1606 IF (o_dvdyn%flag(iff)<=lev_files(iff)) THEN 1336 CALL histwrite_phy(nid_files(iff),o_dvdyn%name,itau_w,d_v_dyn) 1607 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1608 $o_dvdyn%name,itau_w,d_v_dyn) 1337 1609 ENDIF 1338 1610 1339 1611 IF (o_dtcon%flag(iff)<=lev_files(iff)) THEN 1340 1612 zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys 1341 CALL histwrite_phy(nid_files(iff),o_dtcon%name,itau_w,zx_tmp_fi3d) 1342 ENDIF 1613 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1614 $o_dtcon%name,itau_w,zx_tmp_fi3d) 1615 ENDIF 1616 1617 if(iflag_thermals.eq.1)then 1618 IF (o_tntc%flag(iff)<=lev_files(iff)) THEN 1619 zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys + 1620 $ d_t_ajsb(1:klon,1:klev)/pdtphys 1621 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1622 $o_tntc%name,itau_w,zx_tmp_fi3d) 1623 ENDIF 1624 else if(iflag_thermals.gt.1.and.iflag_wake.EQ.1)then 1625 IF (o_tntc%flag(iff)<=lev_files(iff)) THEN 1626 zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys + 1627 $ d_t_ajs(1:klon,1:klev)/pdtphys + 1628 $ d_t_wake(1:klon,1:klev)/pdtphys 1629 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1630 $o_tntc%name,itau_w,zx_tmp_fi3d) 1631 ENDIF 1632 endif 1343 1633 1344 1634 IF (o_ducon%flag(iff)<=lev_files(iff)) THEN 1345 1635 zx_tmp_fi3d(1:klon,1:klev)=d_u_con(1:klon,1:klev)/pdtphys 1346 CALL histwrite_phy(nid_files(iff),o_ducon%name,itau_w,zx_tmp_fi3d) 1636 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1637 $o_ducon%name,itau_w,zx_tmp_fi3d) 1347 1638 ENDIF 1348 1639 1349 1640 IF (o_dqcon%flag(iff)<=lev_files(iff)) THEN 1350 1641 zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys 1351 CALL histwrite_phy(nid_files(iff),o_dqcon%name,itau_w,zx_tmp_fi3d) 1352 ENDIF 1642 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1643 $o_dqcon%name,itau_w,zx_tmp_fi3d) 1644 ENDIF 1645 1646 if(iflag_thermals.eq.1)then 1647 IF (o_tnhusc%flag(iff)<=lev_files(iff)) THEN 1648 zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys 1649 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1650 $o_tnhusc%name,itau_w,zx_tmp_fi3d) 1651 ENDIF 1652 else if(iflag_thermals.gt.1.and.iflag_wake.EQ.1)then 1653 IF (o_tnhusc%flag(iff)<=lev_files(iff)) THEN 1654 zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys + 1655 $ d_q_ajs(1:klon,1:klev)/pdtphys + 1656 $ d_q_wake(1:klon,1:klev)/pdtphys 1657 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1658 $o_tnhusc%name,itau_w,zx_tmp_fi3d) 1659 ENDIF 1660 endif 1353 1661 1354 1662 IF (o_dtlsc%flag(iff)<=lev_files(iff)) THEN 1355 1663 zx_tmp_fi3d(1:klon,1:klev)=d_t_lsc(1:klon,1:klev)/pdtphys 1356 CALL histwrite_phy(nid_files(iff),o_dtlsc%name,itau_w,zx_tmp_fi3d) 1664 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1665 $o_dtlsc%name,itau_w,zx_tmp_fi3d) 1357 1666 ENDIF 1358 1667 … … 1360 1669 zx_tmp_fi3d(1:klon, 1:klev)=(d_t_lsc(1:klon,1:klev)+ 1361 1670 $ d_t_eva(1:klon,1:klev))/pdtphys 1362 CALL histwrite_phy(nid_files(iff), 1671 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1363 1672 s o_dtlschr%name,itau_w,zx_tmp_fi3d) 1364 1673 ENDIF … … 1366 1675 IF (o_dqlsc%flag(iff)<=lev_files(iff)) THEN 1367 1676 zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys 1368 CALL histwrite_phy(nid_files(iff),o_dqlsc%name,itau_w,zx_tmp_fi3d) 1677 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1678 $o_dqlsc%name,itau_w,zx_tmp_fi3d) 1369 1679 ENDIF 1370 1680 … … 1376 1686 zx_tmp_fi3d(1:klon,1:klev)=d_t_lscth(1:klon,1:klev)/pdtphys 1377 1687 CALL 1378 s histwrite_phy(nid_files(iff),o_dtlscth%name,itau_w,zx_tmp_fi3d) 1688 s histwrite_phy(nid_files(iff),clef_stations(iff), 1689 $o_dtlscth%name,itau_w,zx_tmp_fi3d) 1379 1690 ENDIF 1380 1691 … … 1382 1693 zx_tmp_fi3d(1:klon,1:klev)=d_t_lscst(1:klon,1:klev)/pdtphys 1383 1694 CALL 1384 s histwrite_phy(nid_files(iff),o_dtlscst%name,itau_w,zx_tmp_fi3d) 1695 s histwrite_phy(nid_files(iff),clef_stations(iff), 1696 $o_dtlscst%name,itau_w,zx_tmp_fi3d) 1385 1697 ENDIF 1386 1698 … … 1388 1700 zx_tmp_fi3d(1:klon,1:klev)=d_q_lscth(1:klon,1:klev)/pdtphys 1389 1701 CALL 1390 s histwrite_phy(nid_files(iff),o_dqlscth%name,itau_w,zx_tmp_fi3d) 1702 s histwrite_phy(nid_files(iff),clef_stations(iff), 1703 $o_dqlscth%name,itau_w,zx_tmp_fi3d) 1391 1704 ENDIF 1392 1705 … … 1394 1707 zx_tmp_fi3d(1:klon,1:klev)=d_q_lscst(1:klon,1:klev)/pdtphys 1395 1708 CALL 1396 s histwrite_phy(nid_files(iff),o_dqlscst%name,itau_w,zx_tmp_fi3d) 1709 s histwrite_phy(nid_files(iff),clef_stations(iff), 1710 $o_dqlscst%name,itau_w,zx_tmp_fi3d) 1397 1711 ENDIF 1398 1712 1399 1713 IF (o_plulth%flag(iff)<=lev_files(iff)) THEN 1400 CALL histwrite_phy(nid_files(iff),o_plulth%name,itau_w,plul_th) 1714 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1715 $o_plulth%name,itau_w,plul_th) 1401 1716 ENDIF 1402 1717 1403 1718 IF (o_plulst%flag(iff)<=lev_files(iff)) THEN 1404 CALL histwrite_phy(nid_files(iff),o_plulst%name,itau_w,plul_st) 1719 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1720 $o_plulst%name,itau_w,plul_st) 1405 1721 ENDIF 1406 1722 … … 1416 1732 IF (o_ptconvth%flag(iff)<=lev_files(iff)) THEN 1417 1733 CALL 1418 s histwrite_phy(nid_files(iff),o_ptconvth%name,itau_w,zx_tmp_fi3d) 1734 s histwrite_phy(nid_files(iff),clef_stations(iff), 1735 $o_ptconvth%name,itau_w,zx_tmp_fi3d) 1419 1736 ENDIF 1420 1737 … … 1424 1741 IF (o_ptconvth%flag(iff)<=lev_files(iff)) THEN 1425 1742 CALL 1426 s histwrite_phy(nid_files(iff),o_lmaxth%name,itau_w,zx_tmp_fi2d) 1743 s histwrite_phy(nid_files(iff),clef_stations(iff), 1744 $o_lmaxth%name,itau_w,zx_tmp_fi2d) 1427 1745 ENDIF 1428 1746 … … 1433 1751 IF (o_dtvdf%flag(iff)<=lev_files(iff)) THEN 1434 1752 zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)/pdtphys 1435 CALL histwrite_phy(nid_files(iff),o_dtvdf%name,itau_w,zx_tmp_fi3d) 1753 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1754 $o_dtvdf%name,itau_w,zx_tmp_fi3d) 1436 1755 ENDIF 1437 1756 1438 1757 IF (o_dqvdf%flag(iff)<=lev_files(iff)) THEN 1439 1758 zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf(1:klon,1:klev)/pdtphys 1440 CALL histwrite_phy(nid_files(iff),o_dqvdf%name,itau_w,zx_tmp_fi3d) 1759 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1760 $o_dqvdf%name,itau_w,zx_tmp_fi3d) 1441 1761 ENDIF 1442 1762 1443 1763 IF (o_dteva%flag(iff)<=lev_files(iff)) THEN 1444 1764 zx_tmp_fi3d(1:klon,1:klev)=d_t_eva(1:klon,1:klev)/pdtphys 1445 CALL histwrite_phy(nid_files(iff),o_dteva%name,itau_w,zx_tmp_fi3d) 1765 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1766 $o_dteva%name,itau_w,zx_tmp_fi3d) 1446 1767 ENDIF 1447 1768 1448 1769 IF (o_dqeva%flag(iff)<=lev_files(iff)) THEN 1449 1770 zx_tmp_fi3d(1:klon,1:klev)=d_q_eva(1:klon,1:klev)/pdtphys 1450 CALL histwrite_phy(nid_files(iff),o_dqeva%name,itau_w,zx_tmp_fi3d) 1771 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1772 $o_dqeva%name,itau_w,zx_tmp_fi3d) 1451 1773 ENDIF 1452 1774 … … 1454 1776 zpt_conv = 0. 1455 1777 where (ptconv) zpt_conv = 1. 1456 CALL histwrite_phy(nid_files(iff),o_ptconv%name,itau_w,zpt_conv) 1778 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1779 $o_ptconv%name,itau_w,zpt_conv) 1457 1780 ENDIF 1458 1781 1459 1782 IF (o_ratqs%flag(iff)<=lev_files(iff)) THEN 1460 CALL histwrite_phy(nid_files(iff),o_ratqs%name,itau_w,ratqs) 1783 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1784 $o_ratqs%name,itau_w,ratqs) 1461 1785 ENDIF 1462 1786 1463 1787 IF (o_dtthe%flag(iff)<=lev_files(iff)) THEN 1464 zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev)/pdtphys 1465 CALL histwrite_phy(nid_files(iff),o_dtthe%name,itau_w,zx_tmp_fi3d) 1788 zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev)/pdtphys - 1789 $ d_t_ajsb(1:klon,1:klev)/pdtphys 1790 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1791 $o_dtthe%name,itau_w,zx_tmp_fi3d) 1466 1792 ENDIF 1467 1793 … … 1470 1796 ! Pour l instant 0 a y reflichir pour les thermiques 1471 1797 zx_tmp_fi2d=0. 1472 CALL histwrite_phy(nid_files(iff),o_ftime_th%name, 1798 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1799 $o_ftime_th%name, 1473 1800 s itau_w,zx_tmp_fi2d) 1474 1801 ENDIF 1475 1802 1476 1803 IF (o_f_th%flag(iff)<=lev_files(iff)) THEN 1477 CALL histwrite_phy(nid_files(iff),o_f_th%name,itau_w,fm_therm) 1804 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1805 $o_f_th%name,itau_w,fm_therm) 1478 1806 ENDIF 1479 1807 1480 1808 IF (o_e_th%flag(iff)<=lev_files(iff)) THEN 1481 CALL histwrite_phy(nid_files(iff),o_e_th%name,itau_w,entr_therm) 1809 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1810 $o_e_th%name,itau_w,entr_therm) 1482 1811 ENDIF 1483 1812 1484 1813 IF (o_w_th%flag(iff)<=lev_files(iff)) THEN 1485 CALL histwrite_phy(nid_files(iff),o_w_th%name,itau_w,zw2) 1814 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1815 $o_w_th%name,itau_w,zw2) 1486 1816 ENDIF 1487 1817 1488 1818 IF (o_q_th%flag(iff)<=lev_files(iff)) THEN 1489 CALL histwrite_phy(nid_files(iff),o_q_th%name,itau_w,zqasc) 1819 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1820 $o_q_th%name,itau_w,zqasc) 1490 1821 ENDIF 1491 1822 1492 1823 IF (o_lambda_th%flag(iff)<=lev_files(iff)) THEN 1493 CALL histwrite_phy(nid_files(iff), 1824 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1494 1825 s o_lambda_th%name,itau_w,lambda_th) 1495 1826 ENDIF 1496 1827 1497 1828 IF (o_a_th%flag(iff)<=lev_files(iff)) THEN 1498 CALL histwrite_phy(nid_files(iff),o_a_th%name,itau_w,fraca) 1829 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1830 $o_a_th%name,itau_w,fraca) 1499 1831 ENDIF 1500 1832 1501 1833 IF (o_d_th%flag(iff)<=lev_files(iff)) THEN 1502 CALL histwrite_phy(nid_files(iff), o_d_th%name,itau_w,detr_therm)1503 ENDIF1504 1834 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1835 $o_d_th%name,itau_w,detr_therm) 1836 ENDIF 1505 1837 1506 1838 IF (o_f0_th%flag(iff)<=lev_files(iff)) THEN 1507 CALL histwrite_phy(nid_files(iff),o_f0_th%name,itau_w,f0) 1839 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1840 $o_f0_th%name,itau_w,f0) 1508 1841 ENDIF 1509 1842 1510 1843 IF (o_zmax_th%flag(iff)<=lev_files(iff)) THEN 1511 CALL histwrite_phy(nid_files(iff),o_zmax_th%name,itau_w,zmax_th) 1844 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1845 $o_zmax_th%name,itau_w,zmax_th) 1512 1846 ENDIF 1513 1847 1514 1848 IF (o_dqthe%flag(iff)<=lev_files(iff)) THEN 1515 zx_tmp_fi3d(1:klon,1:klev)=d_q_ajs(1:klon,1:klev)/pdtphys 1516 CALL histwrite_phy(nid_files(iff),o_dqthe%name,itau_w,zx_tmp_fi3d) 1849 zx_tmp_fi3d(1:klon,1:klev)=d_q_ajs(1:klon,1:klev)/pdtphys - 1850 $ d_q_ajsb(1:klon,1:klev)/pdtphys 1851 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1852 $o_dqthe%name,itau_w,zx_tmp_fi3d) 1517 1853 ENDIF 1518 1854 … … 1521 1857 IF (o_dtajs%flag(iff)<=lev_files(iff)) THEN 1522 1858 zx_tmp_fi3d(1:klon,1:klev)=d_t_ajsb(1:klon,1:klev)/pdtphys 1523 CALL histwrite_phy(nid_files(iff),o_dtajs%name,itau_w,zx_tmp_fi3d) 1859 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1860 $o_dtajs%name,itau_w,zx_tmp_fi3d) 1524 1861 ENDIF 1525 1862 1526 1863 IF (o_dqajs%flag(iff)<=lev_files(iff)) THEN 1527 1864 zx_tmp_fi3d(1:klon,1:klev)=d_q_ajsb(1:klon,1:klev)/pdtphys 1528 CALL histwrite_phy(nid_files(iff),o_dqajs%name,itau_w,zx_tmp_fi3d) 1865 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1866 $o_dqajs%name,itau_w,zx_tmp_fi3d) 1529 1867 ENDIF 1530 1868 1531 1869 IF (o_dtswr%flag(iff)<=lev_files(iff)) THEN 1532 1870 zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY 1533 CALL histwrite_phy(nid_files(iff),o_dtswr%name,itau_w,zx_tmp_fi3d) 1871 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1872 $o_dtswr%name,itau_w,zx_tmp_fi3d) 1534 1873 ENDIF 1535 1874 1536 1875 IF (o_dtsw0%flag(iff)<=lev_files(iff)) THEN 1537 1876 zx_tmp_fi3d(1:klon,1:klev)=heat0(1:klon,1:klev)/RDAY 1538 CALL histwrite_phy(nid_files(iff),o_dtsw0%name,itau_w,zx_tmp_fi3d) 1877 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1878 $o_dtsw0%name,itau_w,zx_tmp_fi3d) 1539 1879 ENDIF 1540 1880 1541 1881 IF (o_dtlwr%flag(iff)<=lev_files(iff)) THEN 1542 1882 zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)/RDAY 1543 CALL histwrite_phy(nid_files(iff),o_dtlwr%name,itau_w,zx_tmp_fi3d) 1883 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1884 $o_dtlwr%name,itau_w,zx_tmp_fi3d) 1544 1885 ENDIF 1545 1886 1546 1887 IF (o_dtlw0%flag(iff)<=lev_files(iff)) THEN 1547 1888 zx_tmp_fi3d(1:klon,1:klev)=-1.*cool0(1:klon,1:klev)/RDAY 1548 CALL histwrite_phy(nid_files(iff),o_dtlw0%name,itau_w,zx_tmp_fi3d) 1889 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1890 $o_dtlw0%name,itau_w,zx_tmp_fi3d) 1549 1891 ENDIF 1550 1892 1551 1893 IF (o_dtec%flag(iff)<=lev_files(iff)) THEN 1552 1894 zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev) 1553 CALL histwrite_phy(nid_files(iff),o_dtec%name,itau_w,zx_tmp_fi3d) 1895 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1896 $o_dtec%name,itau_w,zx_tmp_fi3d) 1554 1897 ENDIF 1555 1898 1556 1899 IF (o_duvdf%flag(iff)<=lev_files(iff)) THEN 1557 1900 zx_tmp_fi3d(1:klon,1:klev)=d_u_vdf(1:klon,1:klev)/pdtphys 1558 CALL histwrite_phy(nid_files(iff),o_duvdf%name,itau_w,zx_tmp_fi3d) 1901 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1902 $o_duvdf%name,itau_w,zx_tmp_fi3d) 1559 1903 ENDIF 1560 1904 1561 1905 IF (o_dvvdf%flag(iff)<=lev_files(iff)) THEN 1562 1906 zx_tmp_fi3d(1:klon,1:klev)=d_v_vdf(1:klon,1:klev)/pdtphys 1563 CALL histwrite_phy(nid_files(iff),o_dvvdf%name,itau_w,zx_tmp_fi3d) 1907 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1908 $o_dvvdf%name,itau_w,zx_tmp_fi3d) 1564 1909 ENDIF 1565 1910 … … 1567 1912 IF (o_duoro%flag(iff)<=lev_files(iff)) THEN 1568 1913 zx_tmp_fi3d(1:klon,1:klev)=d_u_oro(1:klon,1:klev)/pdtphys 1569 CALL histwrite_phy(nid_files(iff),o_duoro%name,itau_w,zx_tmp_fi3d) 1914 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1915 $o_duoro%name,itau_w,zx_tmp_fi3d) 1570 1916 ENDIF 1571 1917 1572 1918 IF (o_dvoro%flag(iff)<=lev_files(iff)) THEN 1573 1919 zx_tmp_fi3d(1:klon,1:klev)=d_v_oro(1:klon,1:klev)/pdtphys 1574 CALL histwrite_phy(nid_files(iff),o_dvoro%name,itau_w,zx_tmp_fi3d) 1920 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1921 $o_dvoro%name,itau_w,zx_tmp_fi3d) 1922 ENDIF 1923 1924 IF (o_dtoro%flag(iff)<=lev_files(iff)) THEN 1925 zx_tmp_fi3d(1:klon,1:klev)=d_t_oro(1:klon,1:klev)/pdtphys 1926 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1927 $o_dtoro%name,itau_w,zx_tmp_fi3d) 1575 1928 ENDIF 1576 1929 ENDIF … … 1579 1932 IF (o_dulif%flag(iff)<=lev_files(iff)) THEN 1580 1933 zx_tmp_fi3d(1:klon,1:klev)=d_u_lif(1:klon,1:klev)/pdtphys 1581 CALL histwrite_phy(nid_files(iff),o_dulif%name,itau_w,zx_tmp_fi3d) 1934 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1935 $o_dulif%name,itau_w,zx_tmp_fi3d) 1582 1936 ENDIF 1583 1937 1584 1938 IF (o_dvlif%flag(iff)<=lev_files(iff)) THEN 1585 1939 zx_tmp_fi3d(1:klon,1:klev)=d_v_lif(1:klon,1:klev)/pdtphys 1586 CALL histwrite_phy(nid_files(iff),o_dvlif%name,itau_w,zx_tmp_fi3d) 1587 ENDIF 1588 ENDIF 1940 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1941 $o_dvlif%name,itau_w,zx_tmp_fi3d) 1942 ENDIF 1943 1944 IF (o_dtlif%flag(iff)<=lev_files(iff)) THEN 1945 zx_tmp_fi3d(1:klon,1:klev)=d_t_lif(1:klon,1:klev)/pdtphys 1946 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1947 $o_dtlif%name,itau_w,zx_tmp_fi3d) 1948 ENDIF 1949 ENDIF 1950 1951 IF (ok_hines) THEN 1952 IF (o_duhin%flag(iff)<=lev_files(iff)) THEN 1953 zx_tmp_fi3d(1:klon,1:klev)=d_u_hin(1:klon,1:klev)/pdtphys 1954 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1955 $o_duhin%name,itau_w,zx_tmp_fi3d) 1956 ENDIF 1957 1958 1959 IF (o_dvhin%flag(iff)<=lev_files(iff)) THEN 1960 zx_tmp_fi3d(1:klon,1:klev)=d_v_hin(1:klon,1:klev)/pdtphys 1961 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1962 $o_dvhin%name,itau_w,zx_tmp_fi3d) 1963 ENDIF 1964 1965 IF (o_dthin%flag(iff)<=lev_files(iff)) THEN 1966 zx_tmp_fi3d(1:klon,1:klev)=d_t_hin(1:klon,1:klev)/pdtphys 1967 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1968 $o_dthin%name,itau_w,zx_tmp_fi3d) 1969 ENDIF 1970 ENDIF 1971 1972 IF (o_rsu%flag(iff)<=lev_files(iff)) THEN 1973 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1974 $o_rsu%name,itau_w,swup) 1975 ENDIF 1976 IF (o_rsd%flag(iff)<=lev_files(iff)) THEN 1977 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1978 $o_rsd%name,itau_w,swdn) 1979 ENDIF 1980 IF (o_rlu%flag(iff)<=lev_files(iff)) THEN 1981 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1982 $o_rlu%name,itau_w,lwup) 1983 ENDIF 1984 IF (o_rld%flag(iff)<=lev_files(iff)) THEN 1985 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1986 $o_rld%name,itau_w,lwdn) 1987 ENDIF 1988 1989 IF (o_rsucs%flag(iff)<=lev_files(iff)) THEN 1990 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1991 $o_rsucs%name,itau_w,swup0) 1992 ENDIF 1993 IF (o_rsdcs%flag(iff)<=lev_files(iff)) THEN 1994 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1995 $o_rsdcs%name,itau_w,swdn0) 1996 ENDIF 1997 IF (o_rlucs%flag(iff)<=lev_files(iff)) THEN 1998 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1999 $o_rlucs%name,itau_w,lwup0) 2000 ENDIF 2001 IF (o_rldcs%flag(iff)<=lev_files(iff)) THEN 2002 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 2003 $o_rldcs%name,itau_w,lwdn0) 2004 ENDIF 2005 2006 IF (o_tnt%flag(iff)<=lev_files(iff)) THEN 2007 zx_tmp_fi3d(1:klon,1:klev)=d_t(1:klon,1:klev)+ 2008 $d_t_dyn(1:klon,1:klev) 2009 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 2010 $o_tnt%name,itau_w,zx_tmp_fi3d) 2011 ENDIF 2012 2013 IF (o_tntr%flag(iff)<=lev_files(iff)) THEN 2014 zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY - 2015 $cool(1:klon,1:klev)/RDAY 2016 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 2017 $o_tntr%name,itau_w,zx_tmp_fi3d) 2018 ENDIF 2019 2020 IF (o_tntscpbl%flag(iff)<=lev_files(iff)) THEN 2021 zx_tmp_fi3d(1:klon,1:klev)= (d_t_lsc(1:klon,1:klev)+ 2022 $ d_t_eva(1:klon,1:klev)+ 2023 $ d_t_vdf(1:klon,1:klev))/pdtphys 2024 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 2025 $o_tntscpbl%name,itau_w,zx_tmp_fi3d) 2026 ENDIF 2027 2028 IF (o_tnhus%flag(iff)<=lev_files(iff)) THEN 2029 zx_tmp_fi3d(1:klon,1:klev)=d_qx(1:klon,1:klev,ivap)+ 2030 $d_q_dyn(1:klon,1:klev) 2031 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 2032 $o_tnhus%name,itau_w,zx_tmp_fi3d) 2033 ENDIF 2034 2035 IF (o_tnhusscpbl%flag(iff)<=lev_files(iff)) THEN 2036 zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys+ 2037 $ d_q_eva(1:klon,1:klev)/pdtphys 2038 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 2039 $o_tnhusscpbl%name,itau_w,zx_tmp_fi3d) 2040 ENDIF 2041 2042 IF (o_evu%flag(iff)<=lev_files(iff)) THEN 2043 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 2044 $o_evu%name,itau_w,coefm) 2045 ENDIF 2046 2047 IF (o_h2o%flag(iff)<=lev_files(iff)) THEN 2048 zx_tmp_fi3d(1:klon,1:klev)=q_seri(1:klon,1:klev)+ 2049 $ ql_seri(1:klon,1:klev) 2050 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 2051 $o_h2o%name,itau_w,zx_tmp_fi3d) 2052 ENDIF 2053 2054 IF (o_mcd%flag(iff)<=lev_files(iff)) THEN 2055 zx_tmp_fi3d(1:klon,1:klev)=-1 * (dnwd(1:klon,1:klev)+ 2056 $ dnwd0(1:klon,1:klev)) 2057 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 2058 $o_mcd%name,itau_w,zx_tmp_fi3d) 2059 ENDIF 2060 2061 IF (o_dmc%flag(iff)<=lev_files(iff)) THEN 2062 zx_tmp_fi3d(1:klon,1:klev)=upwd(1:klon,1:klev) + 2063 $ dnwd(1:klon,1:klev)+ dnwd0(1:klon,1:klev) 2064 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 2065 $o_dmc%name,itau_w,zx_tmp_fi3d) 2066 ENDIF 2067 2068 IF (o_ref_liq%flag(iff)<=lev_files(iff)) THEN 2069 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 2070 $o_ref_liq%name,itau_w,ref_liq) 2071 ENDIF 2072 2073 IF (o_ref_ice%flag(iff)<=lev_files(iff)) THEN 2074 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 2075 $o_ref_ice%name,itau_w,ref_ice) 2076 ENDIF 2077 2078 if (RCO2_per.NE.RCO2_act.OR.RCH4_per.NE.RCH4_act.OR. 2079 $ RN2O_per.NE.RN2O_act.OR.RCFC11_per.NE.RCFC11_act.OR. 2080 $ RCFC12_per.NE.RCFC12_act) THEN 2081 2082 IF (o_rsut4co2%flag(iff)<=lev_files(iff)) THEN 2083 zx_tmp_fi2d(1 : klon) = swupp ( 1 : klon, klevp1 ) 2084 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 2085 $o_rsut4co2%name,itau_w,zx_tmp_fi2d) 2086 ENDIF 2087 2088 IF (o_rlut4co2%flag(iff)<=lev_files(iff)) THEN 2089 zx_tmp_fi2d(1 : klon) = lwupp ( 1 : klon, klevp1 ) 2090 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 2091 $o_rlut4co2%name,itau_w,zx_tmp_fi2d) 2092 ENDIF 2093 2094 IF (o_rsutcs4co2%flag(iff)<=lev_files(iff)) THEN 2095 zx_tmp_fi2d(1 : klon) = swup0p ( 1 : klon, klevp1 ) 2096 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 2097 $o_rsutcs4co2%name,itau_w,zx_tmp_fi2d) 2098 ENDIF 2099 2100 IF (o_rlutcs4co2%flag(iff)<=lev_files(iff)) THEN 2101 zx_tmp_fi2d(1 : klon) = lwup0p ( 1 : klon, klevp1 ) 2102 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 2103 $o_rlutcs4co2%name,itau_w,zx_tmp_fi2d) 2104 ENDIF 2105 2106 IF (o_rsu4co2%flag(iff)<=lev_files(iff)) THEN 2107 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 2108 $o_rsu4co2%name,itau_w,swupp) 2109 ENDIF 2110 2111 IF (o_rlu4co2%flag(iff)<=lev_files(iff)) THEN 2112 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 2113 $o_rlu4co2%name,itau_w,lwupp) 2114 ENDIF 2115 2116 IF (o_rsucs4co2%flag(iff)<=lev_files(iff)) THEN 2117 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 2118 $o_rsucs4co2%name,itau_w,swup0p) 2119 ENDIF 2120 2121 IF (o_rlucs4co2%flag(iff)<=lev_files(iff)) THEN 2122 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 2123 $o_rlucs4co2%name,itau_w,lwup0p) 2124 ENDIF 2125 2126 IF (o_rsd4co2%flag(iff)<=lev_files(iff)) THEN 2127 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 2128 $o_rsd4co2%name,itau_w,swdnp) 2129 ENDIF 2130 2131 IF (o_rld4co2%flag(iff)<=lev_files(iff)) THEN 2132 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 2133 $o_rld4co2%name,itau_w,lwdnp) 2134 ENDIF 2135 2136 IF (o_rsdcs4co2%flag(iff)<=lev_files(iff)) THEN 2137 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 2138 $o_rsdcs4co2%name,itau_w,swdn0p) 2139 ENDIF 2140 2141 IF (o_rldcs4co2%flag(iff)<=lev_files(iff)) THEN 2142 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 2143 $o_rldcs4co2%name,itau_w,lwdn0p) 2144 ENDIF 2145 2146 endif 1589 2147 1590 2148 if (nqtot.GE.3) THEN 1591 2149 DO iq=3,nqtot 1592 2150 IF (o_trac(iq-2)%flag(iff)<=lev_files(iff)) THEN 1593 CALL histwrite_phy(nid_files(iff), 2151 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 1594 2152 s o_trac(iq-2)%name,itau_w,qx(:,:,iq)) 1595 2153 ENDIF … … 1599 2157 if (ok_sync) then 1600 2158 c$OMP MASTER 1601 2159 call histsync(nid_files(iff)) 1602 2160 c$OMP END MASTER 1603 2161 endif
Note: See TracChangeset
for help on using the changeset viewer.
