Changeset 1419 for LMDZ4/branches/LMDZ4_AR5
- Timestamp:
- Jul 22, 2010, 5:19:06 PM (14 years ago)
- Location:
- LMDZ4/branches/LMDZ4_AR5/libf/phylmd
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/branches/LMDZ4_AR5/libf/phylmd/ini_histhfNMC.h
r1400 r1419 88 88 . "ave(X)", zout,zout) 89 89 c 90 IF (1.EQ.0) THEN 91 c 90 92 CALL histdef(nid_hfnmc, "psbg", 91 93 . "Pressure sfce below ground","%", … … 146 148 . "ave(X)", zout,zout) 147 149 c 150 ENDIF !(1.EQ.0) THEN 148 151 CALL histend(nid_hfnmc) 149 152 c -
LMDZ4/branches/LMDZ4_AR5/libf/phylmd/phys_output_mod.F90
r1405 r1419 66 66 type(ctrl_out),save :: o_t2m_min = ctrl_out((/ 1, 1, 10, 10, 10 /),'t2m_min') 67 67 type(ctrl_out),save :: o_t2m_max = ctrl_out((/ 1, 1, 10, 10, 10 /),'t2m_max') 68 type(ctrl_out),save,dimension(4) :: o_t2m_srf = (/ ctrl_out((/ 10, 4, 10, 10, 10 /),'t2m_ter'), &69 ctrl_out((/ 10, 4, 10, 10, 10 /),'t2m_lic'), &70 ctrl_out((/ 10, 4, 10, 10, 10 /),'t2m_oce'), &71 ctrl_out((/ 10, 4, 10, 10, 10 /),'t2m_sic') /)68 type(ctrl_out),save,dimension(4) :: o_t2m_srf = (/ ctrl_out((/ 10, 6, 10, 10, 10 /),'t2m_ter'), & 69 ctrl_out((/ 10, 6, 10, 10, 10 /),'t2m_lic'), & 70 ctrl_out((/ 10, 6, 10, 10, 10 /),'t2m_oce'), & 71 ctrl_out((/ 10, 6, 10, 10, 10 /),'t2m_sic') /) 72 72 73 73 type(ctrl_out),save :: o_wind10m = ctrl_out((/ 1, 1, 1, 10, 10 /),'wind10m') … … 80 80 type(ctrl_out),save :: o_qsurf = ctrl_out((/ 1, 10, 10, 10, 10 /),'qsurf') 81 81 82 type(ctrl_out),save,dimension(4) :: o_u10m_srf = (/ ctrl_out((/ 10, 4, 10, 10, 10 /),'u10m_ter'), &83 ctrl_out((/ 10, 4, 10, 10, 10 /),'u10m_lic'), &84 ctrl_out((/ 10, 4, 10, 10, 10 /),'u10m_oce'), &85 ctrl_out((/ 10, 4, 10, 10, 10 /),'u10m_sic') /)86 87 type(ctrl_out),save,dimension(4) :: o_v10m_srf = (/ ctrl_out((/ 10, 4, 10, 10, 10 /),'v10m_ter'), &88 ctrl_out((/ 10, 4, 10, 10, 10 /),'v10m_lic'), &89 ctrl_out((/ 10, 4, 10, 10, 10 /),'v10m_oce'), &90 ctrl_out((/ 10, 4, 10, 10, 10 /),'v10m_sic') /)82 type(ctrl_out),save,dimension(4) :: o_u10m_srf = (/ ctrl_out((/ 10, 6, 10, 10, 10 /),'u10m_ter'), & 83 ctrl_out((/ 10, 6, 10, 10, 10 /),'u10m_lic'), & 84 ctrl_out((/ 10, 6, 10, 10, 10 /),'u10m_oce'), & 85 ctrl_out((/ 10, 6, 10, 10, 10 /),'u10m_sic') /) 86 87 type(ctrl_out),save,dimension(4) :: o_v10m_srf = (/ ctrl_out((/ 10, 6, 10, 10, 10 /),'v10m_ter'), & 88 ctrl_out((/ 10, 6, 10, 10, 10 /),'v10m_lic'), & 89 ctrl_out((/ 10, 6, 10, 10, 10 /),'v10m_oce'), & 90 ctrl_out((/ 10, 6, 10, 10, 10 /),'v10m_sic') /) 91 91 92 92 type(ctrl_out),save :: o_qsol = ctrl_out((/ 1, 10, 10, 10, 10 /),'qsol') … … 99 99 type(ctrl_out),save :: o_snow = ctrl_out((/ 1, 1, 10, 10, 5 /),'snow') 100 100 type(ctrl_out),save :: o_evap = ctrl_out((/ 1, 1, 10, 10, 10 /),'evap') 101 type(ctrl_out),save,dimension(4) :: o_evap_srf = (/ ctrl_out((/ 1, 1, 10, 10, 10 /),'evap_ter'), &102 ctrl_out((/ 1, 1, 10, 10, 10 /),'evap_lic'), &103 ctrl_out((/ 1, 1, 10, 10, 10 /),'evap_oce'), &104 ctrl_out((/ 1, 1, 10, 10, 10 /),'evap_sic') /)101 type(ctrl_out),save,dimension(4) :: o_evap_srf = (/ ctrl_out((/ 1, 6, 10, 10, 10 /),'evap_ter'), & 102 ctrl_out((/ 1, 6, 10, 10, 10 /),'evap_lic'), & 103 ctrl_out((/ 1, 6, 10, 10, 10 /),'evap_oce'), & 104 ctrl_out((/ 1, 6, 10, 10, 10 /),'evap_sic') /) 105 105 type(ctrl_out),save :: o_msnow = ctrl_out((/ 1, 10, 10, 10, 10 /),'msnow') 106 106 type(ctrl_out),save :: o_fsnow = ctrl_out((/ 1, 10, 10, 10, 10 /),'fsnow') … … 135 135 type(ctrl_out),save :: o_soll = ctrl_out((/ 1, 1, 10, 10, 10 /),'soll') 136 136 type(ctrl_out),save :: o_soll0 = ctrl_out((/ 1, 5, 10, 10, 10 /),'soll0') 137 type(ctrl_out),save :: o_radsol = ctrl_out((/ 1, 1, 10, 10, 10 /),'radsol')137 type(ctrl_out),save :: o_radsol = ctrl_out((/ 1, 7, 10, 10, 10 /),'radsol') 138 138 type(ctrl_out),save :: o_SWupSFC = ctrl_out((/ 1, 4, 10, 10, 5 /),'SWupSFC') 139 139 type(ctrl_out),save :: o_SWupSFCclr = ctrl_out((/ 1, 4, 10, 10, 5 /),'SWupSFCclr') … … 153 153 type(ctrl_out),save :: o_taux = ctrl_out((/ 1, 10, 10, 10, 10 /),'taux') 154 154 type(ctrl_out),save :: o_tauy = ctrl_out((/ 1, 10, 10, 10, 10 /),'tauy') 155 type(ctrl_out),save,dimension(4) :: o_taux_srf = (/ ctrl_out((/ 1, 4, 10, 10, 10 /),'taux_ter'), &156 ctrl_out((/ 1, 4, 10, 10, 10 /),'taux_lic'), &157 ctrl_out((/ 1, 4, 10, 10, 10 /),'taux_oce'), &158 ctrl_out((/ 1, 4, 10, 10, 10 /),'taux_sic') /)159 160 type(ctrl_out),save,dimension(4) :: o_tauy_srf = (/ ctrl_out((/ 1, 4, 10, 10, 10 /),'tauy_ter'), &161 ctrl_out((/ 1, 4, 10, 10, 10 /),'tauy_lic'), &162 ctrl_out((/ 1, 4, 10, 10, 10 /),'tauy_oce'), &163 ctrl_out((/ 1, 4, 10, 10, 10 /),'tauy_sic') /)164 165 166 type(ctrl_out),save,dimension(4) :: o_pourc_srf = (/ ctrl_out((/ 1, 4, 10, 10, 10 /),'pourc_ter'), &167 ctrl_out((/ 1, 4, 10, 10, 10 /),'pourc_lic'), &168 ctrl_out((/ 1, 4, 10, 10, 10 /),'pourc_oce'), &169 ctrl_out((/ 1, 4, 10, 10, 10 /),'pourc_sic') /)170 171 type(ctrl_out),save,dimension(4) :: o_fract_srf = (/ ctrl_out((/ 1, 4, 10, 10, 10 /),'fract_ter'), &172 ctrl_out((/ 1, 4, 10, 10, 10 /),'fract_lic'), &173 ctrl_out((/ 1, 4, 10, 10, 10 /),'fract_oce'), &174 ctrl_out((/ 1, 4, 10, 10, 10 /),'fract_sic') /)175 176 type(ctrl_out),save,dimension(4) :: o_tsol_srf = (/ ctrl_out((/ 1, 4, 10, 10, 10 /),'tsol_ter'), &177 ctrl_out((/ 1, 4, 10, 10, 10 /),'tsol_lic'), &178 ctrl_out((/ 1, 4, 10, 10, 10 /),'tsol_oce'), &179 ctrl_out((/ 1, 4, 10, 10, 10 /),'tsol_sic') /)180 181 type(ctrl_out),save,dimension(4) :: o_sens_srf = (/ ctrl_out((/ 1, 4, 10, 1, 10 /),'sens_ter'), &182 ctrl_out((/ 1, 4, 10, 1, 10 /),'sens_lic'), &183 ctrl_out((/ 1, 4, 10, 1, 10 /),'sens_oce'), &184 ctrl_out((/ 1, 4, 10, 1, 10 /),'sens_sic') /)185 186 type(ctrl_out),save,dimension(4) :: o_lat_srf = (/ ctrl_out((/ 1, 4, 10, 1, 10 /),'lat_ter'), &187 ctrl_out((/ 1, 4, 10, 1, 10 /),'lat_lic'), &188 ctrl_out((/ 1, 4, 10, 1, 10 /),'lat_oce'), &189 ctrl_out((/ 1, 4, 10, 1, 10 /),'lat_sic') /)155 type(ctrl_out),save,dimension(4) :: o_taux_srf = (/ ctrl_out((/ 1, 6, 10, 10, 10 /),'taux_ter'), & 156 ctrl_out((/ 1, 6, 10, 10, 10 /),'taux_lic'), & 157 ctrl_out((/ 1, 6, 10, 10, 10 /),'taux_oce'), & 158 ctrl_out((/ 1, 6, 10, 10, 10 /),'taux_sic') /) 159 160 type(ctrl_out),save,dimension(4) :: o_tauy_srf = (/ ctrl_out((/ 1, 6, 10, 10, 10 /),'tauy_ter'), & 161 ctrl_out((/ 1, 6, 10, 10, 10 /),'tauy_lic'), & 162 ctrl_out((/ 1, 6, 10, 10, 10 /),'tauy_oce'), & 163 ctrl_out((/ 1, 6, 10, 10, 10 /),'tauy_sic') /) 164 165 166 type(ctrl_out),save,dimension(4) :: o_pourc_srf = (/ ctrl_out((/ 1, 7, 10, 10, 10 /),'pourc_ter'), & 167 ctrl_out((/ 1, 7, 10, 10, 10 /),'pourc_lic'), & 168 ctrl_out((/ 1, 7, 10, 10, 10 /),'pourc_oce'), & 169 ctrl_out((/ 1, 7, 10, 10, 10 /),'pourc_sic') /) 170 171 type(ctrl_out),save,dimension(4) :: o_fract_srf = (/ ctrl_out((/ 1, 6, 10, 10, 10 /),'fract_ter'), & 172 ctrl_out((/ 1, 6, 10, 10, 10 /),'fract_lic'), & 173 ctrl_out((/ 1, 6, 10, 10, 10 /),'fract_oce'), & 174 ctrl_out((/ 1, 6, 10, 10, 10 /),'fract_sic') /) 175 176 type(ctrl_out),save,dimension(4) :: o_tsol_srf = (/ ctrl_out((/ 1, 6, 10, 10, 10 /),'tsol_ter'), & 177 ctrl_out((/ 1, 6, 10, 10, 10 /),'tsol_lic'), & 178 ctrl_out((/ 1, 6, 10, 10, 10 /),'tsol_oce'), & 179 ctrl_out((/ 1, 6, 10, 10, 10 /),'tsol_sic') /) 180 181 type(ctrl_out),save,dimension(4) :: o_sens_srf = (/ ctrl_out((/ 1, 6, 10, 7, 10 /),'sens_ter'), & 182 ctrl_out((/ 1, 6, 10, 7, 10 /),'sens_lic'), & 183 ctrl_out((/ 1, 6, 10, 7, 10 /),'sens_oce'), & 184 ctrl_out((/ 1, 6, 10, 7, 10 /),'sens_sic') /) 185 186 type(ctrl_out),save,dimension(4) :: o_lat_srf = (/ ctrl_out((/ 1, 6, 10, 7, 10 /),'lat_ter'), & 187 ctrl_out((/ 1, 6, 10, 7, 10 /),'lat_lic'), & 188 ctrl_out((/ 1, 6, 10, 7, 10 /),'lat_oce'), & 189 ctrl_out((/ 1, 6, 10, 7, 10 /),'lat_sic') /) 190 190 191 191 type(ctrl_out),save,dimension(4) :: o_flw_srf = (/ ctrl_out((/ 1, 10, 10, 10, 10 /),'flw_ter'), & … … 211 211 212 212 type(ctrl_out),save :: o_cdrm = ctrl_out((/ 1, 10, 10, 10, 10 /),'cdrm') 213 type(ctrl_out),save :: o_cdrh = ctrl_out((/ 1, 10, 10, 1, 10 /),'cdrh')213 type(ctrl_out),save :: o_cdrh = ctrl_out((/ 1, 10, 10, 7, 10 /),'cdrh') 214 214 type(ctrl_out),save :: o_cldl = ctrl_out((/ 1, 1, 10, 10, 10 /),'cldl') 215 215 type(ctrl_out),save :: o_cldm = ctrl_out((/ 1, 1, 10, 10, 10 /),'cldm') … … 257 257 ! Champs interpolles sur des niveaux de pression ??? a faire correctement 258 258 259 type(ctrl_out),save,dimension(6) :: o_uSTDlevs = (/ ctrl_out((/ 1, 1, 3, 10, 10 /),'u850'), &260 ctrl_out((/ 1, 1, 3, 10, 10 /),'u700'), &261 ctrl_out((/ 1, 1, 3, 10, 10 /),'u500'), &262 ctrl_out((/ 1, 1, 3, 10, 10 /),'u200'), &263 ctrl_out((/ 1, 1, 3, 10, 10 /),'u50'), &264 ctrl_out((/ 1, 1, 3, 10, 10 /),'u10') /)259 type(ctrl_out),save,dimension(6) :: o_uSTDlevs = (/ ctrl_out((/ 1, 7, 7, 10, 10 /),'u850'), & 260 ctrl_out((/ 1, 7, 7, 10, 10 /),'u700'), & 261 ctrl_out((/ 1, 7, 7, 10, 10 /),'u500'), & 262 ctrl_out((/ 1, 7, 7, 10, 10 /),'u200'), & 263 ctrl_out((/ 1, 7, 7, 10, 10 /),'u50'), & 264 ctrl_out((/ 1, 7, 7, 10, 10 /),'u10') /) 265 265 266 266 267 type(ctrl_out),save,dimension(6) :: o_vSTDlevs = (/ ctrl_out((/ 1, 1, 3, 10, 10 /),'v850'), &268 ctrl_out((/ 1, 1, 3, 10, 10 /),'v700'), &269 ctrl_out((/ 1, 1, 3, 10, 10 /),'v500'), &270 ctrl_out((/ 1, 1, 3, 10, 10 /),'v200'), &271 ctrl_out((/ 1, 1, 3, 10, 10 /),'v50'), &272 ctrl_out((/ 1, 1, 3, 10, 10 /),'v10') /)273 274 type(ctrl_out),save,dimension(6) :: o_wSTDlevs = (/ ctrl_out((/ 1, 1, 3, 10, 10 /),'w850'), &275 ctrl_out((/ 1, 1, 3, 10, 10 /),'w700'), &276 ctrl_out((/ 1, 1, 3, 10, 10 /),'w500'), &277 ctrl_out((/ 1, 1, 3, 10, 10 /),'w200'), &278 ctrl_out((/ 1, 1, 3, 10, 10 /),'w50'), &279 ctrl_out((/ 1, 1, 3, 10, 10 /),'w10') /)280 281 type(ctrl_out),save,dimension(6) :: o_tSTDlevs = (/ ctrl_out((/ 1, 1, 3, 10, 10 /),'t850'), &282 ctrl_out((/ 1, 1, 3, 10, 10 /),'t700'), &283 ctrl_out((/ 1, 1, 3, 10, 10 /),'t500'), &284 ctrl_out((/ 1, 1, 3, 10, 10 /),'t200'), &285 ctrl_out((/ 1, 1, 3, 10, 10 /),'t50'), &286 ctrl_out((/ 1, 1, 3, 10, 10 /),'t10') /)287 288 type(ctrl_out),save,dimension(6) :: o_qSTDlevs = (/ ctrl_out((/ 1, 1, 3, 10, 10 /),'q850'), &289 ctrl_out((/ 1, 1, 3, 10, 10 /),'q700'), &290 ctrl_out((/ 1, 1, 3, 10, 10 /),'q500'), &291 ctrl_out((/ 1, 1, 3, 10, 10 /),'q200'), &292 ctrl_out((/ 1, 1, 3, 10, 10 /),'q50'), &293 ctrl_out((/ 1, 1, 3, 10, 10 /),'q10') /)294 295 type(ctrl_out),save,dimension(6) :: o_zSTDlevs = (/ ctrl_out((/ 1, 1, 3, 10, 10 /),'z850'), &296 ctrl_out((/ 1, 1, 3, 10, 10 /),'z700'), &297 ctrl_out((/ 1, 1, 3, 10, 10 /),'z500'), &298 ctrl_out((/ 1, 1, 3, 10, 10 /),'z200'), &299 ctrl_out((/ 1, 1, 3, 10, 10 /),'z50'), &300 ctrl_out((/ 1, 1, 3, 10, 10 /),'z10') /)267 type(ctrl_out),save,dimension(6) :: o_vSTDlevs = (/ ctrl_out((/ 1, 7, 7, 10, 10 /),'v850'), & 268 ctrl_out((/ 1, 7, 7, 10, 10 /),'v700'), & 269 ctrl_out((/ 1, 7, 7, 10, 10 /),'v500'), & 270 ctrl_out((/ 1, 7, 7, 10, 10 /),'v200'), & 271 ctrl_out((/ 1, 7, 7, 10, 10 /),'v50'), & 272 ctrl_out((/ 1, 7, 7, 10, 10 /),'v10') /) 273 274 type(ctrl_out),save,dimension(6) :: o_wSTDlevs = (/ ctrl_out((/ 1, 7, 7, 10, 10 /),'w850'), & 275 ctrl_out((/ 1, 7, 7, 10, 10 /),'w700'), & 276 ctrl_out((/ 1, 7, 7, 10, 10 /),'w500'), & 277 ctrl_out((/ 1, 7, 7, 10, 10 /),'w200'), & 278 ctrl_out((/ 1, 7, 7, 10, 10 /),'w50'), & 279 ctrl_out((/ 1, 7, 7, 10, 10 /),'w10') /) 280 281 type(ctrl_out),save,dimension(6) :: o_tSTDlevs = (/ ctrl_out((/ 1, 7, 7, 10, 10 /),'t850'), & 282 ctrl_out((/ 1, 7, 7, 10, 10 /),'t700'), & 283 ctrl_out((/ 1, 7, 7, 10, 10 /),'t500'), & 284 ctrl_out((/ 1, 7, 7, 10, 10 /),'t200'), & 285 ctrl_out((/ 1, 7, 7, 10, 10 /),'t50'), & 286 ctrl_out((/ 1, 7, 7, 10, 10 /),'t10') /) 287 288 type(ctrl_out),save,dimension(6) :: o_qSTDlevs = (/ ctrl_out((/ 1, 7, 7, 10, 10 /),'q850'), & 289 ctrl_out((/ 1, 7, 7, 10, 10 /),'q700'), & 290 ctrl_out((/ 1, 7, 7, 10, 10 /),'q500'), & 291 ctrl_out((/ 1, 7, 7, 10, 10 /),'q200'), & 292 ctrl_out((/ 1, 7, 7, 10, 10 /),'q50'), & 293 ctrl_out((/ 1, 7, 7, 10, 10 /),'q10') /) 294 295 type(ctrl_out),save,dimension(6) :: o_zSTDlevs = (/ ctrl_out((/ 1, 7, 7, 10, 10 /),'z850'), & 296 ctrl_out((/ 1, 7, 7, 10, 10 /),'z700'), & 297 ctrl_out((/ 1, 7, 7, 10, 10 /),'z500'), & 298 ctrl_out((/ 1, 7, 7, 10, 10 /),'z200'), & 299 ctrl_out((/ 1, 7, 7, 10, 10 /),'z50'), & 300 ctrl_out((/ 1, 7, 7, 10, 10 /),'z10') /) 301 301 302 302 … … 419 419 type(ctrl_out),save :: o_theta = ctrl_out((/ 2, 3, 4, 10, 10 /),'theta') 420 420 type(ctrl_out),save :: o_ovap = ctrl_out((/ 2, 3, 4, 10, 10 /),'ovap') 421 type(ctrl_out),save :: o_ovapinit = ctrl_out((/ 2, 3, 10, 10, 10 /),'ovapinit')421 type(ctrl_out),save :: o_ovapinit = ctrl_out((/ 2, 10, 10, 10, 10 /),'ovapinit') 422 422 type(ctrl_out),save :: o_wvapp = ctrl_out((/ 2, 10, 10, 10, 10 /),'wvapp') 423 423 type(ctrl_out),save :: o_geop = ctrl_out((/ 2, 3, 10, 10, 10 /),'geop') … … 449 449 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 450 450 451 type(ctrl_out),save,dimension(4) :: o_albe_srf = (/ ctrl_out((/ 3, 4, 10, 1, 10 /),'albe_ter'), &452 ctrl_out((/ 3, 4, 10, 1, 10 /),'albe_lic'), &453 ctrl_out((/ 3, 4, 10, 1, 10 /),'albe_oce'), &454 ctrl_out((/ 3, 4, 10, 1, 10 /),'albe_sic') /)451 type(ctrl_out),save,dimension(4) :: o_albe_srf = (/ ctrl_out((/ 3, 7, 10, 7, 10 /),'albe_ter'), & 452 ctrl_out((/ 3, 7, 10, 7, 10 /),'albe_lic'), & 453 ctrl_out((/ 3, 7, 10, 7, 10 /),'albe_oce'), & 454 ctrl_out((/ 3, 7, 10, 7, 10 /),'albe_sic') /) 455 455 456 456 type(ctrl_out),save,dimension(4) :: o_ages_srf = (/ ctrl_out((/ 10, 10, 10, 10, 10 /),'ages_ter'), & … … 459 459 ctrl_out((/ 3, 10, 10, 10, 10 /),'ages_sic') /) 460 460 461 type(ctrl_out),save,dimension(4) :: o_rugs_srf = (/ ctrl_out((/ 3, 4, 10, 10, 10 /),'rugs_ter'), &462 ctrl_out((/ 3, 4, 10, 10, 10 /),'rugs_lic'), &463 ctrl_out((/ 3, 4, 10, 10, 10 /),'rugs_oce'), &464 ctrl_out((/ 3, 4, 10, 10, 10 /),'rugs_sic') /)461 type(ctrl_out),save,dimension(4) :: o_rugs_srf = (/ ctrl_out((/ 3, 6, 10, 10, 10 /),'rugs_ter'), & 462 ctrl_out((/ 3, 6, 10, 10, 10 /),'rugs_lic'), & 463 ctrl_out((/ 3, 6, 10, 10, 10 /),'rugs_oce'), & 464 ctrl_out((/ 3, 6, 10, 10, 10 /),'rugs_sic') /) 465 465 466 466 type(ctrl_out),save :: o_alb1 = ctrl_out((/ 3, 10, 10, 10, 10 /),'alb1') … … 477 477 type(ctrl_out),save :: o_dudyn = ctrl_out((/ 4, 10, 10, 10, 10 /),'dudyn') !AXC 478 478 type(ctrl_out),save :: o_dvdyn = ctrl_out((/ 4, 10, 10, 10, 10 /),'dvdyn') !AXC 479 type(ctrl_out),save :: o_dtcon = ctrl_out((/ 4, 5, 10, 10, 10 /),'dtcon')479 type(ctrl_out),save :: o_dtcon = ctrl_out((/ 4, 10, 10, 10, 10 /),'dtcon') 480 480 type(ctrl_out),save :: o_ducon = ctrl_out((/ 4, 10, 10, 10, 10 /),'ducon') 481 type(ctrl_out),save :: o_dqcon = ctrl_out((/ 4, 5, 10, 10, 10 /),'dqcon')481 type(ctrl_out),save :: o_dqcon = ctrl_out((/ 4, 10, 10, 10, 10 /),'dqcon') 482 482 type(ctrl_out),save :: o_dtwak = ctrl_out((/ 4, 5, 10, 10, 10 /),'dtwak') 483 483 type(ctrl_out),save :: o_dqwak = ctrl_out((/ 4, 5, 10, 10, 10 /),'dqwak') … … 524 524 type(ctrl_out),save :: o_dulif = ctrl_out((/ 4, 10, 10, 10, 10 /),'dulif') 525 525 type(ctrl_out),save :: o_dvlif = ctrl_out((/ 4, 10, 10, 10, 10 /),'dvlif') 526 type(ctrl_out),save :: o_duhin = ctrl_out((/ 4, 10, 10, 10, 10 /),'duhin') 527 type(ctrl_out),save :: o_dvhin = ctrl_out((/ 4, 10, 10, 10, 10 /),'dvhin') 526 528 527 529 ! Attention a refaire correctement … … 1224 1226 ENDIF 1225 1227 1228 IF (ok_hines) then 1229 CALL histdef3d(iff,o_duhin%flag,o_duhin%name, "Hines GWD dU", "m/s2") 1230 CALL histdef3d(iff,o_dvhin%flag,o_dvhin%name, "Hines GWD dV", "m/s2") 1231 ENDIF 1232 1226 1233 1227 1234 !Attention : sorties uniquement pour traceurs 3 et 4 -
LMDZ4/branches/LMDZ4_AR5/libf/phylmd/phys_output_write.h
r1405 r1419 1509 1509 ENDIF 1510 1510 1511 IF (ok_hines) THEN 1512 IF (o_duhin%flag(iff)<=lev_files(iff)) THEN 1513 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) 1515 ENDIF 1516 1517 IF (o_dvhin%flag(iff)<=lev_files(iff)) THEN 1518 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 1511 1522 1512 1523 IF (nqtot.GE.3 .AND. o_trac(1)%flag(iff)<=lev_files(iff)) THEN -
LMDZ4/branches/LMDZ4_AR5/libf/phylmd/write_histhfNMC.h
r1398 r1419 30 30 CALL histwrite_phy(nid_hfnmc,"wap",itau_w,wwriteSTD3) 31 31 c 32 IF (1.EQ.0) THEN 33 c 32 34 DO k=1, nlevSTD 33 35 DO i=1, klon … … 60 62 CALL histwrite_phy(nid_hfnmc,"TxT",itau_w,T2sumSTD(:,:,3)) 61 63 c 62 c ENDIF !type_run 64 ENDIF !(1.EQ.0) THEN 63 65 c 64 66 if (ok_sync) then
Note: See TracChangeset
for help on using the changeset viewer.