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