Changeset 1562 for LMDZ5/trunk/libf
- Timestamp:
- Aug 25, 2011, 11:55:45 AM (13 years ago)
- Location:
- LMDZ5/trunk/libf/phylmd
- Files:
-
- 3 edited
-
YOEGWD.h (modified) (1 diff)
-
phys_output_mod.F90 (modified) (20 diffs)
-
phys_output_write.h (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/phylmd/YOEGWD.h
r776 r1562 2 2 ! $Header$ 3 3 ! 4 C ----------------------------------------------------------------- 5 C* *COMMON* *YOEGWD* - PARAMETERS FOR GRAVITY WAVE DRAG CALCULATIONS 6 C ----------------------------------------------------------------- 7 C 4 ! ATTENTION : ce fichier include est compatible format fixe/format libre 5 ! veillez n'utiliser que des ! pour les commentaires 6 ! et bien positionner les & des lignes de continuation 7 ! (les placer en colonne 6 et en colonne 73) 8 ! ----------------------------------------------------------------- 9 !* *COMMON* *YOEGWD* - PARAMETERS FOR GRAVITY WAVE DRAG CALCULATIONS 10 ! ----------------------------------------------------------------- 11 ! 8 12 integer NKTOPG,NSTRA 9 13 real GFRCRIT,GKWAKE,GRCRIT,GVCRIT,GKDRAG,GKLIFT 10 14 real GHMAX,GRAHILO,GSIGCR,GSSEC,GTSEC,GVSEC 11 COMMON/YOEGWD/ GFRCRIT,GKWAKE,GRCRIT,GVCRIT,GKDRAG,GKLIFT 12 * ,GHMAX,GRAHILO,GSIGCR,NKTOPG,NSTRA,GSSEC,GTSEC,GVSEC 13 c$OMP THREADPRIVATE(/YOEGWD/) 14 C 15 16 15 COMMON/YOEGWD/ GFRCRIT,GKWAKE,GRCRIT,GVCRIT,GKDRAG,GKLIFT & 16 & ,GHMAX,GRAHILO,GSIGCR,NKTOPG,NSTRA,GSSEC,GTSEC,GVSEC 17 !$OMP THREADPRIVATE(/YOEGWD/) -
LMDZ5/trunk/libf/phylmd/phys_output_mod.F90
r1560 r1562 17 17 18 18 19 20 21 22 23 24 19 integer, parameter :: nfiles = 6 20 logical, dimension(nfiles), save :: clef_files 21 logical, dimension(nfiles), save :: clef_stations 22 integer, dimension(nfiles), save :: lev_files 23 integer, dimension(nfiles), save :: nid_files 24 integer, dimension(nfiles), save :: nnid_files 25 25 !!$OMP THREADPRIVATE(clef_files, clef_stations, lev_files,nid_files,nnid_files) 26 27 28 29 30 ! integer, dimension(nfiles), private, save :: nvertp031 32 33 34 !$OMP THREADPRIVATE(nhorim, nvertm, zoutm,zdtime,type_ecri)35 36 ! integer, save :: nid_hf3d26 integer, dimension(nfiles), private, save :: nnhorim 27 28 integer, dimension(nfiles), private, save :: nhorim, nvertm 29 integer, dimension(nfiles), private, save :: nvertap, nvertbp, nvertAlt 30 ! integer, dimension(nfiles), private, save :: nvertp0 31 real, dimension(nfiles), private, save :: zoutm 32 real, private, save :: zdtime 33 CHARACTER(len=20), dimension(nfiles), private, save :: type_ecri 34 !$OMP THREADPRIVATE(nhorim, nvertm, zoutm,zdtime,type_ecri) 35 36 ! integer, save :: nid_hf3d 37 37 38 38 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 39 !! Definition pour chaque variable du niveau d ecriture dans chaque fichier39 !! Definition pour chaque variable du niveau d ecriture dans chaque fichier 40 40 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!/ histmth, histday, histhf, histins /),'!!!!!!!!!!!! 41 41 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! … … 45 45 46 46 TYPE ctrl_out 47 integer,dimension(6) :: flag48 character(len=20) :: name47 integer,dimension(6) :: flag 48 character(len=20) :: name 49 49 END TYPE ctrl_out 50 50 … … 61 61 type(ctrl_out),save :: o_contfracOR = ctrl_out((/ 10, 1, 1, 10, 10, 10 /),'contfracOR') 62 62 type(ctrl_out),save :: o_aireTER = ctrl_out((/ 10, 10, 1, 10, 10, 10 /),'aireTER') 63 63 64 64 !!! 2D 65 65 type(ctrl_out),save :: o_flat = ctrl_out((/ 5, 1, 10, 10, 5, 10 /),'flat') … … 70 70 type(ctrl_out),save :: o_t2m_max = ctrl_out((/ 1, 1, 10, 10, 10, 10 /),'t2m_max') 71 71 type(ctrl_out),save,dimension(4) :: o_t2m_srf = (/ ctrl_out((/ 10, 6, 10, 10, 10, 10 /),'t2m_ter'), & 72 73 74 72 ctrl_out((/ 10, 6, 10, 10, 10, 10 /),'t2m_lic'), & 73 ctrl_out((/ 10, 6, 10, 10, 10, 10 /),'t2m_oce'), & 74 ctrl_out((/ 10, 6, 10, 10, 10, 10 /),'t2m_sic') /) 75 75 76 76 type(ctrl_out),save :: o_wind10m = ctrl_out((/ 1, 1, 1, 10, 10, 10 /),'wind10m') … … 84 84 85 85 type(ctrl_out),save,dimension(4) :: o_u10m_srf = (/ ctrl_out((/ 10, 6, 10, 10, 10, 10 /),'u10m_ter'), & 86 87 88 86 ctrl_out((/ 10, 6, 10, 10, 10, 10 /),'u10m_lic'), & 87 ctrl_out((/ 10, 6, 10, 10, 10, 10 /),'u10m_oce'), & 88 ctrl_out((/ 10, 6, 10, 10, 10, 10 /),'u10m_sic') /) 89 89 90 90 type(ctrl_out),save,dimension(4) :: o_v10m_srf = (/ ctrl_out((/ 10, 6, 10, 10, 10, 10 /),'v10m_ter'), & 91 92 93 91 ctrl_out((/ 10, 6, 10, 10, 10, 10 /),'v10m_lic'), & 92 ctrl_out((/ 10, 6, 10, 10, 10, 10 /),'v10m_oce'), & 93 ctrl_out((/ 10, 6, 10, 10, 10, 10 /),'v10m_sic') /) 94 94 95 95 type(ctrl_out),save :: o_qsol = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'qsol') … … 103 103 type(ctrl_out),save :: o_evap = ctrl_out((/ 1, 1, 10, 10, 10, 10 /),'evap') 104 104 type(ctrl_out),save,dimension(4) :: o_evap_srf = (/ ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'evap_ter'), & 105 106 107 105 ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'evap_lic'), & 106 ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'evap_oce'), & 107 ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'evap_sic') /) 108 108 type(ctrl_out),save :: o_msnow = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'msnow') 109 109 type(ctrl_out),save :: o_fsnow = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'fsnow') … … 124 124 type(ctrl_out),save :: o_SWdn200clr = ctrl_out((/ 10, 1, 10, 10, 10, 10 /),'SWdn200clr') 125 125 126 ! arajouter127 ! type(ctrl_out),save :: o_LWupTOA = ctrl_out((/ 1, 4, 10, 10, 10, 10 /),'LWupTOA')128 ! type(ctrl_out),save :: o_LWupTOAclr = ctrl_out((/ 1, 4, 10, 10, 10, 10 /),'LWupTOAclr')129 ! type(ctrl_out),save :: o_LWdnTOA = ctrl_out((/ 1, 4, 10, 10, 10, 10 /),'LWdnTOA')130 ! type(ctrl_out),save :: o_LWdnTOAclr = ctrl_out((/ 1, 4, 10, 10, 10, 10 /),'LWdnTOAclr')126 ! arajouter 127 ! type(ctrl_out),save :: o_LWupTOA = ctrl_out((/ 1, 4, 10, 10, 10, 10 /),'LWupTOA') 128 ! type(ctrl_out),save :: o_LWupTOAclr = ctrl_out((/ 1, 4, 10, 10, 10, 10 /),'LWupTOAclr') 129 ! type(ctrl_out),save :: o_LWdnTOA = ctrl_out((/ 1, 4, 10, 10, 10, 10 /),'LWdnTOA') 130 ! type(ctrl_out),save :: o_LWdnTOAclr = ctrl_out((/ 1, 4, 10, 10, 10, 10 /),'LWdnTOAclr') 131 131 132 132 type(ctrl_out),save :: o_LWup200 = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'LWup200') … … 157 157 type(ctrl_out),save :: o_tauy = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'tauy') 158 158 type(ctrl_out),save,dimension(4) :: o_taux_srf = (/ ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'taux_ter'), & 159 160 161 159 ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'taux_lic'), & 160 ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'taux_oce'), & 161 ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'taux_sic') /) 162 162 163 163 type(ctrl_out),save,dimension(4) :: o_tauy_srf = (/ ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'tauy_ter'), & 164 165 166 164 ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'tauy_lic'), & 165 ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'tauy_oce'), & 166 ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'tauy_sic') /) 167 167 168 168 169 169 type(ctrl_out),save,dimension(4) :: o_pourc_srf = (/ ctrl_out((/ 1, 7, 10, 10, 10, 10 /),'pourc_ter'), & 170 171 172 170 ctrl_out((/ 1, 7, 10, 10, 10, 10 /),'pourc_lic'), & 171 ctrl_out((/ 1, 7, 10, 10, 10, 10 /),'pourc_oce'), & 172 ctrl_out((/ 1, 7, 10, 10, 10, 10 /),'pourc_sic') /) 173 173 174 174 type(ctrl_out),save,dimension(4) :: o_fract_srf = (/ ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'fract_ter'), & 175 176 177 175 ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'fract_lic'), & 176 ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'fract_oce'), & 177 ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'fract_sic') /) 178 178 179 179 type(ctrl_out),save,dimension(4) :: o_tsol_srf = (/ ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'tsol_ter'), & 180 181 182 180 ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'tsol_lic'), & 181 ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'tsol_oce'), & 182 ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'tsol_sic') /) 183 183 184 184 type(ctrl_out),save,dimension(4) :: o_sens_srf = (/ ctrl_out((/ 1, 6, 10, 7, 10, 10 /),'sens_ter'), & 185 186 187 185 ctrl_out((/ 1, 6, 10, 7, 10, 10 /),'sens_lic'), & 186 ctrl_out((/ 1, 6, 10, 7, 10, 10 /),'sens_oce'), & 187 ctrl_out((/ 1, 6, 10, 7, 10, 10 /),'sens_sic') /) 188 188 189 189 type(ctrl_out),save,dimension(4) :: o_lat_srf = (/ ctrl_out((/ 1, 6, 10, 7, 10, 10 /),'lat_ter'), & 190 191 192 190 ctrl_out((/ 1, 6, 10, 7, 10, 10 /),'lat_lic'), & 191 ctrl_out((/ 1, 6, 10, 7, 10, 10 /),'lat_oce'), & 192 ctrl_out((/ 1, 6, 10, 7, 10, 10 /),'lat_sic') /) 193 193 194 194 type(ctrl_out),save,dimension(4) :: o_flw_srf = (/ ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'flw_ter'), & 195 196 197 198 195 ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'flw_lic'), & 196 ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'flw_oce'), & 197 ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'flw_sic') /) 198 199 199 type(ctrl_out),save,dimension(4) :: o_fsw_srf = (/ ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'fsw_ter'), & 200 201 202 200 ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'fsw_lic'), & 201 ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'fsw_oce'), & 202 ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'fsw_sic') /) 203 203 204 204 type(ctrl_out),save,dimension(4) :: o_wbils_srf = (/ ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'wbils_ter'), & 205 206 207 205 ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'wbils_lic'), & 206 ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'wbils_oce'), & 207 ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'wbils_sic') /) 208 208 209 209 type(ctrl_out),save,dimension(4) :: o_wbilo_srf = (/ ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'wbilo_ter'), & 210 211 212 210 ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'wbilo_lic'), & 211 ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'wbilo_oce'), & 212 ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'wbilo_sic') /) 213 213 214 214 … … 226 226 type(ctrl_out),save :: o_uq = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'uq') 227 227 type(ctrl_out),save :: o_vq = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'vq') 228 228 229 229 type(ctrl_out),save :: o_cape = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'cape') 230 230 type(ctrl_out),save :: o_pbase = ctrl_out((/ 1, 5, 10, 10, 10, 10 /),'pbase') … … 240 240 type(ctrl_out),save :: o_s_lcl = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'s_lcl') 241 241 type(ctrl_out),save :: o_s_therm = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'s_therm') 242 !IM : Les champs suivants (s_capCL, s_oliqCL, s_cteiCL, s_trmb1, s_trmb2, s_trmb3) ne sont pas definis dans HBTM.F243 ! type(ctrl_out),save :: o_s_capCL = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'s_capCL')244 ! type(ctrl_out),save :: o_s_oliqCL = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'s_oliqCL')245 ! type(ctrl_out),save :: o_s_cteiCL = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'s_cteiCL')246 ! type(ctrl_out),save :: o_s_trmb1 = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'s_trmb1')247 ! type(ctrl_out),save :: o_s_trmb2 = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'s_trmb2')248 ! type(ctrl_out),save :: o_s_trmb3 = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'s_trmb3')242 !IM : Les champs suivants (s_capCL, s_oliqCL, s_cteiCL, s_trmb1, s_trmb2, s_trmb3) ne sont pas definis dans HBTM.F 243 ! type(ctrl_out),save :: o_s_capCL = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'s_capCL') 244 ! type(ctrl_out),save :: o_s_oliqCL = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'s_oliqCL') 245 ! type(ctrl_out),save :: o_s_cteiCL = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'s_cteiCL') 246 ! type(ctrl_out),save :: o_s_trmb1 = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'s_trmb1') 247 ! type(ctrl_out),save :: o_s_trmb2 = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'s_trmb2') 248 ! type(ctrl_out),save :: o_s_trmb3 = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'s_trmb3') 249 249 250 250 type(ctrl_out),save :: o_slab_bils = ctrl_out((/ 1, 1, 10, 10, 10, 10 /),'slab_bils_oce') … … 261 261 262 262 263 ! Champs interpolles sur des niveaux de pression ??? a faire correctement264 263 ! Champs interpolles sur des niveaux de pression ??? a faire correctement 264 265 265 type(ctrl_out),save,dimension(7) :: o_uSTDlevs = (/ ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'u850'), & 266 267 268 269 270 271 272 266 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'u700'), & 267 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'u500'), & 268 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'u200'), & 269 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'u100'), & 270 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'u50'), & 271 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'u10') /) 272 273 273 274 274 type(ctrl_out),save,dimension(7) :: o_vSTDlevs = (/ ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'v850'), & 275 276 277 278 279 280 275 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'v700'), & 276 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'v500'), & 277 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'v200'), & 278 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'v100'), & 279 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'v50'), & 280 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'v10') /) 281 281 282 282 type(ctrl_out),save,dimension(7) :: o_wSTDlevs = (/ ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'w850'), & 283 284 285 286 287 288 283 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'w700'), & 284 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'w500'), & 285 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'w200'), & 286 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'w100'), & 287 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'w50'), & 288 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'w10') /) 289 289 290 290 type(ctrl_out),save,dimension(7) :: o_tSTDlevs = (/ ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'t850'), & 291 292 293 294 295 296 291 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'t700'), & 292 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'t500'), & 293 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'t200'), & 294 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'t100'), & 295 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'t50'), & 296 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'t10') /) 297 297 298 298 type(ctrl_out),save,dimension(7) :: o_qSTDlevs = (/ ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'q850'), & 299 300 301 302 303 304 299 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'q700'), & 300 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'q500'), & 301 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'q200'), & 302 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'q100'), & 303 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'q50'), & 304 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'q10') /) 305 305 306 306 type(ctrl_out),save,dimension(7) :: o_zSTDlevs = (/ ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'z850'), & 307 308 309 310 311 312 307 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'z700'), & 308 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'z500'), & 309 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'z200'), & 310 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'z100'), & 311 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'z50'), & 312 ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'z10') /) 313 313 314 314 … … 318 318 type(ctrl_out),save :: o_dthmin = ctrl_out((/ 10, 1, 10, 10, 10, 10 /),'dthmin') 319 319 type(ctrl_out),save,dimension(4) :: o_u10_srf = (/ ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'u10_ter'), & 320 321 322 320 ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'u10_lic'), & 321 ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'u10_oce'), & 322 ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'u10_sic') /) 323 323 324 324 type(ctrl_out),save,dimension(4) :: o_v10_srf = (/ ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'v10_ter'), & 325 326 327 325 ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'v10_lic'), & 326 ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'v10_oce'), & 327 ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'v10_sic') /) 328 328 329 329 type(ctrl_out),save :: o_cldtau = ctrl_out((/ 10, 5, 10, 10, 10, 10 /),'cldtau') … … 339 339 340 340 type(ctrl_out),save,dimension(4) :: o_tke_srf = (/ ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'tke_ter'), & 341 342 343 341 ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'tke_lic'), & 342 ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'tke_oce'), & 343 ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'tke_sic') /) 344 344 345 345 type(ctrl_out),save,dimension(4) :: o_tke_max_srf = (/ ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'tke_max_ter'), & 346 347 348 346 ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'tke_max_lic'), & 347 ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'tke_max_oce'), & 348 ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'tke_max_sic') /) 349 349 350 350 type(ctrl_out),save :: o_kz = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'kz') … … 370 370 371 371 type(ctrl_out),save,dimension(10) :: o_tausumaero = (/ ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'OD550_ASBCM'), & 372 373 374 375 376 377 378 379 380 372 ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'OD550_ASPOMM'), & 373 ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'OD550_ASSO4M'), & 374 ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'OD550_CSSO4M'), & 375 ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'OD550_SSSSM'), & 376 ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'OD550_ASSSM'), & 377 ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'OD550_CSSSM'), & 378 ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'OD550_CIDUSTM'), & 379 ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'OD550_AIBCM'), & 380 ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'OD550_AIPOMM') /) 381 381 382 382 type(ctrl_out),save :: o_od550aer = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'od550aer') … … 465 465 466 466 type(ctrl_out),save,dimension(4) :: o_albe_srf = (/ ctrl_out((/ 3, 7, 10, 7, 10, 10 /),'albe_ter'), & 467 468 469 467 ctrl_out((/ 3, 7, 10, 7, 10, 10 /),'albe_lic'), & 468 ctrl_out((/ 3, 7, 10, 7, 10, 10 /),'albe_oce'), & 469 ctrl_out((/ 3, 7, 10, 7, 10, 10 /),'albe_sic') /) 470 470 471 471 type(ctrl_out),save,dimension(4) :: o_ages_srf = (/ ctrl_out((/ 10, 10, 10, 10, 10, 10 /),'ages_ter'), & 472 473 474 472 ctrl_out((/ 3, 10, 10, 10, 10, 10 /),'ages_lic'), & 473 ctrl_out((/ 10, 10, 10, 10, 10, 10 /),'ages_oce'), & 474 ctrl_out((/ 3, 10, 10, 10, 10, 10 /),'ages_sic') /) 475 475 476 476 type(ctrl_out),save,dimension(4) :: o_rugs_srf = (/ ctrl_out((/ 3, 6, 10, 10, 10, 10 /),'rugs_ter'), & 477 478 479 477 ctrl_out((/ 3, 6, 10, 10, 10, 10 /),'rugs_lic'), & 478 ctrl_out((/ 3, 6, 10, 10, 10, 10 /),'rugs_oce'), & 479 ctrl_out((/ 3, 6, 10, 10, 10, 10 /),'rugs_sic') /) 480 480 481 481 type(ctrl_out),save :: o_alb1 = ctrl_out((/ 3, 10, 10, 10, 10, 10 /),'alb1') … … 602 602 603 603 604 604 CONTAINS 605 605 606 606 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 607 607 !!!!!!!!! Ouverture des fichier et definition des variable de sortie !!!!!!!! 608 !! histbeg, histvert et histdef608 !! histbeg, histvert et histdef 609 609 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 610 610 611 611 SUBROUTINE phys_output_open(rlon,rlat,pim,tabij,ipt,jpt,plon,plat, & 612 612 jjmp1,nlevSTD,clevSTD,nbteta, & … … 617 617 new_aod, aerosol_couple) 618 618 619 USE iophy620 USE dimphy621 USE infotrac622 USE ioipsl623 USE mod_phys_lmdz_para624 USE aero_mod, only : naero_spc,name_aero625 626 IMPLICIT NONE627 include "dimensions.h"628 include "temps.h"629 include "indicesol.h"630 include "clesphys.h"631 include "thermcell.h"632 include "comvert.h"619 USE iophy 620 USE dimphy 621 USE infotrac 622 USE ioipsl 623 USE mod_phys_lmdz_para 624 USE aero_mod, only : naero_spc,name_aero 625 626 IMPLICIT NONE 627 include "dimensions.h" 628 include "temps.h" 629 include "indicesol.h" 630 include "clesphys.h" 631 include "thermcell.h" 632 include "comvert.h" 633 633 634 634 real,dimension(klon),intent(in) :: rlon … … 640 640 REAL,dimension(pim,2) :: plat_bounds, plon_bounds 641 641 642 integer :: jjmp1643 integer :: nbteta, nlevSTD, radpas644 logical :: ok_mensuel, ok_journe, ok_hf, ok_instan645 logical :: ok_LES,ok_ade,ok_aie646 logical :: new_aod, aerosol_couple647 integer, intent(in):: read_climoz ! read ozone climatology648 ! Allowed values are 0, 1 and 2649 ! 0: do not read an ozone climatology650 ! 1: read a single ozone climatology that will be used day and night651 ! 2: read two ozone climatologies, the average day and night652 ! climatology and the daylight climatology653 654 real :: dtime655 integer :: idayref656 real :: zjulian657 real, dimension(klev) :: Ahyb, Bhyb, Alt658 character(len=4), dimension(nlevSTD) :: clevSTD659 integer :: nsrf, k, iq, iiq, iff, i, j, ilev660 integer :: naero661 logical :: ok_veget662 integer :: iflag_pbl663 CHARACTER(len=4) :: bb2664 CHARACTER(len=2) :: bb3665 character(len=6) :: type_ocean666 CHARACTER(len=3) :: ctetaSTD(nbteta)667 real, dimension(nfiles) :: ecrit_files668 CHARACTER(len=20), dimension(nfiles) :: phys_out_filenames669 INTEGER, dimension(iim*jjmp1) :: ndex2d670 INTEGER, dimension(iim*jjmp1*klev) :: ndex3d671 integer :: imin_ins, imax_ins672 integer :: jmin_ins, jmax_ins673 integer, dimension(nfiles) :: phys_out_levmin, phys_out_levmax674 integer, dimension(nfiles) :: phys_out_filelevels675 CHARACTER(len=20), dimension(nfiles) :: type_ecri_files, phys_out_filetypes676 character(len=20), dimension(nfiles) :: chtimestep = (/ 'DefFreq', 'DefFreq','DefFreq', 'DefFreq', 'DefFreq', 'DefFreq' /)677 logical, dimension(nfiles) :: phys_out_filekeys678 logical, dimension(nfiles) :: phys_out_filestations642 integer :: jjmp1 643 integer :: nbteta, nlevSTD, radpas 644 logical :: ok_mensuel, ok_journe, ok_hf, ok_instan 645 logical :: ok_LES,ok_ade,ok_aie 646 logical :: new_aod, aerosol_couple 647 integer, intent(in):: read_climoz ! read ozone climatology 648 ! Allowed values are 0, 1 and 2 649 ! 0: do not read an ozone climatology 650 ! 1: read a single ozone climatology that will be used day and night 651 ! 2: read two ozone climatologies, the average day and night 652 ! climatology and the daylight climatology 653 654 real :: dtime 655 integer :: idayref 656 real :: zjulian 657 real, dimension(klev) :: Ahyb, Bhyb, Alt 658 character(len=4), dimension(nlevSTD) :: clevSTD 659 integer :: nsrf, k, iq, iiq, iff, i, j, ilev 660 integer :: naero 661 logical :: ok_veget 662 integer :: iflag_pbl 663 CHARACTER(len=4) :: bb2 664 CHARACTER(len=2) :: bb3 665 character(len=6) :: type_ocean 666 CHARACTER(len=3) :: ctetaSTD(nbteta) 667 real, dimension(nfiles) :: ecrit_files 668 CHARACTER(len=20), dimension(nfiles) :: phys_out_filenames 669 INTEGER, dimension(iim*jjmp1) :: ndex2d 670 INTEGER, dimension(iim*jjmp1*klev) :: ndex3d 671 integer :: imin_ins, imax_ins 672 integer :: jmin_ins, jmax_ins 673 integer, dimension(nfiles) :: phys_out_levmin, phys_out_levmax 674 integer, dimension(nfiles) :: phys_out_filelevels 675 CHARACTER(len=20), dimension(nfiles) :: type_ecri_files, phys_out_filetypes 676 character(len=20), dimension(nfiles) :: chtimestep = (/ 'DefFreq', 'DefFreq','DefFreq', 'DefFreq', 'DefFreq', 'DefFreq' /) 677 logical, dimension(nfiles) :: phys_out_filekeys 678 logical, dimension(nfiles) :: phys_out_filestations 679 679 680 680 !!!!!!!!!! stockage dans une region limitee pour chaque fichier !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 681 ! entre [phys_out_lonmin,phys_out_lonmax] et [phys_out_latmin,phys_out_latmax]682 683 logical, dimension(nfiles), save :: phys_out_regfkey = (/ .false., .false., .false., .false., .false., .false. /)684 real, dimension(nfiles), save :: phys_out_lonmin = (/ -180., -180., -180., -180., -180., -180. /)685 real, dimension(nfiles), save :: phys_out_lonmax = (/ 180., 180., 180., 180., 180., 180. /)686 real, dimension(nfiles), save :: phys_out_latmin = (/ -90., -90., -90., -90., -90., -90. /)687 real, dimension(nfiles), save :: phys_out_latmax = (/ 90., 90., 90., 90., 90., 90. /)688 689 print*,'Debut phys_output_mod.F90'690 ! Initialisations (Valeurs par defaut691 692 if (.not. allocated(o_trac)) ALLOCATE(o_trac(nqtot))693 694 levmax = (/ klev, klev, klev, klev, klev, klev /)695 696 phys_out_filenames(1) = 'histmth'697 phys_out_filenames(2) = 'histday'698 phys_out_filenames(3) = 'histhf'699 phys_out_filenames(4) = 'histins'700 phys_out_filenames(5) = 'histLES'701 phys_out_filenames(6) = 'histstn'702 703 type_ecri(1) = 'ave(X)'704 type_ecri(2) = 'ave(X)'705 type_ecri(3) = 'ave(X)'706 type_ecri(4) = 'inst(X)'707 type_ecri(5) = 'ave(X)'708 type_ecri(6) = 'inst(X)'709 710 clef_files(1) = ok_mensuel711 clef_files(2) = ok_journe712 clef_files(3) = ok_hf713 clef_files(4) = ok_instan714 clef_files(5) = ok_LES715 clef_files(6) = ok_instan716 717 !sortir des fichiers "stations" si clef_stations(:)=.TRUE.718 clef_stations(1) = .FALSE.719 clef_stations(2) = .FALSE.720 clef_stations(3) = .FALSE.721 clef_stations(4) = .FALSE.722 clef_stations(5) = .FALSE.723 clef_stations(6) = .FALSE.724 725 lev_files(1) = lev_histmth726 lev_files(2) = lev_histday727 lev_files(3) = lev_histhf728 lev_files(4) = lev_histins729 lev_files(5) = lev_histLES730 lev_files(6) = lev_histins731 732 ecrit_files(1) = ecrit_mth733 ecrit_files(2) = ecrit_day734 ecrit_files(3) = ecrit_hf735 ecrit_files(4) = ecrit_ins736 ecrit_files(5) = ecrit_LES737 ecrit_files(6) = ecrit_ins738 739 !! Lectures des parametres de sorties dans physiq.def740 741 call getin('phys_out_regfkey',phys_out_regfkey)742 call getin('phys_out_lonmin',phys_out_lonmin)743 call getin('phys_out_lonmax',phys_out_lonmax)744 call getin('phys_out_latmin',phys_out_latmin)745 call getin('phys_out_latmax',phys_out_latmax)746 747 call getin('phys_out_levmin',levmin)748 749 call getin('phys_out_levmax',levmax)750 call getin('phys_out_filenames',phys_out_filenames)751 752 call getin('phys_out_filekeys',clef_files)753 754 call getin('phys_out_filestations',clef_stations)755 756 call getin('phys_out_filelevels',lev_files)757 call getin('phys_out_filetimesteps',chtimestep)758 759 call getin('phys_out_filetypes',type_ecri)760 761 type_ecri_files(:)=type_ecri(:)762 763 print*,'phys_out_lonmin=',phys_out_lonmin764 print*,'phys_out_lonmax=',phys_out_lonmax765 print*,'phys_out_latmin=',phys_out_latmin766 print*,'phys_out_latmax=',phys_out_latmax767 print*,'phys_out_filenames=',phys_out_filenames768 print*,'phys_out_filetypes=',type_ecri769 print*,'phys_out_filekeys=',clef_files770 print*,'phys_out_filestations=',clef_stations771 print*,'phys_out_filelevels=',lev_files681 ! entre [phys_out_lonmin,phys_out_lonmax] et [phys_out_latmin,phys_out_latmax] 682 683 logical, dimension(nfiles), save :: phys_out_regfkey = (/ .false., .false., .false., .false., .false., .false. /) 684 real, dimension(nfiles), save :: phys_out_lonmin = (/ -180., -180., -180., -180., -180., -180. /) 685 real, dimension(nfiles), save :: phys_out_lonmax = (/ 180., 180., 180., 180., 180., 180. /) 686 real, dimension(nfiles), save :: phys_out_latmin = (/ -90., -90., -90., -90., -90., -90. /) 687 real, dimension(nfiles), save :: phys_out_latmax = (/ 90., 90., 90., 90., 90., 90. /) 688 689 print*,'Debut phys_output_mod.F90' 690 ! Initialisations (Valeurs par defaut 691 692 if (.not. allocated(o_trac)) ALLOCATE(o_trac(nqtot)) 693 694 levmax = (/ klev, klev, klev, klev, klev, klev /) 695 696 phys_out_filenames(1) = 'histmth' 697 phys_out_filenames(2) = 'histday' 698 phys_out_filenames(3) = 'histhf' 699 phys_out_filenames(4) = 'histins' 700 phys_out_filenames(5) = 'histLES' 701 phys_out_filenames(6) = 'histstn' 702 703 type_ecri(1) = 'ave(X)' 704 type_ecri(2) = 'ave(X)' 705 type_ecri(3) = 'ave(X)' 706 type_ecri(4) = 'inst(X)' 707 type_ecri(5) = 'ave(X)' 708 type_ecri(6) = 'inst(X)' 709 710 clef_files(1) = ok_mensuel 711 clef_files(2) = ok_journe 712 clef_files(3) = ok_hf 713 clef_files(4) = ok_instan 714 clef_files(5) = ok_LES 715 clef_files(6) = ok_instan 716 717 !sortir des fichiers "stations" si clef_stations(:)=.TRUE. 718 clef_stations(1) = .FALSE. 719 clef_stations(2) = .FALSE. 720 clef_stations(3) = .FALSE. 721 clef_stations(4) = .FALSE. 722 clef_stations(5) = .FALSE. 723 clef_stations(6) = .FALSE. 724 725 lev_files(1) = lev_histmth 726 lev_files(2) = lev_histday 727 lev_files(3) = lev_histhf 728 lev_files(4) = lev_histins 729 lev_files(5) = lev_histLES 730 lev_files(6) = lev_histins 731 732 ecrit_files(1) = ecrit_mth 733 ecrit_files(2) = ecrit_day 734 ecrit_files(3) = ecrit_hf 735 ecrit_files(4) = ecrit_ins 736 ecrit_files(5) = ecrit_LES 737 ecrit_files(6) = ecrit_ins 738 739 !! Lectures des parametres de sorties dans physiq.def 740 741 call getin('phys_out_regfkey',phys_out_regfkey) 742 call getin('phys_out_lonmin',phys_out_lonmin) 743 call getin('phys_out_lonmax',phys_out_lonmax) 744 call getin('phys_out_latmin',phys_out_latmin) 745 call getin('phys_out_latmax',phys_out_latmax) 746 phys_out_levmin(:)=levmin(:) 747 call getin('phys_out_levmin',levmin) 748 phys_out_levmax(:)=levmax(:) 749 call getin('phys_out_levmax',levmax) 750 call getin('phys_out_filenames',phys_out_filenames) 751 phys_out_filekeys(:)=clef_files(:) 752 call getin('phys_out_filekeys',clef_files) 753 phys_out_filestations(:)=clef_stations(:) 754 call getin('phys_out_filestations',clef_stations) 755 phys_out_filelevels(:)=lev_files(:) 756 call getin('phys_out_filelevels',lev_files) 757 call getin('phys_out_filetimesteps',chtimestep) 758 phys_out_filetypes(:)=type_ecri(:) 759 call getin('phys_out_filetypes',type_ecri) 760 761 type_ecri_files(:)=type_ecri(:) 762 763 print*,'phys_out_lonmin=',phys_out_lonmin 764 print*,'phys_out_lonmax=',phys_out_lonmax 765 print*,'phys_out_latmin=',phys_out_latmin 766 print*,'phys_out_latmax=',phys_out_latmax 767 print*,'phys_out_filenames=',phys_out_filenames 768 print*,'phys_out_filetypes=',type_ecri 769 print*,'phys_out_filekeys=',clef_files 770 print*,'phys_out_filestations=',clef_stations 771 print*,'phys_out_filelevels=',lev_files 772 772 773 773 !!!!!!!!!!!!!!!!!!!!!!! Boucle sur les fichiers !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 774 ! Appel de histbeg et histvert pour creer le fichier et les niveaux verticaux !!775 ! Appel des histbeg pour definir les variables (nom, moy ou inst, freq de sortie ..774 ! Appel de histbeg et histvert pour creer le fichier et les niveaux verticaux !! 775 ! Appel des histbeg pour definir les variables (nom, moy ou inst, freq de sortie .. 776 776 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 777 777 778 zdtime = dtime ! Frequence ou l on moyenne779 780 ! Calcul des Ahyb, Bhyb et Alt781 782 783 784 785 786 ! if(prt_level.ge.1) then787 788 789 790 ! endif791 DO iff=1,nfiles792 793 IF (clef_files(iff)) THEN794 795 if ( chtimestep(iff).eq.'DefFreq' ) then796 ! Par defaut ecrit_files = (ecrit_mensuel ecrit_jour ecrit_hf ...)*86400.797 ecrit_files(iff)=ecrit_files(iff)*86400.798 else799 call convers_timesteps(chtimestep(iff),dtime,ecrit_files(iff))800 endif801 print*,'ecrit_files(',iff,')= ',ecrit_files(iff)802 803 zoutm(iff) = ecrit_files(iff) ! Frequence ou l on ecrit en seconde804 805 idayref = day_ref806 CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)778 zdtime = dtime ! Frequence ou l on moyenne 779 780 ! Calcul des Ahyb, Bhyb et Alt 781 do k=1,klev 782 Ahyb(k)=(ap(k)+ap(k+1))/2. 783 Bhyb(k)=(bp(k)+bp(k+1))/2. 784 Alt(k)=log(preff/presnivs(k))*8. 785 enddo 786 ! if(prt_level.ge.1) then 787 print*,'Ap Hybrid = ',Ahyb(1:klev) 788 print*,'Bp Hybrid = ',Bhyb(1:klev) 789 print*,'Alt approx des couches pour une haut d echelle de 8km = ',Alt(1:klev) 790 ! endif 791 DO iff=1,nfiles 792 793 IF (clef_files(iff)) THEN 794 795 if ( chtimestep(iff).eq.'DefFreq' ) then 796 ! Par defaut ecrit_files = (ecrit_mensuel ecrit_jour ecrit_hf ...)*86400. 797 ecrit_files(iff)=ecrit_files(iff)*86400. 798 else 799 call convers_timesteps(chtimestep(iff),dtime,ecrit_files(iff)) 800 endif 801 print*,'ecrit_files(',iff,')= ',ecrit_files(iff) 802 803 zoutm(iff) = ecrit_files(iff) ! Frequence ou l on ecrit en seconde 804 805 idayref = day_ref 806 CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian) 807 807 808 808 !!!!!!!!!!!!!!!!! Traitement dans le cas ou l'on veut stocker sur un domaine limite !! 809 809 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 810 if (phys_out_regfkey(iff)) then811 812 imin_ins=1813 imax_ins=iim814 jmin_ins=1815 jmax_ins=jjmp1816 817 ! correction abderr818 do i=1,iim819 print*,'io_lon(i)=',io_lon(i)820 if (io_lon(i).le.phys_out_lonmin(iff)) imin_ins=i821 if (io_lon(i).le.phys_out_lonmax(iff)) imax_ins=i+1822 enddo823 824 do j=1,jjmp1825 print*,'io_lat(j)=',io_lat(j)826 if (io_lat(j).ge.phys_out_latmin(iff)) jmax_ins=j+1827 if (io_lat(j).ge.phys_out_latmax(iff)) jmin_ins=j828 enddo829 830 print*,'On stoke le fichier histoire numero ',iff,' sur ', &831 imin_ins,imax_ins,jmin_ins,jmax_ins832 print*,'longitudes : ', &833 io_lon(imin_ins),io_lon(imax_ins), &834 'latitudes : ', &835 io_lat(jmax_ins),io_lat(jmin_ins)836 837 CALL histbeg(phys_out_filenames(iff),iim,io_lon,jjmp1,io_lat, &838 imin_ins,imax_ins-imin_ins+1, &839 jmin_ins,jmax_ins-jmin_ins+1, &840 itau_phy,zjulian,dtime,nhorim(iff),nid_files(iff))810 if (phys_out_regfkey(iff)) then 811 812 imin_ins=1 813 imax_ins=iim 814 jmin_ins=1 815 jmax_ins=jjmp1 816 817 ! correction abderr 818 do i=1,iim 819 print*,'io_lon(i)=',io_lon(i) 820 if (io_lon(i).le.phys_out_lonmin(iff)) imin_ins=i 821 if (io_lon(i).le.phys_out_lonmax(iff)) imax_ins=i+1 822 enddo 823 824 do j=1,jjmp1 825 print*,'io_lat(j)=',io_lat(j) 826 if (io_lat(j).ge.phys_out_latmin(iff)) jmax_ins=j+1 827 if (io_lat(j).ge.phys_out_latmax(iff)) jmin_ins=j 828 enddo 829 830 print*,'On stoke le fichier histoire numero ',iff,' sur ', & 831 imin_ins,imax_ins,jmin_ins,jmax_ins 832 print*,'longitudes : ', & 833 io_lon(imin_ins),io_lon(imax_ins), & 834 'latitudes : ', & 835 io_lat(jmax_ins),io_lat(jmin_ins) 836 837 CALL histbeg(phys_out_filenames(iff),iim,io_lon,jjmp1,io_lat, & 838 imin_ins,imax_ins-imin_ins+1, & 839 jmin_ins,jmax_ins-jmin_ins+1, & 840 itau_phy,zjulian,dtime,nhorim(iff),nid_files(iff)) 841 841 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 842 !IM fichiers stations843 else if (clef_stations(iff)) THEN844 845 print*,'phys_output_mod phys_out_filenames=',phys_out_filenames(iff)846 847 call histbeg_phy_all(rlon,rlat,pim,tabij,ipt,jpt,plon,plat,plon_bounds,plat_bounds, &848 849 850 else851 CALL histbeg_phy(phys_out_filenames(iff),itau_phy,zjulian,dtime,nhorim(iff),nid_files(iff))852 endif853 854 CALL histvert(nid_files(iff), "presnivs", "Vertical levels", "Pa", &855 levmax(iff) - levmin(iff) + 1, &856 presnivs(levmin(iff):levmax(iff)), nvertm(iff),"down")842 !IM fichiers stations 843 else if (clef_stations(iff)) THEN 844 845 print*,'phys_output_mod phys_out_filenames=',phys_out_filenames(iff) 846 847 call histbeg_phy_all(rlon,rlat,pim,tabij,ipt,jpt,plon,plat,plon_bounds,plat_bounds, & 848 phys_out_filenames(iff), & 849 itau_phy,zjulian,dtime,nhorim(iff),nid_files(iff)) 850 else 851 CALL histbeg_phy(phys_out_filenames(iff),itau_phy,zjulian,dtime,nhorim(iff),nid_files(iff)) 852 endif 853 854 CALL histvert(nid_files(iff), "presnivs", "Vertical levels", "Pa", & 855 levmax(iff) - levmin(iff) + 1, & 856 presnivs(levmin(iff):levmax(iff)), nvertm(iff),"down") 857 857 858 858 !!!!!!!!!!!!! Traitement des champs 3D pour histhf !!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 859 859 !!!!!!!!!!!!!!! A Revoir plus tard !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 860 ! IF (iff.eq.3.and.lev_files(iff).ge.4) THEN861 ! CALL histbeg_phy("histhf3d",itau_phy, &862 ! & zjulian, dtime, &863 ! & nhorim, nid_hf3d)864 865 ! CALL histvert(nid_hf3d, "presnivs", &866 ! & "Vertical levels", "mb", &867 ! & klev, presnivs/100., nvertm)868 ! ENDIF869 !860 ! IF (iff.eq.3.and.lev_files(iff).ge.4) THEN 861 ! CALL histbeg_phy("histhf3d",itau_phy, & 862 ! & zjulian, dtime, & 863 ! & nhorim, nid_hf3d) 864 865 ! CALL histvert(nid_hf3d, "presnivs", & 866 ! & "Vertical levels", "mb", & 867 ! & klev, presnivs/100., nvertm) 868 ! ENDIF 869 ! 870 870 !!!! Composentes de la coordonnee sigma-hybride 871 CALL histvert(nid_files(iff), "Ahyb","Ahyb comp of Hyb Cord ", "Pa", &872 873 874 CALL histvert(nid_files(iff), "Bhyb","Bhyb comp of Hyb Cord", " ", &875 876 877 CALL histvert(nid_files(iff), "Alt","Height approx for scale heigh of 8km at levels", "Km", &878 879 880 ! CALL histvert(nid_files(iff), "preff","Reference pressure", "Pa", &881 ! 1,preff,nvertp0(iff))871 CALL histvert(nid_files(iff), "Ahyb","Ahyb comp of Hyb Cord ", "Pa", & 872 levmax(iff) - levmin(iff) + 1,Ahyb,nvertap(iff)) 873 874 CALL histvert(nid_files(iff), "Bhyb","Bhyb comp of Hyb Cord", " ", & 875 levmax(iff) - levmin(iff) + 1,Bhyb,nvertbp(iff)) 876 877 CALL histvert(nid_files(iff), "Alt","Height approx for scale heigh of 8km at levels", "Km", & 878 levmax(iff) - levmin(iff) + 1,Alt,nvertAlt(iff)) 879 880 ! CALL histvert(nid_files(iff), "preff","Reference pressure", "Pa", & 881 ! 1,preff,nvertp0(iff)) 882 882 !!! Champs 1D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 883 IF (.NOT.clef_stations(iff)) THEN884 !885 !IM: there is no way to have one single value in a netcdf file886 !887 type_ecri(1) = 'once'888 type_ecri(2) = 'once'889 type_ecri(3) = 'once'890 type_ecri(4) = 'once'891 type_ecri(5) = 'once'892 type_ecri(6) = 'once'893 CALL histdef2d(iff,clef_stations(iff),o_aire%flag,o_aire%name,"Grid area", "-")894 CALL histdef2d(iff,clef_stations(iff),o_contfracATM%flag,o_contfracATM%name,"% sfce ter+lic", "-")895 ENDIF896 type_ecri(:) = type_ecri_files(:)883 IF (.NOT.clef_stations(iff)) THEN 884 ! 885 !IM: there is no way to have one single value in a netcdf file 886 ! 887 type_ecri(1) = 'once' 888 type_ecri(2) = 'once' 889 type_ecri(3) = 'once' 890 type_ecri(4) = 'once' 891 type_ecri(5) = 'once' 892 type_ecri(6) = 'once' 893 CALL histdef2d(iff,clef_stations(iff),o_aire%flag,o_aire%name,"Grid area", "-") 894 CALL histdef2d(iff,clef_stations(iff),o_contfracATM%flag,o_contfracATM%name,"% sfce ter+lic", "-") 895 ENDIF 896 type_ecri(:) = type_ecri_files(:) 897 897 898 898 !!! Champs 2D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 899 CALL histdef2d(iff,clef_stations(iff),o_phis%flag,o_phis%name,"Surface geop.height", "m2/s2" ) 900 CALL histdef2d(iff,clef_stations(iff),o_contfracOR%flag,o_contfracOR%name,"% sfce terre OR", "-" ) 901 CALL histdef2d(iff,clef_stations(iff),o_aireTER%flag,o_aireTER%name,"Grid area CONT", "-" ) 902 CALL histdef2d(iff,clef_stations(iff),o_flat%flag,o_flat%name, "Latent heat flux", "W/m2") 903 CALL histdef2d(iff,clef_stations(iff),o_slp%flag,o_slp%name, "Sea Level Pressure", "Pa" ) 904 CALL histdef2d(iff,clef_stations(iff),o_tsol%flag,o_tsol%name, "Surface Temperature", "K") 905 CALL histdef2d(iff,clef_stations(iff),o_t2m%flag,o_t2m%name, "Temperature 2m", "K" ) 906 IF (.NOT.clef_stations(iff)) THEN 907 ! 908 !IM: there is no way to have one single value in a netcdf file 909 ! 910 type_ecri(1) = 't_min(X)' 911 type_ecri(2) = 't_min(X)' 912 type_ecri(3) = 't_min(X)' 913 type_ecri(4) = 't_min(X)' 914 type_ecri(5) = 't_min(X)' 915 type_ecri(6) = 't_min(X)' 916 CALL histdef2d(iff,clef_stations(iff),o_t2m_min%flag,o_t2m_min%name, "Temp 2m min", "K" ) 917 type_ecri(1) = 't_max(X)' 918 type_ecri(2) = 't_max(X)' 919 type_ecri(3) = 't_max(X)' 920 type_ecri(4) = 't_max(X)' 921 type_ecri(5) = 't_max(X)' 922 type_ecri(6) = 't_max(X)' 923 CALL histdef2d(iff,clef_stations(iff),o_t2m_max%flag,o_t2m_max%name, "Temp 2m max", "K" ) 924 ENDIF 925 type_ecri(:) = type_ecri_files(:) 926 CALL histdef2d(iff,clef_stations(iff),o_wind10m%flag,o_wind10m%name, "10-m wind speed", "m/s") 927 CALL histdef2d(iff,clef_stations(iff),o_wind10max%flag,o_wind10max%name, "10m wind speed max", "m/s") 928 CALL histdef2d(iff,clef_stations(iff),o_sicf%flag,o_sicf%name, "Sea-ice fraction", "-" ) 929 CALL histdef2d(iff,clef_stations(iff),o_q2m%flag,o_q2m%name, "Specific humidity 2m", "kg/kg") 930 CALL histdef2d(iff,clef_stations(iff),o_u10m%flag,o_u10m%name, "Vent zonal 10m", "m/s" ) 931 CALL histdef2d(iff,clef_stations(iff),o_v10m%flag,o_v10m%name, "Vent meridien 10m", "m/s") 932 CALL histdef2d(iff,clef_stations(iff),o_psol%flag,o_psol%name, "Surface Pressure", "Pa" ) 933 CALL histdef2d(iff,clef_stations(iff),o_qsurf%flag,o_qsurf%name, "Surface Air humidity", "kg/kg") 934 935 if (.not. ok_veget) then 936 CALL histdef2d(iff,clef_stations(iff),o_qsol%flag,o_qsol%name, "Soil watter content", "mm" ) 937 endif 938 939 CALL histdef2d(iff,clef_stations(iff),o_ndayrain%flag,o_ndayrain%name, "Number of dayrain(liq+sol)", "-") 940 CALL histdef2d(iff,clef_stations(iff),o_precip%flag,o_precip%name, "Precip Totale liq+sol", "kg/(s*m2)" ) 941 CALL histdef2d(iff,clef_stations(iff),o_plul%flag,o_plul%name, "Large-scale Precip.", "kg/(s*m2)") 942 CALL histdef2d(iff,clef_stations(iff),o_pluc%flag,o_pluc%name, "Convective Precip.", "kg/(s*m2)") 943 CALL histdef2d(iff,clef_stations(iff),o_snow%flag,o_snow%name, "Snow fall", "kg/(s*m2)" ) 944 CALL histdef2d(iff,clef_stations(iff),o_msnow%flag,o_msnow%name, "Surface snow amount", "kg/m2" ) 945 CALL histdef2d(iff,clef_stations(iff),o_fsnow%flag,o_fsnow%name, "Surface snow area fraction", "-" ) 946 CALL histdef2d(iff,clef_stations(iff),o_evap%flag,o_evap%name, "Evaporat", "kg/(s*m2)" ) 947 CALL histdef2d(iff,clef_stations(iff),o_tops%flag,o_tops%name, "Solar rad. at TOA", "W/m2") 948 CALL histdef2d(iff,clef_stations(iff),o_tops0%flag,o_tops0%name, "CS Solar rad. at TOA", "W/m2") 949 CALL histdef2d(iff,clef_stations(iff),o_topl%flag,o_topl%name, "IR rad. at TOA", "W/m2" ) 950 CALL histdef2d(iff,clef_stations(iff),o_topl0%flag,o_topl0%name, "IR rad. at TOA", "W/m2") 951 CALL histdef2d(iff,clef_stations(iff),o_SWupTOA%flag,o_SWupTOA%name, "SWup at TOA", "W/m2") 952 CALL histdef2d(iff,clef_stations(iff),o_SWupTOAclr%flag,o_SWupTOAclr%name, "SWup clear sky at TOA", "W/m2") 953 CALL histdef2d(iff,clef_stations(iff),o_SWdnTOA%flag,o_SWdnTOA%name, "SWdn at TOA", "W/m2" ) 954 CALL histdef2d(iff,clef_stations(iff),o_SWdnTOAclr%flag,o_SWdnTOAclr%name, "SWdn clear sky at TOA", "W/m2") 955 CALL histdef2d(iff,clef_stations(iff),o_nettop%flag,o_nettop%name, "Net dn radiatif flux at TOA", "W/m2") 956 CALL histdef2d(iff,clef_stations(iff),o_SWup200%flag,o_SWup200%name, "SWup at 200mb", "W/m2" ) 957 CALL histdef2d(iff,clef_stations(iff),o_SWup200clr%flag,o_SWup200clr%name, "SWup clear sky at 200mb", "W/m2") 958 CALL histdef2d(iff,clef_stations(iff),o_SWdn200%flag,o_SWdn200%name, "SWdn at 200mb", "W/m2" ) 959 CALL histdef2d(iff,clef_stations(iff),o_SWdn200clr%flag,o_SWdn200clr%name, "SWdn clear sky at 200mb", "W/m2") 960 CALL histdef2d(iff,clef_stations(iff),o_LWup200%flag,o_LWup200%name, "LWup at 200mb", "W/m2") 961 CALL histdef2d(iff,clef_stations(iff),o_LWup200clr%flag,o_LWup200clr%name, "LWup clear sky at 200mb", "W/m2") 962 CALL histdef2d(iff,clef_stations(iff),o_LWdn200%flag,o_LWdn200%name, "LWdn at 200mb", "W/m2") 963 CALL histdef2d(iff,clef_stations(iff),o_LWdn200clr%flag,o_LWdn200clr%name, "LWdn clear sky at 200mb", "W/m2") 964 CALL histdef2d(iff,clef_stations(iff),o_sols%flag,o_sols%name, "Solar rad. at surf.", "W/m2") 965 CALL histdef2d(iff,clef_stations(iff),o_sols0%flag,o_sols0%name, "Solar rad. at surf.", "W/m2") 966 CALL histdef2d(iff,clef_stations(iff),o_soll%flag,o_soll%name, "IR rad. at surface", "W/m2") 967 CALL histdef2d(iff,clef_stations(iff),o_radsol%flag,o_radsol%name, "Rayonnement au sol", "W/m2") 968 CALL histdef2d(iff,clef_stations(iff),o_soll0%flag,o_soll0%name, "IR rad. at surface", "W/m2") 969 CALL histdef2d(iff,clef_stations(iff),o_SWupSFC%flag,o_SWupSFC%name, "SWup at surface", "W/m2") 970 CALL histdef2d(iff,clef_stations(iff),o_SWupSFCclr%flag,o_SWupSFCclr%name, "SWup clear sky at surface", "W/m2") 971 CALL histdef2d(iff,clef_stations(iff),o_SWdnSFC%flag,o_SWdnSFC%name, "SWdn at surface", "W/m2") 972 CALL histdef2d(iff,clef_stations(iff),o_SWdnSFCclr%flag,o_SWdnSFCclr%name, "SWdn clear sky at surface", "W/m2") 973 CALL histdef2d(iff,clef_stations(iff),o_LWupSFC%flag,o_LWupSFC%name, "Upwd. IR rad. at surface", "W/m2") 974 CALL histdef2d(iff,clef_stations(iff),o_LWdnSFC%flag,o_LWdnSFC%name, "Down. IR rad. at surface", "W/m2") 975 CALL histdef2d(iff,clef_stations(iff),o_LWupSFCclr%flag,o_LWupSFCclr%name, "CS Upwd. IR rad. at surface", "W/m2") 976 CALL histdef2d(iff,clef_stations(iff),o_LWdnSFCclr%flag,o_LWdnSFCclr%name, "Down. CS IR rad. at surface", "W/m2") 977 CALL histdef2d(iff,clef_stations(iff),o_bils%flag,o_bils%name, "Surf. total heat flux", "W/m2") 978 CALL histdef2d(iff,clef_stations(iff),o_sens%flag,o_sens%name, "Sensible heat flux", "W/m2") 979 CALL histdef2d(iff,clef_stations(iff),o_fder%flag,o_fder%name, "Heat flux derivation", "W/m2") 980 CALL histdef2d(iff,clef_stations(iff),o_ffonte%flag,o_ffonte%name, "Thermal flux for snow melting", "W/m2") 981 CALL histdef2d(iff,clef_stations(iff),o_fqcalving%flag,o_fqcalving%name, "Ice Calving", "kg/m2/s") 982 CALL histdef2d(iff,clef_stations(iff),o_fqfonte%flag,o_fqfonte%name, "Land ice melt", "kg/m2/s") 983 984 CALL histdef2d(iff,clef_stations(iff),o_taux%flag,o_taux%name, "Zonal wind stress","Pa") 985 CALL histdef2d(iff,clef_stations(iff),o_tauy%flag,o_tauy%name, "Meridional wind stress","Pa") 986 987 DO nsrf = 1, nbsrf 988 CALL histdef2d(iff,clef_stations(iff),o_pourc_srf(nsrf)%flag,o_pourc_srf(nsrf)%name,"% "//clnsurf(nsrf),"%") 989 CALL histdef2d(iff,clef_stations(iff),o_fract_srf(nsrf)%flag,o_fract_srf(nsrf)%name,"Fraction "//clnsurf(nsrf),"1") 990 CALL histdef2d(iff,clef_stations(iff), & 991 o_taux_srf(nsrf)%flag,o_taux_srf(nsrf)%name,"Zonal wind stress"//clnsurf(nsrf),"Pa") 992 CALL histdef2d(iff,clef_stations(iff), & 993 o_tauy_srf(nsrf)%flag,o_tauy_srf(nsrf)%name,"Meridional wind stress "//clnsurf(nsrf),"Pa") 994 CALL histdef2d(iff,clef_stations(iff), & 995 o_tsol_srf(nsrf)%flag,o_tsol_srf(nsrf)%name,"Temperature "//clnsurf(nsrf),"K") 996 CALL histdef2d(iff,clef_stations(iff), & 997 o_u10m_srf(nsrf)%flag,o_u10m_srf(nsrf)%name,"Vent Zonal 10m "//clnsurf(nsrf),"m/s") 998 CALL histdef2d(iff,clef_stations(iff), & 999 o_evap_srf(nsrf)%flag,o_evap_srf(nsrf)%name,"evaporation at surface "//clnsurf(nsrf),"kg/(s*m2)") 1000 CALL histdef2d(iff,clef_stations(iff), & 1001 o_v10m_srf(nsrf)%flag,o_v10m_srf(nsrf)%name,"Vent meredien 10m "//clnsurf(nsrf),"m/s") 1002 CALL histdef2d(iff,clef_stations(iff), & 1003 o_t2m_srf(nsrf)%flag,o_t2m_srf(nsrf)%name,"Temp 2m "//clnsurf(nsrf),"K") 1004 CALL histdef2d(iff,clef_stations(iff), & 1005 o_sens_srf(nsrf)%flag,o_sens_srf(nsrf)%name,"Sensible heat flux "//clnsurf(nsrf),"W/m2") 1006 CALL histdef2d(iff,clef_stations(iff), & 1007 o_lat_srf(nsrf)%flag,o_lat_srf(nsrf)%name,"Latent heat flux "//clnsurf(nsrf),"W/m2") 1008 CALL histdef2d(iff,clef_stations(iff), & 1009 o_flw_srf(nsrf)%flag,o_flw_srf(nsrf)%name,"LW "//clnsurf(nsrf),"W/m2") 1010 CALL histdef2d(iff,clef_stations(iff), & 1011 o_fsw_srf(nsrf)%flag,o_fsw_srf(nsrf)%name,"SW "//clnsurf(nsrf),"W/m2") 1012 CALL histdef2d(iff,clef_stations(iff), & 1013 o_wbils_srf(nsrf)%flag,o_wbils_srf(nsrf)%name,"Bilan sol "//clnsurf(nsrf),"W/m2" ) 1014 CALL histdef2d(iff,clef_stations(iff), & 1015 o_wbilo_srf(nsrf)%flag,o_wbilo_srf(nsrf)%name,"Bilan eau "//clnsurf(nsrf),"kg/(m2*s)") 1016 if (iflag_pbl>1 .and. lev_files(iff).gt.10 ) then 1017 CALL histdef2d(iff,clef_stations(iff), & 1018 o_tke_srf(nsrf)%flag,o_tke_srf(nsrf)%name,"Max Turb. Kinetic Energy "//clnsurf(nsrf),"-") 1019 1020 IF (.NOT.clef_stations(iff)) THEN 1021 ! 1022 !IM: there is no way to have one single value in a netcdf file 1023 ! 1024 type_ecri(1) = 't_max(X)' 1025 type_ecri(2) = 't_max(X)' 1026 type_ecri(3) = 't_max(X)' 1027 type_ecri(4) = 't_max(X)' 1028 type_ecri(5) = 't_max(X)' 1029 type_ecri(6) = 't_max(X)' 1030 CALL histdef2d(iff,clef_stations(iff), & 1031 o_tke_max_srf(nsrf)%flag,o_tke_max_srf(nsrf)%name,"Max Turb. Kinetic Energy "//clnsurf(nsrf),"-") 1032 type_ecri(:) = type_ecri_files(:) 1033 ENDIF 1034 1035 endif 1036 1037 CALL histdef2d(iff,clef_stations(iff), & 1038 o_albe_srf(nsrf)%flag,o_albe_srf(nsrf)%name,"Albedo VIS surf. "//clnsurf(nsrf),"-") 1039 CALL histdef2d(iff,clef_stations(iff), & 1040 o_rugs_srf(nsrf)%flag,o_rugs_srf(nsrf)%name,"Surface roughness "//clnsurf(nsrf),"m") 1041 CALL histdef2d(iff,clef_stations(iff), & 1042 o_ages_srf(nsrf)%flag,o_ages_srf(nsrf)%name,"Snow age", "day") 1043 END DO 1044 1045 IF (new_aod .AND. (.NOT. aerosol_couple)) THEN 1046 IF (ok_ade.OR.ok_aie) THEN 1047 1048 CALL histdef2d(iff,clef_stations(iff), & 1049 o_od550aer%flag,o_od550aer%name, "Total aerosol optical depth at 550nm", "-") 1050 CALL histdef2d(iff,clef_stations(iff), & 1051 o_od865aer%flag,o_od865aer%name, "Total aerosol optical depth at 870nm", "-") 1052 CALL histdef2d(iff,clef_stations(iff), & 1053 o_absvisaer%flag,o_absvisaer%name, "Absorption aerosol visible optical depth", "-") 1054 CALL histdef2d(iff,clef_stations(iff), & 1055 o_od550lt1aer%flag,o_od550lt1aer%name, "Fine mode optical depth", "-") 1056 1057 1058 CALL histdef2d(iff,clef_stations(iff), & 1059 o_sconcso4%flag,o_sconcso4%name,"Surface Concentration of Sulfate ","kg/m3") 1060 CALL histdef2d(iff,clef_stations(iff), & 1061 o_sconcoa%flag,o_sconcoa%name,"Surface Concentration of Organic Aerosol ","kg/m3") 1062 CALL histdef2d(iff,clef_stations(iff), & 1063 o_sconcbc%flag,o_sconcbc%name,"Surface Concentration of Black Carbon ","kg/m3") 1064 CALL histdef2d(iff,clef_stations(iff), & 1065 o_sconcss%flag,o_sconcss%name,"Surface Concentration of Sea Salt ","kg/m3") 1066 CALL histdef2d(iff,clef_stations(iff), & 1067 o_sconcdust%flag,o_sconcdust%name,"Surface Concentration of Dust ","kg/m3") 1068 CALL histdef3d(iff,clef_stations(iff), & 1069 o_concso4%flag,o_concso4%name,"Concentration of Sulfate ","kg/m3") 1070 CALL histdef3d(iff,clef_stations(iff), & 1071 o_concoa%flag,o_concoa%name,"Concentration of Organic Aerosol ","kg/m3") 1072 CALL histdef3d(iff,clef_stations(iff), & 1073 o_concbc%flag,o_concbc%name,"Concentration of Black Carbon ","kg/m3") 1074 CALL histdef3d(iff,clef_stations(iff), & 1075 o_concss%flag,o_concss%name,"Concentration of Sea Salt ","kg/m3") 1076 CALL histdef3d(iff,clef_stations(iff), & 1077 o_concdust%flag,o_concdust%name,"Concentration of Dust ","kg/m3") 1078 CALL histdef2d(iff,clef_stations(iff), & 1079 o_loadso4%flag,o_loadso4%name,"Column Load of Sulfate ","kg/m2") 1080 CALL histdef2d(iff,clef_stations(iff), & 1081 o_loadoa%flag,o_loadoa%name,"Column Load of Organic Aerosol ","kg/m2") 1082 CALL histdef2d(iff,clef_stations(iff), & 1083 o_loadbc%flag,o_loadbc%name,"Column Load of Black Carbon ","kg/m2") 1084 CALL histdef2d(iff,clef_stations(iff), & 1085 o_loadss%flag,o_loadss%name,"Column Load of Sea Salt ","kg/m2") 1086 CALL histdef2d(iff,clef_stations(iff), & 1087 o_loaddust%flag,o_loaddust%name,"Column Load of Dust ","kg/m2") 1088 1089 DO naero = 1, naero_spc 1090 CALL histdef2d(iff,clef_stations(iff), & 1091 o_tausumaero(naero)%flag,o_tausumaero(naero)%name,"Aerosol Optical depth at 550 nm "//name_aero(naero),"1") 1092 END DO 1093 ENDIF 1094 ENDIF 1095 1096 IF (ok_ade) THEN 1097 CALL histdef2d(iff,clef_stations(iff), & 1098 o_topswad%flag,o_topswad%name, "ADE at TOA", "W/m2") 1099 CALL histdef2d(iff,clef_stations(iff), & 1100 o_solswad%flag,o_solswad%name, "ADE at SRF", "W/m2") 1101 1102 CALL histdef2d(iff,clef_stations(iff), & 1103 o_swtoaas_nat%flag,o_swtoaas_nat%name, "Natural aerosol radiative forcing all-sky at TOA", "W/m2") 1104 CALL histdef2d(iff,clef_stations(iff), & 1105 o_swsrfas_nat%flag,o_swsrfas_nat%name, "Natural aerosol radiative forcing all-sky at SRF", "W/m2") 1106 CALL histdef2d(iff,clef_stations(iff), & 1107 o_swtoacs_nat%flag,o_swtoacs_nat%name, "Natural aerosol radiative forcing clear-sky at TOA", "W/m2") 1108 CALL histdef2d(iff,clef_stations(iff), & 1109 o_swsrfcs_nat%flag,o_swsrfcs_nat%name, "Natural aerosol radiative forcing clear-sky at SRF", "W/m2") 1110 1111 CALL histdef2d(iff,clef_stations(iff), & 1112 o_swtoaas_ant%flag,o_swtoaas_ant%name, "Anthropogenic aerosol radiative forcing all-sky at TOA", "W/m2") 1113 CALL histdef2d(iff,clef_stations(iff), & 1114 o_swsrfas_ant%flag,o_swsrfas_ant%name, "Anthropogenic aerosol radiative forcing all-sky at SRF", "W/m2") 1115 CALL histdef2d(iff,clef_stations(iff), & 1116 o_swtoacs_ant%flag,o_swtoacs_ant%name, "Anthropogenic aerosol radiative forcing clear-sky at TOA", "W/m2") 1117 CALL histdef2d(iff,clef_stations(iff), & 1118 o_swsrfcs_ant%flag,o_swsrfcs_ant%name, "Anthropogenic aerosol radiative forcing clear-sky at SRF", "W/m2") 1119 1120 IF (.NOT. aerosol_couple) THEN 1121 CALL histdef2d(iff,clef_stations(iff), & 1122 o_swtoacf_nat%flag,o_swtoacf_nat%name, "Natural aerosol impact on cloud radiative forcing at TOA", "W/m2") 1123 CALL histdef2d(iff,clef_stations(iff), & 1124 o_swsrfcf_nat%flag,o_swsrfcf_nat%name, "Natural aerosol impact on cloud radiative forcing at SRF", "W/m2") 1125 CALL histdef2d(iff,clef_stations(iff), & 1126 o_swtoacf_ant%flag,o_swtoacf_ant%name, "Anthropogenic aerosol impact on cloud radiative forcing at TOA", "W/m2") 1127 CALL histdef2d(iff,clef_stations(iff), & 1128 o_swsrfcf_ant%flag,o_swsrfcf_ant%name, "Anthropogenic aerosol impact on cloud radiative forcing at SRF", "W/m2") 1129 CALL histdef2d(iff,clef_stations(iff), & 1130 o_swtoacf_zero%flag,o_swtoacf_zero%name, "Cloud radiative forcing (allsky-clearsky fluxes) at TOA", "W/m2") 1131 CALL histdef2d(iff,clef_stations(iff), & 1132 o_swsrfcf_zero%flag,o_swsrfcf_zero%name, "Cloud radiative forcing (allsky-clearsky fluxes) at SRF", "W/m2") 1133 ENDIF 1134 1135 ENDIF 1136 1137 IF (ok_aie) THEN 1138 CALL histdef2d(iff,clef_stations(iff), & 1139 o_topswai%flag,o_topswai%name, "AIE at TOA", "W/m2") 1140 CALL histdef2d(iff,clef_stations(iff), & 1141 o_solswai%flag,o_solswai%name, "AIE at SFR", "W/m2") 1142 !Cloud droplet number concentration 1143 CALL histdef3d(iff,clef_stations(iff), & 1144 o_scdnc%flag,o_scdnc%name, "Cloud droplet number concentration","m-3") 1145 CALL histdef2d(iff,clef_stations(iff), & 1146 o_cldncl%flag,o_cldncl%name, "CDNC at top of liquid water cloud", "m-3") 1147 CALL histdef3d(iff,clef_stations(iff), & 1148 o_reffclws%flag,o_reffclws%name, "Stratiform Cloud Droplet Effective Radius (aerosol diags.)","m") 1149 CALL histdef3d(iff,clef_stations(iff), & 1150 o_reffclwc%flag,o_reffclwc%name, "Convective Cloud Droplet Effective Radius (aerosol diags.)","m") 1151 CALL histdef2d(iff,clef_stations(iff), & 1152 o_cldnvi%flag,o_cldnvi%name, "Column Integrated Cloud Droplet Number", "m-2") 1153 CALL histdef3d(iff,clef_stations(iff), & 1154 o_lcc3d%flag,o_lcc3d%name, "Cloud liquid fraction","1") 1155 CALL histdef3d(iff,clef_stations(iff), & 1156 o_lcc3dcon%flag,o_lcc3dcon%name, "Convective cloud liquid fraction","1") 1157 CALL histdef3d(iff,clef_stations(iff), & 1158 o_lcc3dstra%flag,o_lcc3dstra%name, "Stratiform cloud liquid fraction","1") 1159 CALL histdef2d(iff,clef_stations(iff), & 1160 o_lcc%flag,o_lcc%name, "Cloud liquid fraction at top of cloud","1") 1161 CALL histdef2d(iff,clef_stations(iff), & 1162 o_reffclwtop%flag,o_reffclwtop%name, "Droplet effective radius at top of liquid water cloud", "m") 1163 ENDIF 1164 1165 1166 CALL histdef2d(iff,clef_stations(iff), & 1167 o_alb1%flag,o_alb1%name, "Surface VIS albedo", "-") 1168 CALL histdef2d(iff,clef_stations(iff), & 1169 o_alb2%flag,o_alb2%name, "Surface Near IR albedo", "-") 1170 CALL histdef2d(iff,clef_stations(iff), & 1171 o_cdrm%flag,o_cdrm%name, "Momentum drag coef.", "-") 1172 CALL histdef2d(iff,clef_stations(iff), & 1173 o_cdrh%flag,o_cdrh%name, "Heat drag coef.", "-" ) 1174 CALL histdef2d(iff,clef_stations(iff), & 1175 o_cldl%flag,o_cldl%name, "Low-level cloudiness", "-") 1176 CALL histdef2d(iff,clef_stations(iff), & 1177 o_cldm%flag,o_cldm%name, "Mid-level cloudiness", "-") 1178 CALL histdef2d(iff,clef_stations(iff), & 1179 o_cldh%flag,o_cldh%name, "High-level cloudiness", "-") 1180 CALL histdef2d(iff,clef_stations(iff), & 1181 o_cldt%flag,o_cldt%name, "Total cloudiness", "-") 1182 CALL histdef2d(iff,clef_stations(iff), & 1183 o_cldq%flag,o_cldq%name, "Cloud liquid water path", "kg/m2") 1184 CALL histdef2d(iff,clef_stations(iff), & 1185 o_lwp%flag,o_lwp%name, "Cloud water path", "kg/m2") 1186 CALL histdef2d(iff,clef_stations(iff), & 1187 o_iwp%flag,o_iwp%name, "Cloud ice water path", "kg/m2" ) 1188 CALL histdef2d(iff,clef_stations(iff), & 1189 o_ue%flag,o_ue%name, "Zonal energy transport", "-") 1190 CALL histdef2d(iff,clef_stations(iff), & 1191 o_ve%flag,o_ve%name, "Merid energy transport", "-") 1192 CALL histdef2d(iff,clef_stations(iff), & 1193 o_uq%flag,o_uq%name, "Zonal humidity transport", "-") 1194 CALL histdef2d(iff,clef_stations(iff), & 1195 o_vq%flag,o_vq%name, "Merid humidity transport", "-") 1196 1197 IF(iflag_con.GE.3) THEN ! sb 1198 CALL histdef2d(iff,clef_stations(iff), & 1199 o_cape%flag,o_cape%name, "Conv avlbl pot ener", "J/kg") 1200 CALL histdef2d(iff,clef_stations(iff), & 1201 o_pbase%flag,o_pbase%name, "Cld base pressure", "Pa") 1202 CALL histdef2d(iff,clef_stations(iff), & 1203 o_ptop%flag,o_ptop%name, "Cld top pressure", "Pa") 1204 CALL histdef2d(iff,clef_stations(iff), & 1205 o_fbase%flag,o_fbase%name, "Cld base mass flux", "kg/m2/s") 1206 CALL histdef2d(iff,clef_stations(iff), & 1207 o_plcl%flag,o_plcl%name, "Lifting Condensation Level", "hPa") 1208 CALL histdef2d(iff,clef_stations(iff), & 1209 o_plfc%flag,o_plfc%name, "Level of Free Convection", "hPa") 1210 CALL histdef2d(iff,clef_stations(iff), & 1211 o_wbeff%flag,o_wbeff%name, "Conv. updraft velocity at LFC", "m/s") 1212 CALL histdef2d(iff,clef_stations(iff), & 1213 o_prw%flag,o_prw%name, "Precipitable water", "kg/m2") 1214 IF (.NOT.clef_stations(iff)) THEN 1215 ! 1216 !IM: there is no way to have one single value in a netcdf file 1217 ! 1218 type_ecri(1) = 't_max(X)' 1219 type_ecri(2) = 't_max(X)' 1220 type_ecri(3) = 't_max(X)' 1221 type_ecri(4) = 't_max(X)' 1222 type_ecri(5) = 't_max(X)' 1223 type_ecri(6) = 't_max(X)' 1224 CALL histdef2d(iff,clef_stations(iff), & 1225 o_cape_max%flag,o_cape_max%name, "CAPE max.", "J/kg") 1226 ENDIF 1227 type_ecri(:) = type_ecri_files(:) 1228 CALL histdef3d(iff,clef_stations(iff), & 1229 o_upwd%flag,o_upwd%name, "saturated updraft", "kg/m2/s") 1230 CALL histdef3d(iff,clef_stations(iff), & 1231 o_Ma%flag,o_Ma%name, "undilute adiab updraft", "kg/m2/s") 1232 CALL histdef3d(iff,clef_stations(iff), & 1233 o_dnwd%flag,o_dnwd%name, "saturated downdraft", "kg/m2/s") 1234 CALL histdef3d(iff,clef_stations(iff), & 1235 o_dnwd0%flag,o_dnwd0%name, "unsat. downdraft", "kg/m2/s") 1236 CALL histdef3d(iff,clef_stations(iff), & 1237 o_mc%flag,o_mc%name, "Convective mass flux", "kg/m2/s") 1238 type_ecri(1) = 'inst(X)' 1239 type_ecri(2) = 'inst(X)' 1240 type_ecri(3) = 'inst(X)' 1241 type_ecri(4) = 'inst(X)' 1242 type_ecri(5) = 'inst(X)' 1243 type_ecri(6) = 'inst(X)' 1244 CALL histdef2d(iff,clef_stations(iff), & 1245 o_ftime_con%flag,o_ftime_con%name, "Fraction of time convection Occurs", " ") 1246 type_ecri(:) = type_ecri_files(:) 1247 ENDIF !iflag_con .GE. 3 1248 1249 CALL histdef2d(iff,clef_stations(iff), & 1250 o_s_pblh%flag,o_s_pblh%name, "Boundary Layer Height", "m") 1251 CALL histdef2d(iff,clef_stations(iff), & 1252 o_s_pblt%flag,o_s_pblt%name, "t at Boundary Layer Height", "K") 1253 CALL histdef2d(iff,clef_stations(iff), & 1254 o_s_lcl%flag,o_s_lcl%name, "Condensation level", "m") 1255 CALL histdef2d(iff,clef_stations(iff), & 1256 o_s_therm%flag,o_s_therm%name, "Exces du thermique", "K") 1257 !IM : Les champs suivants (s_oliqCL, s_cteiCL, s_trmb1, s_trmb2, s_trmb3) ne sont pas definis dans HBTM.F 1258 !CALL histdef2d(iff,clef_stations(iff), & 1259 !o_s_capCL%flag,o_s_capCL%name, "Conv avlbl pot enerfor ABL", "J/m2" ) 1260 !CALL histdef2d(iff,clef_stations(iff), & 1261 !o_s_oliqCL%flag,o_s_oliqCL%name, "Liq Water in BL", "kg/m2") 1262 !CALL histdef2d(iff,clef_stations(iff), & 1263 !o_s_cteiCL%flag,o_s_cteiCL%name, "Instability criteria(ABL)", "K") 1264 !CALL histdef2d(iff,clef_stations(iff), & 1265 !o_s_trmb1%flag,o_s_trmb1%name, "deep_cape(HBTM2)", "J/m2") 1266 !CALL histdef2d(iff,clef_stations(iff), & 1267 !o_s_trmb2%flag,o_s_trmb2%name, "inhibition (HBTM2)", "J/m2") 1268 !CALL histdef2d(iff,clef_stations(iff), & 1269 !o_s_trmb3%flag,o_s_trmb3%name, "Point Omega (HBTM2)", "m") 1270 1271 ! Champs interpolles sur des niveaux de pression 1272 1273 type_ecri(1) = 'inst(X)' 1274 type_ecri(2) = 'inst(X)' 1275 type_ecri(3) = 'inst(X)' 1276 type_ecri(4) = 'inst(X)' 1277 type_ecri(5) = 'inst(X)' 1278 type_ecri(6) = 'inst(X)' 1279 1280 ! Attention a reverifier 1281 1282 ilev=0 1283 DO k=1, nlevSTD 1284 bb2=clevSTD(k) 1285 IF(bb2.EQ."850".OR.bb2.EQ."700".OR.bb2.EQ."500".OR.bb2.EQ."200" & 1286 .OR.bb2.EQ."100".OR.bb2.EQ."50".OR.bb2.EQ."10")THEN 1287 ilev=ilev+1 1288 ! print*,'ilev k bb2 flag name ',ilev,k, bb2,o_uSTDlevs(ilev)%flag,o_uSTDlevs(ilev)%name 1289 CALL histdef2d(iff,clef_stations(iff), & 1290 o_uSTDlevs(ilev)%flag,o_uSTDlevs(ilev)%name,"Zonal wind "//bb2//"hPa", "m/s") 1291 CALL histdef2d(iff,clef_stations(iff), & 1292 o_vSTDlevs(ilev)%flag,o_vSTDlevs(ilev)%name,"Meridional wind "//bb2//"hPa", "m/s") 1293 CALL histdef2d(iff,clef_stations(iff), & 1294 o_wSTDlevs(ilev)%flag,o_wSTDlevs(ilev)%name,"Vertical wind "//bb2//"hPa", "Pa/s") 1295 CALL histdef2d(iff,clef_stations(iff), & 1296 o_zSTDlevs(ilev)%flag,o_zSTDlevs(ilev)%name,"Geopotential height "//bb2//"hPa", "m") 1297 CALL histdef2d(iff,clef_stations(iff), & 1298 o_qSTDlevs(ilev)%flag,o_qSTDlevs(ilev)%name,"Specific humidity "//bb2//"hPa", "kg/kg" ) 1299 CALL histdef2d(iff,clef_stations(iff), & 1300 o_tSTDlevs(ilev)%flag,o_tSTDlevs(ilev)%name,"Temperature "//bb2//"hPa", "K") 1301 ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR."500".OR.bb2.EQ."200".OR.bb2.EQ."50".OR.bb2.EQ."10") 1302 ENDDO 1303 type_ecri(:) = type_ecri_files(:) 1304 1305 CALL histdef2d(iff,clef_stations(iff), & 1306 o_t_oce_sic%flag,o_t_oce_sic%name, "Temp mixte oce-sic", "K") 1307 1308 IF (type_ocean=='slab') & 1309 CALL histdef2d(iff,clef_stations(iff), & 1310 o_slab_bils%flag, o_slab_bils%name,"Bilan au sol sur ocean slab", "W/m2") 1311 1312 ! Couplage conv-CL 1313 IF (iflag_con.GE.3) THEN 1314 IF (iflag_coupl>=1) THEN 1315 CALL histdef2d(iff,clef_stations(iff), & 1316 o_ale_bl%flag,o_ale_bl%name, "ALE BL", "m2/s2") 1317 CALL histdef2d(iff,clef_stations(iff), & 1318 o_alp_bl%flag,o_alp_bl%name, "ALP BL", "m2/s2") 1319 ENDIF 1320 ENDIF !(iflag_con.GE.3) 1321 1322 CALL histdef2d(iff,clef_stations(iff), & 1323 o_weakinv%flag,o_weakinv%name, "Weak inversion", "-") 1324 CALL histdef2d(iff,clef_stations(iff), & 1325 o_dthmin%flag,o_dthmin%name, "dTheta mini", "K/m") 1326 CALL histdef2d(iff,clef_stations(iff), & 1327 o_rh2m%flag,o_rh2m%name, "Relative humidity at 2m", "%" ) 1328 1329 IF (.NOT.clef_stations(iff)) THEN 1330 ! 1331 !IM: there is no way to have one single value in a netcdf file 1332 ! 1333 type_ecri(1) = 't_min(X)' 1334 type_ecri(2) = 't_min(X)' 1335 type_ecri(3) = 't_min(X)' 1336 type_ecri(4) = 't_min(X)' 1337 type_ecri(5) = 't_min(X)' 1338 type_ecri(6) = 't_min(X)' 1339 CALL histdef2d(iff,clef_stations(iff),o_rh2m_min%flag,o_rh2m_min%name, "Min Relative humidity at 2m", "%" ) 1340 type_ecri(1) = 't_max(X)' 1341 type_ecri(2) = 't_max(X)' 1342 type_ecri(3) = 't_max(X)' 1343 type_ecri(4) = 't_max(X)' 1344 type_ecri(5) = 't_max(X)' 1345 type_ecri(6) = 't_max(X)' 1346 CALL histdef2d(iff,clef_stations(iff),o_rh2m_max%flag,o_rh2m_max%name, "Max Relative humidity at 2m", "%" ) 1347 ENDIF 1348 1349 type_ecri(:) = type_ecri_files(:) 1350 CALL histdef2d(iff,clef_stations(iff),o_qsat2m%flag,o_qsat2m%name, "Saturant humidity at 2m", "%") 1351 CALL histdef2d(iff,clef_stations(iff),o_tpot%flag,o_tpot%name, "Surface air potential temperature", "K") 1352 CALL histdef2d(iff,clef_stations(iff), & 1353 o_tpote%flag,o_tpote%name, "Surface air equivalent potential temperature", "K") 1354 CALL histdef2d(iff,clef_stations(iff),o_SWnetOR%flag,o_SWnetOR%name, "Sfce net SW radiation OR", "W/m2") 1355 CALL histdef2d(iff,clef_stations(iff),o_SWdownOR%flag,o_SWdownOR%name, "Sfce incident SW radiation OR", "W/m2") 1356 CALL histdef2d(iff,clef_stations(iff),o_LWdownOR%flag,o_LWdownOR%name, "Sfce incident LW radiation OR", "W/m2") 1357 CALL histdef2d(iff,clef_stations(iff),o_snowl%flag,o_snowl%name, "Solid Large-scale Precip.", "kg/(m2*s)") 1358 1359 CALL histdef2d(iff,clef_stations(iff),o_solldown%flag,o_solldown%name, "Down. IR rad. at surface", "W/m2") 1360 CALL histdef2d(iff,clef_stations(iff),o_dtsvdfo%flag,o_dtsvdfo%name, "Boundary-layer dTs(o)", "K/s") 1361 CALL histdef2d(iff,clef_stations(iff),o_dtsvdft%flag,o_dtsvdft%name, "Boundary-layer dTs(t)", "K/s") 1362 CALL histdef2d(iff,clef_stations(iff),o_dtsvdfg%flag,o_dtsvdfg%name, "Boundary-layer dTs(g)", "K/s") 1363 CALL histdef2d(iff,clef_stations(iff),o_dtsvdfi%flag,o_dtsvdfi%name, "Boundary-layer dTs(g)", "K/s") 1364 CALL histdef2d(iff,clef_stations(iff),o_rugs%flag,o_rugs%name, "rugosity", "-" ) 1365 1366 ! Champs 3D: 1367 CALL histdef3d(iff,clef_stations(iff),o_ec550aer%flag,o_ec550aer%name, "Extinction at 550nm", "m^-1") 1368 CALL histdef3d(iff,clef_stations(iff),o_lwcon%flag,o_lwcon%name, "Cloud liquid water content", "kg/kg") 1369 CALL histdef3d(iff,clef_stations(iff),o_iwcon%flag,o_iwcon%name, "Cloud ice water content", "kg/kg") 1370 CALL histdef3d(iff,clef_stations(iff),o_temp%flag,o_temp%name, "Air temperature", "K" ) 1371 CALL histdef3d(iff,clef_stations(iff),o_theta%flag,o_theta%name, "Potential air temperature", "K" ) 1372 CALL histdef3d(iff,clef_stations(iff),o_ovap%flag,o_ovap%name, "Specific humidity", "kg/kg" ) 1373 CALL histdef3d(iff,clef_stations(iff), & 1374 o_ovapinit%flag,o_ovapinit%name, "Specific humidity (begin of timestep)", "kg/kg" ) 1375 CALL histdef3d(iff,clef_stations(iff), & 1376 o_geop%flag,o_geop%name, "Geopotential height", "m2/s2") 1377 CALL histdef3d(iff,clef_stations(iff), & 1378 o_vitu%flag,o_vitu%name, "Zonal wind", "m/s" ) 1379 CALL histdef3d(iff,clef_stations(iff), & 1380 o_vitv%flag,o_vitv%name, "Meridional wind", "m/s" ) 1381 CALL histdef3d(iff,clef_stations(iff), & 1382 o_vitw%flag,o_vitw%name, "Vertical wind", "Pa/s" ) 1383 CALL histdef3d(iff,clef_stations(iff), & 1384 o_pres%flag,o_pres%name, "Air pressure", "Pa" ) 1385 CALL histdef3d(iff,clef_stations(iff), & 1386 o_paprs%flag,o_paprs%name, "Air pressure Inter-Couches", "Pa" ) 1387 CALL histdef3d(iff,clef_stations(iff), & 1388 o_mass%flag,o_mass%name, "Masse Couches", "kg/m2" ) 1389 CALL histdef3d(iff,clef_stations(iff), & 1390 o_zfull%flag,o_zfull%name, "Altitude of full pressure levels", "m" ) 1391 CALL histdef3d(iff,clef_stations(iff), & 1392 o_zhalf%flag,o_zhalf%name, "Altitude of half pressure levels", "m" ) 1393 CALL histdef3d(iff,clef_stations(iff), & 1394 o_rneb%flag,o_rneb%name, "Cloud fraction", "-") 1395 CALL histdef3d(iff,clef_stations(iff), & 1396 o_rnebcon%flag,o_rnebcon%name, "Convective Cloud Fraction", "-") 1397 CALL histdef3d(iff,clef_stations(iff), & 1398 o_rhum%flag,o_rhum%name, "Relative humidity", "-") 1399 CALL histdef3d(iff,clef_stations(iff), & 1400 o_ozone%flag,o_ozone%name, "Ozone mole fraction", "-") 1401 if (read_climoz == 2) & 1402 CALL histdef3d(iff,clef_stations(iff), & 1403 o_ozone_light%flag,o_ozone_light%name, & 1404 "Daylight ozone mole fraction", "-") 1405 CALL histdef3d(iff,clef_stations(iff), & 1406 o_dtphy%flag,o_dtphy%name, "Physics dT", "K/s") 1407 CALL histdef3d(iff,clef_stations(iff), & 1408 o_dqphy%flag,o_dqphy%name, "Physics dQ", "(kg/kg)/s") 1409 CALL histdef3d(iff,clef_stations(iff), & 1410 o_cldtau%flag,o_cldtau%name, "Cloud optical thickness", "1") 1411 CALL histdef3d(iff,clef_stations(iff), & 1412 o_cldemi%flag,o_cldemi%name, "Cloud optical emissivity", "1") 1413 !IM: bug ?? dimensionnement variables (klon,klev+1) pmflxr, pmflxs, prfl, psfl 1414 CALL histdef3d(iff,clef_stations(iff), & 1415 o_pr_con_l%flag,o_pr_con_l%name, "Convective precipitation lic", " ") 1416 CALL histdef3d(iff,clef_stations(iff), & 1417 o_pr_con_i%flag,o_pr_con_i%name, "Convective precipitation ice", " ") 1418 CALL histdef3d(iff,clef_stations(iff), & 1419 o_pr_lsc_l%flag,o_pr_lsc_l%name, "Large scale precipitation lic", " ") 1420 CALL histdef3d(iff,clef_stations(iff), & 1421 o_pr_lsc_i%flag,o_pr_lsc_i%name, "Large scale precipitation ice", " ") 1422 !Cloud droplet effective radius 1423 CALL histdef3d(iff,clef_stations(iff), & 1424 o_re%flag,o_re%name, "Cloud droplet effective radius","um") 1425 CALL histdef3d(iff,clef_stations(iff), & 1426 o_fl%flag,o_fl%name, "Denominator of Cloud droplet effective radius"," ") 1427 !FH Sorties pour la couche limite 1428 if (iflag_pbl>1) then 1429 CALL histdef3d(iff,clef_stations(iff), & 1430 o_tke%flag,o_tke%name, "TKE", "m2/s2") 1431 IF (.NOT.clef_stations(iff)) THEN 1432 ! 1433 !IM: there is no way to have one single value in a netcdf file 1434 ! 1435 type_ecri(1) = 't_max(X)' 1436 type_ecri(2) = 't_max(X)' 1437 type_ecri(3) = 't_max(X)' 1438 type_ecri(4) = 't_max(X)' 1439 type_ecri(5) = 't_max(X)' 1440 type_ecri(6) = 't_max(X)' 1441 CALL histdef3d(iff,clef_stations(iff), & 1442 o_tke_max%flag,o_tke_max%name, "TKE max", "m2/s2") 1443 ENDIF 1444 type_ecri(:) = type_ecri_files(:) 1445 endif 1446 1447 CALL histdef3d(iff,clef_stations(iff), & 1448 o_kz%flag,o_kz%name, "Kz melange", "m2/s") 1449 IF (.NOT.clef_stations(iff)) THEN 1450 ! 1451 !IM: there is no way to have one single value in a netcdf file 1452 ! 1453 type_ecri(1) = 't_max(X)' 1454 type_ecri(2) = 't_max(X)' 1455 type_ecri(3) = 't_max(X)' 1456 type_ecri(4) = 't_max(X)' 1457 type_ecri(5) = 't_max(X)' 1458 type_ecri(6) = 't_max(X)' 1459 CALL histdef3d(iff,clef_stations(iff), & 1460 o_kz_max%flag,o_kz_max%name, "Kz melange max", "m2/s" ) 1461 ENDIF 1462 type_ecri(:) = type_ecri_files(:) 1463 CALL histdef3d(iff,clef_stations(iff), & 1464 o_clwcon%flag,o_clwcon%name, "Convective Cloud Liquid water content", "kg/kg") 1465 CALL histdef3d(iff,clef_stations(iff), & 1466 o_dtdyn%flag,o_dtdyn%name, "Dynamics dT", "K/s") 1467 CALL histdef3d(iff,clef_stations(iff), & 1468 o_dqdyn%flag,o_dqdyn%name, "Dynamics dQ", "(kg/kg)/s") 1469 CALL histdef3d(iff,clef_stations(iff), & 1470 o_dudyn%flag,o_dudyn%name, "Dynamics dU", "m/s2") 1471 CALL histdef3d(iff,clef_stations(iff), & 1472 o_dvdyn%flag,o_dvdyn%name, "Dynamics dV", "m/s2") 1473 CALL histdef3d(iff,clef_stations(iff), & 1474 o_dtcon%flag,o_dtcon%name, "Convection dT", "K/s") 1475 CALL histdef3d(iff,clef_stations(iff), & 1476 o_ducon%flag,o_ducon%name, "Convection du", "m/s2") 1477 CALL histdef3d(iff,clef_stations(iff), & 1478 o_dqcon%flag,o_dqcon%name, "Convection dQ", "(kg/kg)/s") 1479 1480 ! Wakes 1481 IF(iflag_con.EQ.3) THEN 1482 IF (iflag_wake >= 1) THEN 1483 CALL histdef2d(iff,clef_stations(iff), & 1484 o_ale_wk%flag,o_ale_wk%name, "ALE WK", "m2/s2") 1485 CALL histdef2d(iff,clef_stations(iff), & 1486 o_alp_wk%flag,o_alp_wk%name, "ALP WK", "m2/s2") 1487 CALL histdef2d(iff,clef_stations(iff), & 1488 o_ale%flag,o_ale%name, "ALE", "m2/s2") 1489 CALL histdef2d(iff,clef_stations(iff), & 1490 o_alp%flag,o_alp%name, "ALP", "W/m2") 1491 CALL histdef2d(iff,clef_stations(iff),o_cin%flag,o_cin%name, "Convective INhibition", "m2/s2") 1492 CALL histdef2d(iff,clef_stations(iff),o_wape%flag,o_WAPE%name, "WAPE", "m2/s2") 1493 CALL histdef2d(iff,clef_stations(iff),o_wake_h%flag,o_wake_h%name, "wake_h", "-") 1494 CALL histdef2d(iff,clef_stations(iff),o_wake_s%flag,o_wake_s%name, "wake_s", "-") 1495 CALL histdef3d(iff,clef_stations(iff),o_dtwak%flag,o_dtwak%name, "Wake dT", "K/s") 1496 CALL histdef3d(iff,clef_stations(iff),o_dqwak%flag,o_dqwak%name, "Wake dQ", "(kg/kg)/s") 1497 CALL histdef3d(iff,clef_stations(iff),o_wake_deltat%flag,o_wake_deltat%name, "wake_deltat", " ") 1498 CALL histdef3d(iff,clef_stations(iff),o_wake_deltaq%flag,o_wake_deltaq%name, "wake_deltaq", " ") 1499 CALL histdef3d(iff,clef_stations(iff),o_wake_omg%flag,o_wake_omg%name, "wake_omg", " ") 1500 ENDIF 1501 CALL histdef3d(iff,clef_stations(iff),o_Vprecip%flag,o_Vprecip%name, "precipitation vertical profile", "-") 1502 CALL histdef3d(iff,clef_stations(iff),o_ftd%flag,o_ftd%name, "tend temp due aux descentes precip", "-") 1503 CALL histdef3d(iff,clef_stations(iff),o_fqd%flag,o_fqd%name,"tend vap eau due aux descentes precip", "-") 1504 ENDIF !(iflag_con.EQ.3) 1505 1506 CALL histdef3d(iff,clef_stations(iff),o_dtlsc%flag,o_dtlsc%name, "Condensation dT", "K/s") 1507 CALL histdef3d(iff,clef_stations(iff),o_dtlschr%flag,o_dtlschr%name,"Large-scale condensational heating rate","K/s") 1508 CALL histdef3d(iff,clef_stations(iff),o_dqlsc%flag,o_dqlsc%name, "Condensation dQ", "(kg/kg)/s") 1509 CALL histdef3d(iff,clef_stations(iff),o_dtvdf%flag,o_dtvdf%name, "Boundary-layer dT", "K/s") 1510 CALL histdef3d(iff,clef_stations(iff),o_dqvdf%flag,o_dqvdf%name, "Boundary-layer dQ", "(kg/kg)/s") 1511 CALL histdef3d(iff,clef_stations(iff),o_dteva%flag,o_dteva%name, "Reevaporation dT", "K/s") 1512 CALL histdef3d(iff,clef_stations(iff),o_dqeva%flag,o_dqeva%name, "Reevaporation dQ", "(kg/kg)/s") 1513 CALL histdef3d(iff,clef_stations(iff),o_ptconv%flag,o_ptconv%name, "POINTS CONVECTIFS", " ") 1514 CALL histdef3d(iff,clef_stations(iff),o_ratqs%flag,o_ratqs%name, "RATQS", " ") 1515 CALL histdef3d(iff,clef_stations(iff),o_dtthe%flag,o_dtthe%name, "Thermal dT", "K/s") 1516 1517 if(iflag_thermals.gt.1) THEN 1518 CALL histdef3d(iff,clef_stations(iff),o_dqlscth%flag,o_dqlscth%name, "dQ therm.", "(kg/kg)/s") 1519 CALL histdef3d(iff,clef_stations(iff),o_dqlscst%flag,o_dqlscst%name, "dQ strat.", "(kg/kg)/s") 1520 CALL histdef3d(iff,clef_stations(iff),o_dtlscth%flag,o_dtlscth%name, "dQ therm.", "K/s") 1521 CALL histdef3d(iff,clef_stations(iff),o_dtlscst%flag,o_dtlscst%name, "dQ strat.", "K/s") 1522 CALL histdef2d(iff,clef_stations(iff),o_plulth%flag,o_plulth%name, "Rainfall therm.", "K/s") 1523 CALL histdef2d(iff,clef_stations(iff),o_plulst%flag,o_plulst%name, "Rainfall strat.", "K/s") 1524 CALL histdef2d(iff,clef_stations(iff),o_lmaxth%flag,o_lmaxth%name, "Upper level thermals", "") 1525 CALL histdef3d(iff,clef_stations(iff),o_ptconvth%flag,o_ptconvth%name, "POINTS CONVECTIFS therm.", " ") 1526 CALL histdef3d(iff,clef_stations(iff),o_f_th%flag,o_f_th%name, "Thermal plume mass flux", "kg/(m2*s)") 1527 CALL histdef3d(iff,clef_stations(iff),o_e_th%flag,o_e_th%name,"Thermal plume entrainment","K/s") 1528 CALL histdef3d(iff,clef_stations(iff),o_w_th%flag,o_w_th%name,"Thermal plume vertical velocity","m/s") 1529 CALL histdef3d(iff,clef_stations(iff), & 1530 o_lambda_th%flag,o_lambda_th%name,"Thermal plume vertical velocity","m/s") 1531 CALL histdef2d(iff,clef_stations(iff), & 1532 o_ftime_th%flag,o_ftime_th%name,"Fraction of time Shallow convection occurs"," ") 1533 CALL histdef3d(iff,clef_stations(iff), & 1534 o_q_th%flag,o_q_th%name, "Thermal plume total humidity", "kg/kg") 1535 CALL histdef3d(iff,clef_stations(iff), & 1536 o_a_th%flag,o_a_th%name, "Thermal plume fraction", "") 1537 CALL histdef3d(iff,clef_stations(iff), & 1538 o_d_th%flag,o_d_th%name, "Thermal plume detrainment", "K/s") 1539 1540 CALL histdef2d(iff,clef_stations(iff), & 1541 o_f0_th%flag,o_f0_th%name, "Thermal closure mass flux", "K/s") 1542 CALL histdef2d(iff,clef_stations(iff), & 1543 o_zmax_th%flag,o_zmax_th%name, "Thermal plume height", "K/s") 1544 CALL histdef3d(iff,clef_stations(iff), & 1545 o_dqthe%flag,o_dqthe%name, "Thermal dQ", "(kg/kg)/s") 1546 endif !iflag_thermals.gt.1 1547 CALL histdef3d(iff,clef_stations(iff), & 1548 o_dtajs%flag,o_dtajs%name, "Dry adjust. dT", "K/s") 1549 CALL histdef3d(iff,clef_stations(iff), & 1550 o_dqajs%flag,o_dqajs%name, "Dry adjust. dQ", "(kg/kg)/s") 1551 CALL histdef3d(iff,clef_stations(iff), & 1552 o_dtswr%flag,o_dtswr%name, "SW radiation dT", "K/s") 1553 CALL histdef3d(iff,clef_stations(iff), & 1554 o_dtsw0%flag,o_dtsw0%name, "CS SW radiation dT", "K/s") 1555 CALL histdef3d(iff,clef_stations(iff), & 1556 o_dtlwr%flag,o_dtlwr%name, "LW radiation dT", "K/s") 1557 CALL histdef3d(iff,clef_stations(iff), & 1558 o_dtlw0%flag,o_dtlw0%name, "CS LW radiation dT", "K/s") 1559 CALL histdef3d(iff,clef_stations(iff), & 1560 o_dtec%flag,o_dtec%name, "Cinetic dissip dT", "K/s") 1561 CALL histdef3d(iff,clef_stations(iff), & 1562 o_duvdf%flag,o_duvdf%name, "Boundary-layer dU", "m/s2") 1563 CALL histdef3d(iff,clef_stations(iff), & 1564 o_dvvdf%flag,o_dvvdf%name, "Boundary-layer dV", "m/s2") 1565 1566 IF (ok_orodr) THEN 1567 CALL histdef3d(iff,clef_stations(iff), & 1568 o_duoro%flag,o_duoro%name, "Orography dU", "m/s2") 1569 CALL histdef3d(iff,clef_stations(iff), & 1570 o_dvoro%flag,o_dvoro%name, "Orography dV", "m/s2") 1571 CALL histdef3d(iff,clef_stations(iff), & 1572 o_dtoro%flag,o_dtoro%name, "Orography dT", "K/s") 1573 ENDIF 1574 1575 IF (ok_orolf) THEN 1576 CALL histdef3d(iff,clef_stations(iff), & 1577 o_dulif%flag,o_dulif%name, "Orography dU", "m/s2") 1578 CALL histdef3d(iff,clef_stations(iff), & 1579 o_dvlif%flag,o_dvlif%name, "Orography dV", "m/s2") 1580 CALL histdef3d(iff,clef_stations(iff), & 1581 o_dtlif%flag,o_dtlif%name, "Orography dT", "K/s") 1582 ENDIF 1583 1584 IF (ok_hines) then 1585 CALL histdef3d(iff,clef_stations(iff), & 1586 o_duhin%flag,o_duhin%name, "Hines GWD dU", "m/s2") 1587 CALL histdef3d(iff,clef_stations(iff), & 1588 o_dvhin%flag,o_dvhin%name, "Hines GWD dV", "m/s2") 1589 1590 CALL histdef3d(iff,clef_stations(iff), & 1591 o_dthin%flag,o_dthin%name, "Hines GWD dT", "K/s") 1592 ENDIF 1593 1594 CALL histdef3d(iff,clef_stations(iff), & 1595 o_rsu%flag,o_rsu%name, "SW upward radiation", "W m-2") 1596 CALL histdef3d(iff,clef_stations(iff), & 1597 o_rsd%flag,o_rsd%name, "SW downward radiation", "W m-2") 1598 CALL histdef3d(iff,clef_stations(iff), & 1599 o_rlu%flag,o_rlu%name, "LW upward radiation", "W m-2") 1600 CALL histdef3d(iff,clef_stations(iff), & 1601 o_rld%flag,o_rld%name, "LW downward radiation", "W m-2") 1602 1603 CALL histdef3d(iff,clef_stations(iff), & 1604 o_rsucs%flag,o_rsucs%name, "SW CS upward radiation", "W m-2") 1605 CALL histdef3d(iff,clef_stations(iff), & 1606 o_rsdcs%flag,o_rsdcs%name, "SW CS downward radiation", "W m-2") 1607 CALL histdef3d(iff,clef_stations(iff), & 1608 o_rlucs%flag,o_rlucs%name, "LW CS upward radiation", "W m-2") 1609 CALL histdef3d(iff,clef_stations(iff), & 1610 o_rldcs%flag,o_rldcs%name, "LW CS downward radiation", "W m-2") 1611 1612 CALL histdef3d(iff,clef_stations(iff), & 1613 o_tnt%flag,o_tnt%name, "Tendency of air temperature", "K s-1") 1614 1615 CALL histdef3d(iff,clef_stations(iff), & 1616 o_tntc%flag,o_tntc%name, "Tendency of air temperature due to Moist Convection", & 1617 "K s-1") 1618 1619 CALL histdef3d(iff,clef_stations(iff), & 1620 o_tntr%flag,o_tntr%name, "Air temperature tendency due to Radiative heating", & 1621 "K s-1") 1622 1623 CALL histdef3d(iff,clef_stations(iff), & 1624 o_tntscpbl%flag,o_tntscpbl%name, "Air temperature tendency due to St cloud and precipitation and BL mixing", & 1625 "K s-1") 1626 1627 CALL histdef3d(iff,clef_stations(iff), & 1628 o_tnhus%flag,o_tnhus%name, "Tendency of specific humidity", "s-1") 1629 1630 CALL histdef3d(iff,clef_stations(iff), & 1631 o_tnhusc%flag,o_tnhusc%name, "Tendency of specific humidity due to convection", "s-1") 1632 1633 CALL histdef3d(iff,clef_stations(iff), & 1634 o_tnhusscpbl%flag,o_tnhusscpbl%name, "Tendency of Specific humidity due to ST cl, precip and BL mixing", & 1635 "s-1") 1636 1637 CALL histdef3d(iff,clef_stations(iff), & 1638 o_evu%flag,o_evu%name, "Eddy viscosity coefficient for Momentum Variables", "m2 s-1") 1639 1640 CALL histdef3d(iff,clef_stations(iff), & 1641 o_h2o%flag,o_h2o%name, "Mass Fraction of Water", "1") 1642 1643 CALL histdef3d(iff,clef_stations(iff), & 1644 o_mcd%flag,o_mcd%name, "Downdraft COnvective Mass Flux", "kg/(m2*s)") 1645 1646 CALL histdef3d(iff,clef_stations(iff), & 1647 o_dmc%flag,o_dmc%name, "Deep COnvective Mass Flux", "kg/(m2*s)") 1648 1649 CALL histdef3d(iff,clef_stations(iff), & 1650 o_ref_liq%flag,o_ref_liq%name, "Effective radius of convective cloud liquid water particle", "m") 1651 1652 CALL histdef3d(iff,clef_stations(iff), & 1653 o_ref_ice%flag,o_ref_ice%name, "Effective radius of startiform cloud ice particle", "m") 1654 1655 if (RCO2_per.NE.RCO2_act.OR.RCH4_per.NE.RCH4_act.OR. & 1656 RN2O_per.NE.RN2O_act.OR.RCFC11_per.NE.RCFC11_act.OR. & 1657 RCFC12_per.NE.RCFC12_act) THEN 1658 1659 CALL histdef2d(iff,clef_stations(iff),o_rsut4co2%flag,o_rsut4co2%name, & 1660 "TOA Out SW in 4xCO2 atmosphere", "W/m2") 1661 CALL histdef2d(iff,clef_stations(iff),o_rlut4co2%flag,o_rlut4co2%name, & 1662 "TOA Out LW in 4xCO2 atmosphere", "W/m2") 1663 CALL histdef2d(iff,clef_stations(iff),o_rsutcs4co2%flag,o_rsutcs4co2%name, & 1664 "TOA Out CS SW in 4xCO2 atmosphere", "W/m2") 1665 CALL histdef2d(iff,clef_stations(iff),o_rlutcs4co2%flag,o_rlutcs4co2%name, & 1666 "TOA Out CS LW in 4xCO2 atmosphere", "W/m2") 1667 1668 CALL histdef3d(iff,clef_stations(iff),o_rsu4co2%flag,o_rsu4co2%name, & 1669 "Upwelling SW 4xCO2 atmosphere", "W/m2") 1670 CALL histdef3d(iff,clef_stations(iff),o_rlu4co2%flag,o_rlu4co2%name, & 1671 "Upwelling LW 4xCO2 atmosphere", "W/m2") 1672 CALL histdef3d(iff,clef_stations(iff),o_rsucs4co2%flag,o_rsucs4co2%name, & 1673 "Upwelling CS SW 4xCO2 atmosphere", "W/m2") 1674 CALL histdef3d(iff,clef_stations(iff),o_rlucs4co2%flag,o_rlucs4co2%name, & 1675 "Upwelling CS LW 4xCO2 atmosphere", "W/m2") 1676 1677 CALL histdef3d(iff,clef_stations(iff),o_rsd4co2%flag,o_rsd4co2%name, & 1678 "Downwelling SW 4xCO2 atmosphere", "W/m2") 1679 CALL histdef3d(iff,clef_stations(iff),o_rld4co2%flag,o_rld4co2%name, & 1680 "Downwelling LW 4xCO2 atmosphere", "W/m2") 1681 CALL histdef3d(iff,clef_stations(iff),o_rsdcs4co2%flag,o_rsdcs4co2%name, & 1682 "Downwelling CS SW 4xCO2 atmosphere", "W/m2") 1683 CALL histdef3d(iff,clef_stations(iff),o_rldcs4co2%flag,o_rldcs4co2%name, & 1684 "Downwelling CS LW 4xCO2 atmosphere", "W/m2") 1685 1686 endif 1687 1688 1689 IF (nqtot>=3) THEN 1690 DO iq=3,nqtot 1691 iiq=niadv(iq) 1692 o_trac(iq-2) = ctrl_out((/ 4, 5, 1, 1, 1, 10 /),tname(iiq)) 1693 CALL histdef3d (iff,clef_stations(iff), & 1694 o_trac(iq-2)%flag,o_trac(iq-2)%name,'Tracer '//ttext(iiq), "-" ) 1695 ENDDO 1696 ENDIF 1697 1698 CALL histend(nid_files(iff)) 1699 1700 ndex2d = 0 1701 ndex3d = 0 1702 1703 ENDIF ! clef_files 1704 1705 ENDDO ! iff 1706 print*,'Fin phys_output_mod.F90' 1707 end subroutine phys_output_open 1708 1709 SUBROUTINE histdef2d (iff,lpoint,flag_var,nomvar,titrevar,unitvar) 1710 1711 use ioipsl 1712 USE dimphy 1713 USE mod_phys_lmdz_para 1714 USE iophy 1715 1716 IMPLICIT NONE 1717 1718 include "dimensions.h" 1719 include "temps.h" 1720 include "indicesol.h" 1721 include "clesphys.h" 1722 1723 integer :: iff 1724 logical :: lpoint 1725 integer, dimension(nfiles) :: flag_var 1726 character(len=20) :: nomvar 1727 character(len=*) :: titrevar 1728 character(len=*) :: unitvar 1729 1730 real zstophym 1731 1732 if (type_ecri(iff)=='inst(X)'.OR.type_ecri(iff)=='once') then 1733 zstophym=zoutm(iff) 1734 else 1735 zstophym=zdtime 1736 endif 1737 1738 ! Appel a la lecture des noms et niveau d'ecriture des variables dans output.def 1739 call conf_physoutputs(nomvar,flag_var) 1740 1741 if(.NOT.lpoint) THEN 899 CALL histdef2d(iff,clef_stations(iff),o_phis%flag,o_phis%name,"Surface geop.height", "m2/s2" ) 900 CALL histdef2d(iff,clef_stations(iff),o_contfracOR%flag,o_contfracOR%name,"% sfce terre OR", "-" ) 901 CALL histdef2d(iff,clef_stations(iff),o_aireTER%flag,o_aireTER%name,"Grid area CONT", "-" ) 902 CALL histdef2d(iff,clef_stations(iff),o_flat%flag,o_flat%name, "Latent heat flux", "W/m2") 903 CALL histdef2d(iff,clef_stations(iff),o_slp%flag,o_slp%name, "Sea Level Pressure", "Pa" ) 904 CALL histdef2d(iff,clef_stations(iff),o_tsol%flag,o_tsol%name, "Surface Temperature", "K") 905 CALL histdef2d(iff,clef_stations(iff),o_t2m%flag,o_t2m%name, "Temperature 2m", "K" ) 906 IF (.NOT.clef_stations(iff)) THEN 907 ! 908 !IM: there is no way to have one single value in a netcdf file 909 ! 910 type_ecri(1) = 't_min(X)' 911 type_ecri(2) = 't_min(X)' 912 type_ecri(3) = 't_min(X)' 913 type_ecri(4) = 't_min(X)' 914 type_ecri(5) = 't_min(X)' 915 type_ecri(6) = 't_min(X)' 916 CALL histdef2d(iff,clef_stations(iff),o_t2m_min%flag,o_t2m_min%name, "Temp 2m min", "K" ) 917 type_ecri(1) = 't_max(X)' 918 type_ecri(2) = 't_max(X)' 919 type_ecri(3) = 't_max(X)' 920 type_ecri(4) = 't_max(X)' 921 type_ecri(5) = 't_max(X)' 922 type_ecri(6) = 't_max(X)' 923 CALL histdef2d(iff,clef_stations(iff),o_t2m_max%flag,o_t2m_max%name, "Temp 2m max", "K" ) 924 ENDIF 925 type_ecri(:) = type_ecri_files(:) 926 CALL histdef2d(iff,clef_stations(iff),o_wind10m%flag,o_wind10m%name, "10-m wind speed", "m/s") 927 CALL histdef2d(iff,clef_stations(iff),o_wind10max%flag,o_wind10max%name, "10m wind speed max", "m/s") 928 CALL histdef2d(iff,clef_stations(iff),o_sicf%flag,o_sicf%name, "Sea-ice fraction", "-" ) 929 CALL histdef2d(iff,clef_stations(iff),o_q2m%flag,o_q2m%name, "Specific humidity 2m", "kg/kg") 930 CALL histdef2d(iff,clef_stations(iff),o_u10m%flag,o_u10m%name, "Vent zonal 10m", "m/s" ) 931 CALL histdef2d(iff,clef_stations(iff),o_v10m%flag,o_v10m%name, "Vent meridien 10m", "m/s") 932 CALL histdef2d(iff,clef_stations(iff),o_psol%flag,o_psol%name, "Surface Pressure", "Pa" ) 933 CALL histdef2d(iff,clef_stations(iff),o_qsurf%flag,o_qsurf%name, "Surface Air humidity", "kg/kg") 934 935 if (.not. ok_veget) then 936 CALL histdef2d(iff,clef_stations(iff),o_qsol%flag,o_qsol%name, "Soil watter content", "mm" ) 937 endif 938 939 CALL histdef2d(iff,clef_stations(iff),o_ndayrain%flag,o_ndayrain%name, "Number of dayrain(liq+sol)", "-") 940 CALL histdef2d(iff,clef_stations(iff),o_precip%flag,o_precip%name, "Precip Totale liq+sol", "kg/(s*m2)" ) 941 CALL histdef2d(iff,clef_stations(iff),o_plul%flag,o_plul%name, "Large-scale Precip.", "kg/(s*m2)") 942 CALL histdef2d(iff,clef_stations(iff),o_pluc%flag,o_pluc%name, "Convective Precip.", "kg/(s*m2)") 943 CALL histdef2d(iff,clef_stations(iff),o_snow%flag,o_snow%name, "Snow fall", "kg/(s*m2)" ) 944 CALL histdef2d(iff,clef_stations(iff),o_msnow%flag,o_msnow%name, "Surface snow amount", "kg/m2" ) 945 CALL histdef2d(iff,clef_stations(iff),o_fsnow%flag,o_fsnow%name, "Surface snow area fraction", "-" ) 946 CALL histdef2d(iff,clef_stations(iff),o_evap%flag,o_evap%name, "Evaporat", "kg/(s*m2)" ) 947 CALL histdef2d(iff,clef_stations(iff),o_tops%flag,o_tops%name, "Solar rad. at TOA", "W/m2") 948 CALL histdef2d(iff,clef_stations(iff),o_tops0%flag,o_tops0%name, "CS Solar rad. at TOA", "W/m2") 949 CALL histdef2d(iff,clef_stations(iff),o_topl%flag,o_topl%name, "IR rad. at TOA", "W/m2" ) 950 CALL histdef2d(iff,clef_stations(iff),o_topl0%flag,o_topl0%name, "IR rad. at TOA", "W/m2") 951 CALL histdef2d(iff,clef_stations(iff),o_SWupTOA%flag,o_SWupTOA%name, "SWup at TOA", "W/m2") 952 CALL histdef2d(iff,clef_stations(iff),o_SWupTOAclr%flag,o_SWupTOAclr%name, "SWup clear sky at TOA", "W/m2") 953 CALL histdef2d(iff,clef_stations(iff),o_SWdnTOA%flag,o_SWdnTOA%name, "SWdn at TOA", "W/m2" ) 954 CALL histdef2d(iff,clef_stations(iff),o_SWdnTOAclr%flag,o_SWdnTOAclr%name, "SWdn clear sky at TOA", "W/m2") 955 CALL histdef2d(iff,clef_stations(iff),o_nettop%flag,o_nettop%name, "Net dn radiatif flux at TOA", "W/m2") 956 CALL histdef2d(iff,clef_stations(iff),o_SWup200%flag,o_SWup200%name, "SWup at 200mb", "W/m2" ) 957 CALL histdef2d(iff,clef_stations(iff),o_SWup200clr%flag,o_SWup200clr%name, "SWup clear sky at 200mb", "W/m2") 958 CALL histdef2d(iff,clef_stations(iff),o_SWdn200%flag,o_SWdn200%name, "SWdn at 200mb", "W/m2" ) 959 CALL histdef2d(iff,clef_stations(iff),o_SWdn200clr%flag,o_SWdn200clr%name, "SWdn clear sky at 200mb", "W/m2") 960 CALL histdef2d(iff,clef_stations(iff),o_LWup200%flag,o_LWup200%name, "LWup at 200mb", "W/m2") 961 CALL histdef2d(iff,clef_stations(iff),o_LWup200clr%flag,o_LWup200clr%name, "LWup clear sky at 200mb", "W/m2") 962 CALL histdef2d(iff,clef_stations(iff),o_LWdn200%flag,o_LWdn200%name, "LWdn at 200mb", "W/m2") 963 CALL histdef2d(iff,clef_stations(iff),o_LWdn200clr%flag,o_LWdn200clr%name, "LWdn clear sky at 200mb", "W/m2") 964 CALL histdef2d(iff,clef_stations(iff),o_sols%flag,o_sols%name, "Solar rad. at surf.", "W/m2") 965 CALL histdef2d(iff,clef_stations(iff),o_sols0%flag,o_sols0%name, "Solar rad. at surf.", "W/m2") 966 CALL histdef2d(iff,clef_stations(iff),o_soll%flag,o_soll%name, "IR rad. at surface", "W/m2") 967 CALL histdef2d(iff,clef_stations(iff),o_radsol%flag,o_radsol%name, "Rayonnement au sol", "W/m2") 968 CALL histdef2d(iff,clef_stations(iff),o_soll0%flag,o_soll0%name, "IR rad. at surface", "W/m2") 969 CALL histdef2d(iff,clef_stations(iff),o_SWupSFC%flag,o_SWupSFC%name, "SWup at surface", "W/m2") 970 CALL histdef2d(iff,clef_stations(iff),o_SWupSFCclr%flag,o_SWupSFCclr%name, "SWup clear sky at surface", "W/m2") 971 CALL histdef2d(iff,clef_stations(iff),o_SWdnSFC%flag,o_SWdnSFC%name, "SWdn at surface", "W/m2") 972 CALL histdef2d(iff,clef_stations(iff),o_SWdnSFCclr%flag,o_SWdnSFCclr%name, "SWdn clear sky at surface", "W/m2") 973 CALL histdef2d(iff,clef_stations(iff),o_LWupSFC%flag,o_LWupSFC%name, "Upwd. IR rad. at surface", "W/m2") 974 CALL histdef2d(iff,clef_stations(iff),o_LWdnSFC%flag,o_LWdnSFC%name, "Down. IR rad. at surface", "W/m2") 975 CALL histdef2d(iff,clef_stations(iff),o_LWupSFCclr%flag,o_LWupSFCclr%name, "CS Upwd. IR rad. at surface", "W/m2") 976 CALL histdef2d(iff,clef_stations(iff),o_LWdnSFCclr%flag,o_LWdnSFCclr%name, "Down. CS IR rad. at surface", "W/m2") 977 CALL histdef2d(iff,clef_stations(iff),o_bils%flag,o_bils%name, "Surf. total heat flux", "W/m2") 978 CALL histdef2d(iff,clef_stations(iff),o_sens%flag,o_sens%name, "Sensible heat flux", "W/m2") 979 CALL histdef2d(iff,clef_stations(iff),o_fder%flag,o_fder%name, "Heat flux derivation", "W/m2") 980 CALL histdef2d(iff,clef_stations(iff),o_ffonte%flag,o_ffonte%name, "Thermal flux for snow melting", "W/m2") 981 CALL histdef2d(iff,clef_stations(iff),o_fqcalving%flag,o_fqcalving%name, "Ice Calving", "kg/m2/s") 982 CALL histdef2d(iff,clef_stations(iff),o_fqfonte%flag,o_fqfonte%name, "Land ice melt", "kg/m2/s") 983 984 CALL histdef2d(iff,clef_stations(iff),o_taux%flag,o_taux%name, "Zonal wind stress","Pa") 985 CALL histdef2d(iff,clef_stations(iff),o_tauy%flag,o_tauy%name, "Meridional wind stress","Pa") 986 987 DO nsrf = 1, nbsrf 988 CALL histdef2d(iff,clef_stations(iff),o_pourc_srf(nsrf)%flag,o_pourc_srf(nsrf)%name,"% "//clnsurf(nsrf),"%") 989 CALL histdef2d(iff,clef_stations(iff),o_fract_srf(nsrf)%flag,o_fract_srf(nsrf)%name,"Fraction "//clnsurf(nsrf),"1") 990 CALL histdef2d(iff,clef_stations(iff), & 991 o_taux_srf(nsrf)%flag,o_taux_srf(nsrf)%name,"Zonal wind stress"//clnsurf(nsrf),"Pa") 992 CALL histdef2d(iff,clef_stations(iff), & 993 o_tauy_srf(nsrf)%flag,o_tauy_srf(nsrf)%name,"Meridional wind stress "//clnsurf(nsrf),"Pa") 994 CALL histdef2d(iff,clef_stations(iff), & 995 o_tsol_srf(nsrf)%flag,o_tsol_srf(nsrf)%name,"Temperature "//clnsurf(nsrf),"K") 996 CALL histdef2d(iff,clef_stations(iff), & 997 o_u10m_srf(nsrf)%flag,o_u10m_srf(nsrf)%name,"Vent Zonal 10m "//clnsurf(nsrf),"m/s") 998 CALL histdef2d(iff,clef_stations(iff), & 999 o_evap_srf(nsrf)%flag,o_evap_srf(nsrf)%name,"evaporation at surface "//clnsurf(nsrf),"kg/(s*m2)") 1000 CALL histdef2d(iff,clef_stations(iff), & 1001 o_v10m_srf(nsrf)%flag,o_v10m_srf(nsrf)%name,"Vent meredien 10m "//clnsurf(nsrf),"m/s") 1002 CALL histdef2d(iff,clef_stations(iff), & 1003 o_t2m_srf(nsrf)%flag,o_t2m_srf(nsrf)%name,"Temp 2m "//clnsurf(nsrf),"K") 1004 CALL histdef2d(iff,clef_stations(iff), & 1005 o_sens_srf(nsrf)%flag,o_sens_srf(nsrf)%name,"Sensible heat flux "//clnsurf(nsrf),"W/m2") 1006 CALL histdef2d(iff,clef_stations(iff), & 1007 o_lat_srf(nsrf)%flag,o_lat_srf(nsrf)%name,"Latent heat flux "//clnsurf(nsrf),"W/m2") 1008 CALL histdef2d(iff,clef_stations(iff), & 1009 o_flw_srf(nsrf)%flag,o_flw_srf(nsrf)%name,"LW "//clnsurf(nsrf),"W/m2") 1010 CALL histdef2d(iff,clef_stations(iff), & 1011 o_fsw_srf(nsrf)%flag,o_fsw_srf(nsrf)%name,"SW "//clnsurf(nsrf),"W/m2") 1012 CALL histdef2d(iff,clef_stations(iff), & 1013 o_wbils_srf(nsrf)%flag,o_wbils_srf(nsrf)%name,"Bilan sol "//clnsurf(nsrf),"W/m2" ) 1014 CALL histdef2d(iff,clef_stations(iff), & 1015 o_wbilo_srf(nsrf)%flag,o_wbilo_srf(nsrf)%name,"Bilan eau "//clnsurf(nsrf),"kg/(m2*s)") 1016 if (iflag_pbl>1 .and. lev_files(iff).gt.10 ) then 1017 CALL histdef2d(iff,clef_stations(iff), & 1018 o_tke_srf(nsrf)%flag,o_tke_srf(nsrf)%name,"Max Turb. Kinetic Energy "//clnsurf(nsrf),"-") 1019 1020 IF (.NOT.clef_stations(iff)) THEN 1021 ! 1022 !IM: there is no way to have one single value in a netcdf file 1023 ! 1024 type_ecri(1) = 't_max(X)' 1025 type_ecri(2) = 't_max(X)' 1026 type_ecri(3) = 't_max(X)' 1027 type_ecri(4) = 't_max(X)' 1028 type_ecri(5) = 't_max(X)' 1029 type_ecri(6) = 't_max(X)' 1030 CALL histdef2d(iff,clef_stations(iff), & 1031 o_tke_max_srf(nsrf)%flag,o_tke_max_srf(nsrf)%name,"Max Turb. Kinetic Energy "//clnsurf(nsrf),"-") 1032 type_ecri(:) = type_ecri_files(:) 1033 ENDIF 1034 1035 endif 1036 1037 CALL histdef2d(iff,clef_stations(iff), & 1038 o_albe_srf(nsrf)%flag,o_albe_srf(nsrf)%name,"Albedo VIS surf. "//clnsurf(nsrf),"-") 1039 CALL histdef2d(iff,clef_stations(iff), & 1040 o_rugs_srf(nsrf)%flag,o_rugs_srf(nsrf)%name,"Surface roughness "//clnsurf(nsrf),"m") 1041 CALL histdef2d(iff,clef_stations(iff), & 1042 o_ages_srf(nsrf)%flag,o_ages_srf(nsrf)%name,"Snow age", "day") 1043 END DO 1044 1045 IF (new_aod .AND. (.NOT. aerosol_couple)) THEN 1046 IF (ok_ade.OR.ok_aie) THEN 1047 1048 CALL histdef2d(iff,clef_stations(iff), & 1049 o_od550aer%flag,o_od550aer%name, "Total aerosol optical depth at 550nm", "-") 1050 CALL histdef2d(iff,clef_stations(iff), & 1051 o_od865aer%flag,o_od865aer%name, "Total aerosol optical depth at 870nm", "-") 1052 CALL histdef2d(iff,clef_stations(iff), & 1053 o_absvisaer%flag,o_absvisaer%name, "Absorption aerosol visible optical depth", "-") 1054 CALL histdef2d(iff,clef_stations(iff), & 1055 o_od550lt1aer%flag,o_od550lt1aer%name, "Fine mode optical depth", "-") 1056 1057 1058 CALL histdef2d(iff,clef_stations(iff), & 1059 o_sconcso4%flag,o_sconcso4%name,"Surface Concentration of Sulfate ","kg/m3") 1060 CALL histdef2d(iff,clef_stations(iff), & 1061 o_sconcoa%flag,o_sconcoa%name,"Surface Concentration of Organic Aerosol ","kg/m3") 1062 CALL histdef2d(iff,clef_stations(iff), & 1063 o_sconcbc%flag,o_sconcbc%name,"Surface Concentration of Black Carbon ","kg/m3") 1064 CALL histdef2d(iff,clef_stations(iff), & 1065 o_sconcss%flag,o_sconcss%name,"Surface Concentration of Sea Salt ","kg/m3") 1066 CALL histdef2d(iff,clef_stations(iff), & 1067 o_sconcdust%flag,o_sconcdust%name,"Surface Concentration of Dust ","kg/m3") 1068 CALL histdef3d(iff,clef_stations(iff), & 1069 o_concso4%flag,o_concso4%name,"Concentration of Sulfate ","kg/m3") 1070 CALL histdef3d(iff,clef_stations(iff), & 1071 o_concoa%flag,o_concoa%name,"Concentration of Organic Aerosol ","kg/m3") 1072 CALL histdef3d(iff,clef_stations(iff), & 1073 o_concbc%flag,o_concbc%name,"Concentration of Black Carbon ","kg/m3") 1074 CALL histdef3d(iff,clef_stations(iff), & 1075 o_concss%flag,o_concss%name,"Concentration of Sea Salt ","kg/m3") 1076 CALL histdef3d(iff,clef_stations(iff), & 1077 o_concdust%flag,o_concdust%name,"Concentration of Dust ","kg/m3") 1078 CALL histdef2d(iff,clef_stations(iff), & 1079 o_loadso4%flag,o_loadso4%name,"Column Load of Sulfate ","kg/m2") 1080 CALL histdef2d(iff,clef_stations(iff), & 1081 o_loadoa%flag,o_loadoa%name,"Column Load of Organic Aerosol ","kg/m2") 1082 CALL histdef2d(iff,clef_stations(iff), & 1083 o_loadbc%flag,o_loadbc%name,"Column Load of Black Carbon ","kg/m2") 1084 CALL histdef2d(iff,clef_stations(iff), & 1085 o_loadss%flag,o_loadss%name,"Column Load of Sea Salt ","kg/m2") 1086 CALL histdef2d(iff,clef_stations(iff), & 1087 o_loaddust%flag,o_loaddust%name,"Column Load of Dust ","kg/m2") 1088 1089 DO naero = 1, naero_spc 1090 CALL histdef2d(iff,clef_stations(iff), & 1091 o_tausumaero(naero)%flag,o_tausumaero(naero)%name,"Aerosol Optical depth at 550 nm "//name_aero(naero),"1") 1092 END DO 1093 ENDIF 1094 ENDIF 1095 1096 IF (ok_ade) THEN 1097 CALL histdef2d(iff,clef_stations(iff), & 1098 o_topswad%flag,o_topswad%name, "ADE at TOA", "W/m2") 1099 CALL histdef2d(iff,clef_stations(iff), & 1100 o_solswad%flag,o_solswad%name, "ADE at SRF", "W/m2") 1101 1102 CALL histdef2d(iff,clef_stations(iff), & 1103 o_swtoaas_nat%flag,o_swtoaas_nat%name, "Natural aerosol radiative forcing all-sky at TOA", "W/m2") 1104 CALL histdef2d(iff,clef_stations(iff), & 1105 o_swsrfas_nat%flag,o_swsrfas_nat%name, "Natural aerosol radiative forcing all-sky at SRF", "W/m2") 1106 CALL histdef2d(iff,clef_stations(iff), & 1107 o_swtoacs_nat%flag,o_swtoacs_nat%name, "Natural aerosol radiative forcing clear-sky at TOA", "W/m2") 1108 CALL histdef2d(iff,clef_stations(iff), & 1109 o_swsrfcs_nat%flag,o_swsrfcs_nat%name, "Natural aerosol radiative forcing clear-sky at SRF", "W/m2") 1110 1111 CALL histdef2d(iff,clef_stations(iff), & 1112 o_swtoaas_ant%flag,o_swtoaas_ant%name, "Anthropogenic aerosol radiative forcing all-sky at TOA", "W/m2") 1113 CALL histdef2d(iff,clef_stations(iff), & 1114 o_swsrfas_ant%flag,o_swsrfas_ant%name, "Anthropogenic aerosol radiative forcing all-sky at SRF", "W/m2") 1115 CALL histdef2d(iff,clef_stations(iff), & 1116 o_swtoacs_ant%flag,o_swtoacs_ant%name, "Anthropogenic aerosol radiative forcing clear-sky at TOA", "W/m2") 1117 CALL histdef2d(iff,clef_stations(iff), & 1118 o_swsrfcs_ant%flag,o_swsrfcs_ant%name, "Anthropogenic aerosol radiative forcing clear-sky at SRF", "W/m2") 1119 1120 IF (.NOT. aerosol_couple) THEN 1121 CALL histdef2d(iff,clef_stations(iff), & 1122 o_swtoacf_nat%flag,o_swtoacf_nat%name, "Natural aerosol impact on cloud radiative forcing at TOA", "W/m2") 1123 CALL histdef2d(iff,clef_stations(iff), & 1124 o_swsrfcf_nat%flag,o_swsrfcf_nat%name, "Natural aerosol impact on cloud radiative forcing at SRF", "W/m2") 1125 CALL histdef2d(iff, clef_stations(iff), o_swtoacf_ant%flag, & 1126 o_swtoacf_ant%name, & 1127 "Anthropogenic aerosol impact on cloud radiative forcing at TOA", & 1128 "W/m2") 1129 CALL histdef2d(iff, clef_stations(iff), o_swsrfcf_ant%flag, & 1130 o_swsrfcf_ant%name, & 1131 "Anthropogenic aerosol impact on cloud radiative forcing at SRF", & 1132 "W/m2") 1133 CALL histdef2d(iff,clef_stations(iff), & 1134 o_swtoacf_zero%flag,o_swtoacf_zero%name, "Cloud radiative forcing (allsky-clearsky fluxes) at TOA", "W/m2") 1135 CALL histdef2d(iff,clef_stations(iff), & 1136 o_swsrfcf_zero%flag,o_swsrfcf_zero%name, "Cloud radiative forcing (allsky-clearsky fluxes) at SRF", "W/m2") 1137 ENDIF 1138 ENDIF 1139 1140 IF (ok_aie) THEN 1141 CALL histdef2d(iff,clef_stations(iff), & 1142 o_topswai%flag,o_topswai%name, "AIE at TOA", "W/m2") 1143 CALL histdef2d(iff,clef_stations(iff), & 1144 o_solswai%flag,o_solswai%name, "AIE at SFR", "W/m2") 1145 !Cloud droplet number concentration 1146 CALL histdef3d(iff,clef_stations(iff), & 1147 o_scdnc%flag,o_scdnc%name, "Cloud droplet number concentration","m-3") 1148 CALL histdef2d(iff,clef_stations(iff), & 1149 o_cldncl%flag,o_cldncl%name, "CDNC at top of liquid water cloud", "m-3") 1150 CALL histdef3d(iff,clef_stations(iff), & 1151 o_reffclws%flag,o_reffclws%name, "Stratiform Cloud Droplet Effective Radius (aerosol diags.)","m") 1152 CALL histdef3d(iff,clef_stations(iff), & 1153 o_reffclwc%flag,o_reffclwc%name, "Convective Cloud Droplet Effective Radius (aerosol diags.)","m") 1154 CALL histdef2d(iff,clef_stations(iff), & 1155 o_cldnvi%flag,o_cldnvi%name, "Column Integrated Cloud Droplet Number", "m-2") 1156 CALL histdef3d(iff,clef_stations(iff), & 1157 o_lcc3d%flag,o_lcc3d%name, "Cloud liquid fraction","1") 1158 CALL histdef3d(iff,clef_stations(iff), & 1159 o_lcc3dcon%flag,o_lcc3dcon%name, "Convective cloud liquid fraction","1") 1160 CALL histdef3d(iff,clef_stations(iff), & 1161 o_lcc3dstra%flag,o_lcc3dstra%name, "Stratiform cloud liquid fraction","1") 1162 CALL histdef2d(iff,clef_stations(iff), & 1163 o_lcc%flag,o_lcc%name, "Cloud liquid fraction at top of cloud","1") 1164 CALL histdef2d(iff,clef_stations(iff), & 1165 o_reffclwtop%flag,o_reffclwtop%name, "Droplet effective radius at top of liquid water cloud", "m") 1166 ENDIF 1167 1168 1169 CALL histdef2d(iff,clef_stations(iff), & 1170 o_alb1%flag,o_alb1%name, "Surface VIS albedo", "-") 1171 CALL histdef2d(iff,clef_stations(iff), & 1172 o_alb2%flag,o_alb2%name, "Surface Near IR albedo", "-") 1173 CALL histdef2d(iff,clef_stations(iff), & 1174 o_cdrm%flag,o_cdrm%name, "Momentum drag coef.", "-") 1175 CALL histdef2d(iff,clef_stations(iff), & 1176 o_cdrh%flag,o_cdrh%name, "Heat drag coef.", "-" ) 1177 CALL histdef2d(iff,clef_stations(iff), & 1178 o_cldl%flag,o_cldl%name, "Low-level cloudiness", "-") 1179 CALL histdef2d(iff,clef_stations(iff), & 1180 o_cldm%flag,o_cldm%name, "Mid-level cloudiness", "-") 1181 CALL histdef2d(iff,clef_stations(iff), & 1182 o_cldh%flag,o_cldh%name, "High-level cloudiness", "-") 1183 CALL histdef2d(iff,clef_stations(iff), & 1184 o_cldt%flag,o_cldt%name, "Total cloudiness", "-") 1185 CALL histdef2d(iff,clef_stations(iff), & 1186 o_cldq%flag,o_cldq%name, "Cloud liquid water path", "kg/m2") 1187 CALL histdef2d(iff,clef_stations(iff), & 1188 o_lwp%flag,o_lwp%name, "Cloud water path", "kg/m2") 1189 CALL histdef2d(iff,clef_stations(iff), & 1190 o_iwp%flag,o_iwp%name, "Cloud ice water path", "kg/m2" ) 1191 CALL histdef2d(iff,clef_stations(iff), & 1192 o_ue%flag,o_ue%name, "Zonal energy transport", "-") 1193 CALL histdef2d(iff,clef_stations(iff), & 1194 o_ve%flag,o_ve%name, "Merid energy transport", "-") 1195 CALL histdef2d(iff,clef_stations(iff), & 1196 o_uq%flag,o_uq%name, "Zonal humidity transport", "-") 1197 CALL histdef2d(iff,clef_stations(iff), & 1198 o_vq%flag,o_vq%name, "Merid humidity transport", "-") 1199 1200 IF(iflag_con.GE.3) THEN ! sb 1201 CALL histdef2d(iff,clef_stations(iff), & 1202 o_cape%flag,o_cape%name, "Conv avlbl pot ener", "J/kg") 1203 CALL histdef2d(iff,clef_stations(iff), & 1204 o_pbase%flag,o_pbase%name, "Cld base pressure", "Pa") 1205 CALL histdef2d(iff,clef_stations(iff), & 1206 o_ptop%flag,o_ptop%name, "Cld top pressure", "Pa") 1207 CALL histdef2d(iff,clef_stations(iff), & 1208 o_fbase%flag,o_fbase%name, "Cld base mass flux", "kg/m2/s") 1209 if (iflag_con /= 30) then 1210 CALL histdef2d(iff,clef_stations(iff), & 1211 o_plcl%flag,o_plcl%name, "Lifting Condensation Level", "hPa") 1212 CALL histdef2d(iff,clef_stations(iff), & 1213 o_plfc%flag,o_plfc%name, "Level of Free Convection", "hPa") 1214 CALL histdef2d(iff,clef_stations(iff), & 1215 o_wbeff%flag,o_wbeff%name, "Conv. updraft velocity at LFC", "m/s") 1216 end if 1217 CALL histdef2d(iff,clef_stations(iff), & 1218 o_prw%flag,o_prw%name, "Precipitable water", "kg/m2") 1219 IF (.NOT.clef_stations(iff)) THEN 1220 ! 1221 !IM: there is no way to have one single value in a netcdf file 1222 ! 1223 type_ecri(1) = 't_max(X)' 1224 type_ecri(2) = 't_max(X)' 1225 type_ecri(3) = 't_max(X)' 1226 type_ecri(4) = 't_max(X)' 1227 type_ecri(5) = 't_max(X)' 1228 type_ecri(6) = 't_max(X)' 1229 CALL histdef2d(iff,clef_stations(iff), & 1230 o_cape_max%flag,o_cape_max%name, "CAPE max.", "J/kg") 1231 ENDIF 1232 type_ecri(:) = type_ecri_files(:) 1233 CALL histdef3d(iff,clef_stations(iff), & 1234 o_upwd%flag,o_upwd%name, "saturated updraft", "kg/m2/s") 1235 CALL histdef3d(iff,clef_stations(iff), & 1236 o_Ma%flag,o_Ma%name, "undilute adiab updraft", "kg/m2/s") 1237 CALL histdef3d(iff,clef_stations(iff), & 1238 o_dnwd%flag,o_dnwd%name, "saturated downdraft", "kg/m2/s") 1239 CALL histdef3d(iff,clef_stations(iff), & 1240 o_dnwd0%flag,o_dnwd0%name, "unsat. downdraft", "kg/m2/s") 1241 CALL histdef3d(iff,clef_stations(iff), & 1242 o_mc%flag,o_mc%name, "Convective mass flux", "kg/m2/s") 1243 type_ecri(1) = 'inst(X)' 1244 type_ecri(2) = 'inst(X)' 1245 type_ecri(3) = 'inst(X)' 1246 type_ecri(4) = 'inst(X)' 1247 type_ecri(5) = 'inst(X)' 1248 type_ecri(6) = 'inst(X)' 1249 CALL histdef2d(iff,clef_stations(iff), & 1250 o_ftime_con%flag,o_ftime_con%name, "Fraction of time convection Occurs", " ") 1251 type_ecri(:) = type_ecri_files(:) 1252 ENDIF !iflag_con .GE. 3 1253 1254 CALL histdef2d(iff,clef_stations(iff), & 1255 o_s_pblh%flag,o_s_pblh%name, "Boundary Layer Height", "m") 1256 CALL histdef2d(iff,clef_stations(iff), & 1257 o_s_pblt%flag,o_s_pblt%name, "t at Boundary Layer Height", "K") 1258 CALL histdef2d(iff,clef_stations(iff), & 1259 o_s_lcl%flag,o_s_lcl%name, "Condensation level", "m") 1260 CALL histdef2d(iff,clef_stations(iff), & 1261 o_s_therm%flag,o_s_therm%name, "Exces du thermique", "K") 1262 !IM : Les champs suivants (s_oliqCL, s_cteiCL, s_trmb1, s_trmb2, s_trmb3) ne sont pas definis dans HBTM.F 1263 !CALL histdef2d(iff,clef_stations(iff), & 1264 !o_s_capCL%flag,o_s_capCL%name, "Conv avlbl pot enerfor ABL", "J/m2" ) 1265 !CALL histdef2d(iff,clef_stations(iff), & 1266 !o_s_oliqCL%flag,o_s_oliqCL%name, "Liq Water in BL", "kg/m2") 1267 !CALL histdef2d(iff,clef_stations(iff), & 1268 !o_s_cteiCL%flag,o_s_cteiCL%name, "Instability criteria(ABL)", "K") 1269 !CALL histdef2d(iff,clef_stations(iff), & 1270 !o_s_trmb1%flag,o_s_trmb1%name, "deep_cape(HBTM2)", "J/m2") 1271 !CALL histdef2d(iff,clef_stations(iff), & 1272 !o_s_trmb2%flag,o_s_trmb2%name, "inhibition (HBTM2)", "J/m2") 1273 !CALL histdef2d(iff,clef_stations(iff), & 1274 !o_s_trmb3%flag,o_s_trmb3%name, "Point Omega (HBTM2)", "m") 1275 1276 ! Champs interpolles sur des niveaux de pression 1277 1278 type_ecri(1) = 'inst(X)' 1279 type_ecri(2) = 'inst(X)' 1280 type_ecri(3) = 'inst(X)' 1281 type_ecri(4) = 'inst(X)' 1282 type_ecri(5) = 'inst(X)' 1283 type_ecri(6) = 'inst(X)' 1284 1285 ! Attention a reverifier 1286 1287 ilev=0 1288 DO k=1, nlevSTD 1289 bb2=clevSTD(k) 1290 IF(bb2.EQ."850".OR.bb2.EQ."700".OR.bb2.EQ."500".OR.bb2.EQ."200" & 1291 .OR.bb2.EQ."100".OR.bb2.EQ."50".OR.bb2.EQ."10")THEN 1292 ilev=ilev+1 1293 ! print*,'ilev k bb2 flag name ',ilev,k, bb2,o_uSTDlevs(ilev)%flag,o_uSTDlevs(ilev)%name 1294 CALL histdef2d(iff,clef_stations(iff), & 1295 o_uSTDlevs(ilev)%flag,o_uSTDlevs(ilev)%name,"Zonal wind "//bb2//"hPa", "m/s") 1296 CALL histdef2d(iff,clef_stations(iff), & 1297 o_vSTDlevs(ilev)%flag,o_vSTDlevs(ilev)%name,"Meridional wind "//bb2//"hPa", "m/s") 1298 CALL histdef2d(iff,clef_stations(iff), & 1299 o_wSTDlevs(ilev)%flag,o_wSTDlevs(ilev)%name,"Vertical wind "//bb2//"hPa", "Pa/s") 1300 CALL histdef2d(iff,clef_stations(iff), & 1301 o_zSTDlevs(ilev)%flag,o_zSTDlevs(ilev)%name,"Geopotential height "//bb2//"hPa", "m") 1302 CALL histdef2d(iff,clef_stations(iff), & 1303 o_qSTDlevs(ilev)%flag,o_qSTDlevs(ilev)%name,"Specific humidity "//bb2//"hPa", "kg/kg" ) 1304 CALL histdef2d(iff,clef_stations(iff), & 1305 o_tSTDlevs(ilev)%flag,o_tSTDlevs(ilev)%name,"Temperature "//bb2//"hPa", "K") 1306 ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR."500".OR.bb2.EQ."200".OR.bb2.EQ."50".OR.bb2.EQ."10") 1307 ENDDO 1308 type_ecri(:) = type_ecri_files(:) 1309 1310 CALL histdef2d(iff,clef_stations(iff), & 1311 o_t_oce_sic%flag,o_t_oce_sic%name, "Temp mixte oce-sic", "K") 1312 1313 IF (type_ocean=='slab') & 1314 CALL histdef2d(iff,clef_stations(iff), & 1315 o_slab_bils%flag, o_slab_bils%name,"Bilan au sol sur ocean slab", "W/m2") 1316 1317 ! Couplage conv-CL 1318 IF (iflag_con.GE.3) THEN 1319 IF (iflag_coupl>=1) THEN 1320 CALL histdef2d(iff,clef_stations(iff), & 1321 o_ale_bl%flag,o_ale_bl%name, "ALE BL", "m2/s2") 1322 CALL histdef2d(iff,clef_stations(iff), & 1323 o_alp_bl%flag,o_alp_bl%name, "ALP BL", "m2/s2") 1324 ENDIF 1325 ENDIF !(iflag_con.GE.3) 1326 1327 CALL histdef2d(iff,clef_stations(iff), & 1328 o_weakinv%flag,o_weakinv%name, "Weak inversion", "-") 1329 CALL histdef2d(iff,clef_stations(iff), & 1330 o_dthmin%flag,o_dthmin%name, "dTheta mini", "K/m") 1331 CALL histdef2d(iff,clef_stations(iff), & 1332 o_rh2m%flag,o_rh2m%name, "Relative humidity at 2m", "%" ) 1333 1334 IF (.NOT.clef_stations(iff)) THEN 1335 ! 1336 !IM: there is no way to have one single value in a netcdf file 1337 ! 1338 type_ecri(1) = 't_min(X)' 1339 type_ecri(2) = 't_min(X)' 1340 type_ecri(3) = 't_min(X)' 1341 type_ecri(4) = 't_min(X)' 1342 type_ecri(5) = 't_min(X)' 1343 type_ecri(6) = 't_min(X)' 1344 CALL histdef2d(iff,clef_stations(iff),o_rh2m_min%flag,o_rh2m_min%name, "Min Relative humidity at 2m", "%" ) 1345 type_ecri(1) = 't_max(X)' 1346 type_ecri(2) = 't_max(X)' 1347 type_ecri(3) = 't_max(X)' 1348 type_ecri(4) = 't_max(X)' 1349 type_ecri(5) = 't_max(X)' 1350 type_ecri(6) = 't_max(X)' 1351 CALL histdef2d(iff,clef_stations(iff),o_rh2m_max%flag,o_rh2m_max%name, "Max Relative humidity at 2m", "%" ) 1352 ENDIF 1353 1354 type_ecri(:) = type_ecri_files(:) 1355 CALL histdef2d(iff,clef_stations(iff),o_qsat2m%flag,o_qsat2m%name, "Saturant humidity at 2m", "%") 1356 CALL histdef2d(iff,clef_stations(iff),o_tpot%flag,o_tpot%name, "Surface air potential temperature", "K") 1357 CALL histdef2d(iff,clef_stations(iff), & 1358 o_tpote%flag,o_tpote%name, "Surface air equivalent potential temperature", "K") 1359 CALL histdef2d(iff,clef_stations(iff),o_SWnetOR%flag,o_SWnetOR%name, "Sfce net SW radiation OR", "W/m2") 1360 CALL histdef2d(iff,clef_stations(iff),o_SWdownOR%flag,o_SWdownOR%name, "Sfce incident SW radiation OR", "W/m2") 1361 CALL histdef2d(iff,clef_stations(iff),o_LWdownOR%flag,o_LWdownOR%name, "Sfce incident LW radiation OR", "W/m2") 1362 CALL histdef2d(iff,clef_stations(iff),o_snowl%flag,o_snowl%name, "Solid Large-scale Precip.", "kg/(m2*s)") 1363 1364 CALL histdef2d(iff,clef_stations(iff),o_solldown%flag,o_solldown%name, "Down. IR rad. at surface", "W/m2") 1365 CALL histdef2d(iff,clef_stations(iff),o_dtsvdfo%flag,o_dtsvdfo%name, "Boundary-layer dTs(o)", "K/s") 1366 CALL histdef2d(iff,clef_stations(iff),o_dtsvdft%flag,o_dtsvdft%name, "Boundary-layer dTs(t)", "K/s") 1367 CALL histdef2d(iff,clef_stations(iff),o_dtsvdfg%flag,o_dtsvdfg%name, "Boundary-layer dTs(g)", "K/s") 1368 CALL histdef2d(iff,clef_stations(iff),o_dtsvdfi%flag,o_dtsvdfi%name, "Boundary-layer dTs(g)", "K/s") 1369 CALL histdef2d(iff,clef_stations(iff),o_rugs%flag,o_rugs%name, "rugosity", "-" ) 1370 1371 ! Champs 3D: 1372 CALL histdef3d(iff,clef_stations(iff),o_ec550aer%flag,o_ec550aer%name, "Extinction at 550nm", "m^-1") 1373 CALL histdef3d(iff,clef_stations(iff),o_lwcon%flag,o_lwcon%name, "Cloud liquid water content", "kg/kg") 1374 CALL histdef3d(iff,clef_stations(iff),o_iwcon%flag,o_iwcon%name, "Cloud ice water content", "kg/kg") 1375 CALL histdef3d(iff,clef_stations(iff),o_temp%flag,o_temp%name, "Air temperature", "K" ) 1376 CALL histdef3d(iff,clef_stations(iff),o_theta%flag,o_theta%name, "Potential air temperature", "K" ) 1377 CALL histdef3d(iff,clef_stations(iff),o_ovap%flag,o_ovap%name, "Specific humidity", "kg/kg" ) 1378 CALL histdef3d(iff,clef_stations(iff), & 1379 o_ovapinit%flag,o_ovapinit%name, "Specific humidity (begin of timestep)", "kg/kg" ) 1380 CALL histdef3d(iff,clef_stations(iff), & 1381 o_geop%flag,o_geop%name, "Geopotential height", "m2/s2") 1382 CALL histdef3d(iff,clef_stations(iff), & 1383 o_vitu%flag,o_vitu%name, "Zonal wind", "m/s" ) 1384 CALL histdef3d(iff,clef_stations(iff), & 1385 o_vitv%flag,o_vitv%name, "Meridional wind", "m/s" ) 1386 CALL histdef3d(iff,clef_stations(iff), & 1387 o_vitw%flag,o_vitw%name, "Vertical wind", "Pa/s" ) 1388 CALL histdef3d(iff,clef_stations(iff), & 1389 o_pres%flag,o_pres%name, "Air pressure", "Pa" ) 1390 CALL histdef3d(iff,clef_stations(iff), & 1391 o_paprs%flag,o_paprs%name, "Air pressure Inter-Couches", "Pa" ) 1392 CALL histdef3d(iff,clef_stations(iff), & 1393 o_mass%flag,o_mass%name, "Masse Couches", "kg/m2" ) 1394 CALL histdef3d(iff,clef_stations(iff), & 1395 o_zfull%flag,o_zfull%name, "Altitude of full pressure levels", "m" ) 1396 CALL histdef3d(iff,clef_stations(iff), & 1397 o_zhalf%flag,o_zhalf%name, "Altitude of half pressure levels", "m" ) 1398 CALL histdef3d(iff,clef_stations(iff), & 1399 o_rneb%flag,o_rneb%name, "Cloud fraction", "-") 1400 CALL histdef3d(iff,clef_stations(iff), & 1401 o_rnebcon%flag,o_rnebcon%name, "Convective Cloud Fraction", "-") 1402 CALL histdef3d(iff,clef_stations(iff), & 1403 o_rhum%flag,o_rhum%name, "Relative humidity", "-") 1404 CALL histdef3d(iff,clef_stations(iff), & 1405 o_ozone%flag,o_ozone%name, "Ozone mole fraction", "-") 1406 if (read_climoz == 2) & 1407 CALL histdef3d(iff,clef_stations(iff), & 1408 o_ozone_light%flag,o_ozone_light%name, & 1409 "Daylight ozone mole fraction", "-") 1410 CALL histdef3d(iff,clef_stations(iff), & 1411 o_dtphy%flag,o_dtphy%name, "Physics dT", "K/s") 1412 CALL histdef3d(iff,clef_stations(iff), & 1413 o_dqphy%flag,o_dqphy%name, "Physics dQ", "(kg/kg)/s") 1414 CALL histdef3d(iff,clef_stations(iff), & 1415 o_cldtau%flag,o_cldtau%name, "Cloud optical thickness", "1") 1416 CALL histdef3d(iff,clef_stations(iff), & 1417 o_cldemi%flag,o_cldemi%name, "Cloud optical emissivity", "1") 1418 !IM: bug ?? dimensionnement variables (klon,klev+1) pmflxr, pmflxs, prfl, psfl 1419 CALL histdef3d(iff,clef_stations(iff), & 1420 o_pr_con_l%flag,o_pr_con_l%name, "Convective precipitation lic", " ") 1421 CALL histdef3d(iff,clef_stations(iff), & 1422 o_pr_con_i%flag,o_pr_con_i%name, "Convective precipitation ice", " ") 1423 CALL histdef3d(iff,clef_stations(iff), & 1424 o_pr_lsc_l%flag,o_pr_lsc_l%name, "Large scale precipitation lic", " ") 1425 CALL histdef3d(iff,clef_stations(iff), & 1426 o_pr_lsc_i%flag,o_pr_lsc_i%name, "Large scale precipitation ice", " ") 1427 !Cloud droplet effective radius 1428 CALL histdef3d(iff,clef_stations(iff), & 1429 o_re%flag,o_re%name, "Cloud droplet effective radius","um") 1430 CALL histdef3d(iff,clef_stations(iff), & 1431 o_fl%flag,o_fl%name, "Denominator of Cloud droplet effective radius"," ") 1432 !FH Sorties pour la couche limite 1433 if (iflag_pbl>1) then 1434 CALL histdef3d(iff,clef_stations(iff), & 1435 o_tke%flag,o_tke%name, "TKE", "m2/s2") 1436 IF (.NOT.clef_stations(iff)) THEN 1437 ! 1438 !IM: there is no way to have one single value in a netcdf file 1439 ! 1440 type_ecri(1) = 't_max(X)' 1441 type_ecri(2) = 't_max(X)' 1442 type_ecri(3) = 't_max(X)' 1443 type_ecri(4) = 't_max(X)' 1444 type_ecri(5) = 't_max(X)' 1445 type_ecri(6) = 't_max(X)' 1446 CALL histdef3d(iff,clef_stations(iff), & 1447 o_tke_max%flag,o_tke_max%name, "TKE max", "m2/s2") 1448 ENDIF 1449 type_ecri(:) = type_ecri_files(:) 1450 endif 1451 1452 CALL histdef3d(iff,clef_stations(iff), & 1453 o_kz%flag,o_kz%name, "Kz melange", "m2/s") 1454 IF (.NOT.clef_stations(iff)) THEN 1455 ! 1456 !IM: there is no way to have one single value in a netcdf file 1457 ! 1458 type_ecri(1) = 't_max(X)' 1459 type_ecri(2) = 't_max(X)' 1460 type_ecri(3) = 't_max(X)' 1461 type_ecri(4) = 't_max(X)' 1462 type_ecri(5) = 't_max(X)' 1463 type_ecri(6) = 't_max(X)' 1464 CALL histdef3d(iff,clef_stations(iff), & 1465 o_kz_max%flag,o_kz_max%name, "Kz melange max", "m2/s" ) 1466 ENDIF 1467 type_ecri(:) = type_ecri_files(:) 1468 CALL histdef3d(iff,clef_stations(iff), & 1469 o_clwcon%flag,o_clwcon%name, "Convective Cloud Liquid water content", "kg/kg") 1470 CALL histdef3d(iff,clef_stations(iff), & 1471 o_dtdyn%flag,o_dtdyn%name, "Dynamics dT", "K/s") 1472 CALL histdef3d(iff,clef_stations(iff), & 1473 o_dqdyn%flag,o_dqdyn%name, "Dynamics dQ", "(kg/kg)/s") 1474 CALL histdef3d(iff,clef_stations(iff), & 1475 o_dudyn%flag,o_dudyn%name, "Dynamics dU", "m/s2") 1476 CALL histdef3d(iff,clef_stations(iff), & 1477 o_dvdyn%flag,o_dvdyn%name, "Dynamics dV", "m/s2") 1478 CALL histdef3d(iff,clef_stations(iff), & 1479 o_dtcon%flag,o_dtcon%name, "Convection dT", "K/s") 1480 CALL histdef3d(iff,clef_stations(iff), & 1481 o_ducon%flag,o_ducon%name, "Convection du", "m/s2") 1482 CALL histdef3d(iff,clef_stations(iff), & 1483 o_dqcon%flag,o_dqcon%name, "Convection dQ", "(kg/kg)/s") 1484 1485 ! Wakes 1486 IF(iflag_con.EQ.3) THEN 1487 IF (iflag_wake >= 1) THEN 1488 CALL histdef2d(iff,clef_stations(iff), & 1489 o_ale_wk%flag,o_ale_wk%name, "ALE WK", "m2/s2") 1490 CALL histdef2d(iff,clef_stations(iff), & 1491 o_alp_wk%flag,o_alp_wk%name, "ALP WK", "m2/s2") 1492 CALL histdef2d(iff,clef_stations(iff), & 1493 o_ale%flag,o_ale%name, "ALE", "m2/s2") 1494 CALL histdef2d(iff,clef_stations(iff), & 1495 o_alp%flag,o_alp%name, "ALP", "W/m2") 1496 CALL histdef2d(iff,clef_stations(iff),o_cin%flag,o_cin%name, "Convective INhibition", "m2/s2") 1497 CALL histdef2d(iff,clef_stations(iff),o_wape%flag,o_WAPE%name, "WAPE", "m2/s2") 1498 CALL histdef2d(iff,clef_stations(iff),o_wake_h%flag,o_wake_h%name, "wake_h", "-") 1499 CALL histdef2d(iff,clef_stations(iff),o_wake_s%flag,o_wake_s%name, "wake_s", "-") 1500 CALL histdef3d(iff,clef_stations(iff),o_dtwak%flag,o_dtwak%name, "Wake dT", "K/s") 1501 CALL histdef3d(iff,clef_stations(iff),o_dqwak%flag,o_dqwak%name, "Wake dQ", "(kg/kg)/s") 1502 CALL histdef3d(iff,clef_stations(iff),o_wake_deltat%flag,o_wake_deltat%name, "wake_deltat", " ") 1503 CALL histdef3d(iff,clef_stations(iff),o_wake_deltaq%flag,o_wake_deltaq%name, "wake_deltaq", " ") 1504 CALL histdef3d(iff,clef_stations(iff),o_wake_omg%flag,o_wake_omg%name, "wake_omg", " ") 1505 ENDIF 1506 CALL histdef3d(iff,clef_stations(iff),o_Vprecip%flag,o_Vprecip%name, "precipitation vertical profile", "-") 1507 CALL histdef3d(iff,clef_stations(iff),o_ftd%flag,o_ftd%name, "tend temp due aux descentes precip", "-") 1508 CALL histdef3d(iff,clef_stations(iff),o_fqd%flag,o_fqd%name,"tend vap eau due aux descentes precip", "-") 1509 ENDIF !(iflag_con.EQ.3) 1510 1511 CALL histdef3d(iff,clef_stations(iff),o_dtlsc%flag,o_dtlsc%name, "Condensation dT", "K/s") 1512 CALL histdef3d(iff,clef_stations(iff),o_dtlschr%flag,o_dtlschr%name,"Large-scale condensational heating rate","K/s") 1513 CALL histdef3d(iff,clef_stations(iff),o_dqlsc%flag,o_dqlsc%name, "Condensation dQ", "(kg/kg)/s") 1514 CALL histdef3d(iff,clef_stations(iff),o_dtvdf%flag,o_dtvdf%name, "Boundary-layer dT", "K/s") 1515 CALL histdef3d(iff,clef_stations(iff),o_dqvdf%flag,o_dqvdf%name, "Boundary-layer dQ", "(kg/kg)/s") 1516 CALL histdef3d(iff,clef_stations(iff),o_dteva%flag,o_dteva%name, "Reevaporation dT", "K/s") 1517 CALL histdef3d(iff,clef_stations(iff),o_dqeva%flag,o_dqeva%name, "Reevaporation dQ", "(kg/kg)/s") 1518 CALL histdef3d(iff,clef_stations(iff),o_ptconv%flag,o_ptconv%name, "POINTS CONVECTIFS", " ") 1519 CALL histdef3d(iff,clef_stations(iff),o_ratqs%flag,o_ratqs%name, "RATQS", " ") 1520 CALL histdef3d(iff,clef_stations(iff),o_dtthe%flag,o_dtthe%name, "Thermal dT", "K/s") 1521 1522 if(iflag_thermals.gt.1) THEN 1523 CALL histdef3d(iff,clef_stations(iff),o_dqlscth%flag,o_dqlscth%name, "dQ therm.", "(kg/kg)/s") 1524 CALL histdef3d(iff,clef_stations(iff),o_dqlscst%flag,o_dqlscst%name, "dQ strat.", "(kg/kg)/s") 1525 CALL histdef3d(iff,clef_stations(iff),o_dtlscth%flag,o_dtlscth%name, "dQ therm.", "K/s") 1526 CALL histdef3d(iff,clef_stations(iff),o_dtlscst%flag,o_dtlscst%name, "dQ strat.", "K/s") 1527 CALL histdef2d(iff,clef_stations(iff),o_plulth%flag,o_plulth%name, "Rainfall therm.", "K/s") 1528 CALL histdef2d(iff,clef_stations(iff),o_plulst%flag,o_plulst%name, "Rainfall strat.", "K/s") 1529 CALL histdef2d(iff,clef_stations(iff),o_lmaxth%flag,o_lmaxth%name, "Upper level thermals", "") 1530 CALL histdef3d(iff,clef_stations(iff),o_ptconvth%flag,o_ptconvth%name, "POINTS CONVECTIFS therm.", " ") 1531 CALL histdef3d(iff,clef_stations(iff),o_f_th%flag,o_f_th%name, "Thermal plume mass flux", "kg/(m2*s)") 1532 CALL histdef3d(iff,clef_stations(iff),o_e_th%flag,o_e_th%name,"Thermal plume entrainment","K/s") 1533 CALL histdef3d(iff,clef_stations(iff),o_w_th%flag,o_w_th%name,"Thermal plume vertical velocity","m/s") 1534 CALL histdef3d(iff,clef_stations(iff), & 1535 o_lambda_th%flag,o_lambda_th%name,"Thermal plume vertical velocity","m/s") 1536 CALL histdef2d(iff,clef_stations(iff), & 1537 o_ftime_th%flag,o_ftime_th%name,"Fraction of time Shallow convection occurs"," ") 1538 CALL histdef3d(iff,clef_stations(iff), & 1539 o_q_th%flag,o_q_th%name, "Thermal plume total humidity", "kg/kg") 1540 CALL histdef3d(iff,clef_stations(iff), & 1541 o_a_th%flag,o_a_th%name, "Thermal plume fraction", "") 1542 CALL histdef3d(iff,clef_stations(iff), & 1543 o_d_th%flag,o_d_th%name, "Thermal plume detrainment", "K/s") 1544 1545 CALL histdef2d(iff,clef_stations(iff), & 1546 o_f0_th%flag,o_f0_th%name, "Thermal closure mass flux", "K/s") 1547 CALL histdef2d(iff,clef_stations(iff), & 1548 o_zmax_th%flag,o_zmax_th%name, "Thermal plume height", "K/s") 1549 CALL histdef3d(iff,clef_stations(iff), & 1550 o_dqthe%flag,o_dqthe%name, "Thermal dQ", "(kg/kg)/s") 1551 endif !iflag_thermals.gt.1 1552 CALL histdef3d(iff,clef_stations(iff), & 1553 o_dtajs%flag,o_dtajs%name, "Dry adjust. dT", "K/s") 1554 CALL histdef3d(iff,clef_stations(iff), & 1555 o_dqajs%flag,o_dqajs%name, "Dry adjust. dQ", "(kg/kg)/s") 1556 CALL histdef3d(iff,clef_stations(iff), & 1557 o_dtswr%flag,o_dtswr%name, "SW radiation dT", "K/s") 1558 CALL histdef3d(iff,clef_stations(iff), & 1559 o_dtsw0%flag,o_dtsw0%name, "CS SW radiation dT", "K/s") 1560 CALL histdef3d(iff,clef_stations(iff), & 1561 o_dtlwr%flag,o_dtlwr%name, "LW radiation dT", "K/s") 1562 CALL histdef3d(iff,clef_stations(iff), & 1563 o_dtlw0%flag,o_dtlw0%name, "CS LW radiation dT", "K/s") 1564 CALL histdef3d(iff,clef_stations(iff), & 1565 o_dtec%flag,o_dtec%name, "Cinetic dissip dT", "K/s") 1566 CALL histdef3d(iff,clef_stations(iff), & 1567 o_duvdf%flag,o_duvdf%name, "Boundary-layer dU", "m/s2") 1568 CALL histdef3d(iff,clef_stations(iff), & 1569 o_dvvdf%flag,o_dvvdf%name, "Boundary-layer dV", "m/s2") 1570 1571 IF (ok_orodr) THEN 1572 CALL histdef3d(iff,clef_stations(iff), & 1573 o_duoro%flag,o_duoro%name, "Orography dU", "m/s2") 1574 CALL histdef3d(iff,clef_stations(iff), & 1575 o_dvoro%flag,o_dvoro%name, "Orography dV", "m/s2") 1576 CALL histdef3d(iff,clef_stations(iff), & 1577 o_dtoro%flag,o_dtoro%name, "Orography dT", "K/s") 1578 ENDIF 1579 1580 IF (ok_orolf) THEN 1581 CALL histdef3d(iff,clef_stations(iff), & 1582 o_dulif%flag,o_dulif%name, "Orography dU", "m/s2") 1583 CALL histdef3d(iff,clef_stations(iff), & 1584 o_dvlif%flag,o_dvlif%name, "Orography dV", "m/s2") 1585 CALL histdef3d(iff,clef_stations(iff), & 1586 o_dtlif%flag,o_dtlif%name, "Orography dT", "K/s") 1587 ENDIF 1588 1589 IF (ok_hines) then 1590 CALL histdef3d(iff,clef_stations(iff), & 1591 o_duhin%flag,o_duhin%name, "Hines GWD dU", "m/s2") 1592 CALL histdef3d(iff,clef_stations(iff), & 1593 o_dvhin%flag,o_dvhin%name, "Hines GWD dV", "m/s2") 1594 1595 CALL histdef3d(iff,clef_stations(iff), & 1596 o_dthin%flag,o_dthin%name, "Hines GWD dT", "K/s") 1597 ENDIF 1598 1599 CALL histdef3d(iff,clef_stations(iff), & 1600 o_rsu%flag,o_rsu%name, "SW upward radiation", "W m-2") 1601 CALL histdef3d(iff,clef_stations(iff), & 1602 o_rsd%flag,o_rsd%name, "SW downward radiation", "W m-2") 1603 CALL histdef3d(iff,clef_stations(iff), & 1604 o_rlu%flag,o_rlu%name, "LW upward radiation", "W m-2") 1605 CALL histdef3d(iff,clef_stations(iff), & 1606 o_rld%flag,o_rld%name, "LW downward radiation", "W m-2") 1607 1608 CALL histdef3d(iff,clef_stations(iff), & 1609 o_rsucs%flag,o_rsucs%name, "SW CS upward radiation", "W m-2") 1610 CALL histdef3d(iff,clef_stations(iff), & 1611 o_rsdcs%flag,o_rsdcs%name, "SW CS downward radiation", "W m-2") 1612 CALL histdef3d(iff,clef_stations(iff), & 1613 o_rlucs%flag,o_rlucs%name, "LW CS upward radiation", "W m-2") 1614 CALL histdef3d(iff,clef_stations(iff), & 1615 o_rldcs%flag,o_rldcs%name, "LW CS downward radiation", "W m-2") 1616 1617 CALL histdef3d(iff,clef_stations(iff), & 1618 o_tnt%flag,o_tnt%name, "Tendency of air temperature", "K s-1") 1619 1620 CALL histdef3d(iff,clef_stations(iff), & 1621 o_tntc%flag,o_tntc%name, "Tendency of air temperature due to Moist Convection", & 1622 "K s-1") 1623 1624 CALL histdef3d(iff,clef_stations(iff), & 1625 o_tntr%flag,o_tntr%name, "Air temperature tendency due to Radiative heating", & 1626 "K s-1") 1627 1628 CALL histdef3d(iff,clef_stations(iff), & 1629 o_tntscpbl%flag,o_tntscpbl%name, "Air temperature tendency due to St cloud and precipitation and BL mixing", & 1630 "K s-1") 1631 1632 CALL histdef3d(iff,clef_stations(iff), & 1633 o_tnhus%flag,o_tnhus%name, "Tendency of specific humidity", "s-1") 1634 1635 CALL histdef3d(iff,clef_stations(iff), & 1636 o_tnhusc%flag,o_tnhusc%name, "Tendency of specific humidity due to convection", "s-1") 1637 1638 CALL histdef3d(iff,clef_stations(iff), & 1639 o_tnhusscpbl%flag,o_tnhusscpbl%name, "Tendency of Specific humidity due to ST cl, precip and BL mixing", & 1640 "s-1") 1641 1642 CALL histdef3d(iff,clef_stations(iff), & 1643 o_evu%flag,o_evu%name, "Eddy viscosity coefficient for Momentum Variables", "m2 s-1") 1644 1645 CALL histdef3d(iff,clef_stations(iff), & 1646 o_h2o%flag,o_h2o%name, "Mass Fraction of Water", "1") 1647 1648 CALL histdef3d(iff,clef_stations(iff), & 1649 o_mcd%flag,o_mcd%name, "Downdraft COnvective Mass Flux", "kg/(m2*s)") 1650 1651 CALL histdef3d(iff,clef_stations(iff), & 1652 o_dmc%flag,o_dmc%name, "Deep COnvective Mass Flux", "kg/(m2*s)") 1653 1654 CALL histdef3d(iff,clef_stations(iff), & 1655 o_ref_liq%flag,o_ref_liq%name, "Effective radius of convective cloud liquid water particle", "m") 1656 1657 CALL histdef3d(iff,clef_stations(iff), & 1658 o_ref_ice%flag,o_ref_ice%name, "Effective radius of startiform cloud ice particle", "m") 1659 1660 if (RCO2_per.NE.RCO2_act.OR.RCH4_per.NE.RCH4_act.OR. & 1661 RN2O_per.NE.RN2O_act.OR.RCFC11_per.NE.RCFC11_act.OR. & 1662 RCFC12_per.NE.RCFC12_act) THEN 1663 1664 CALL histdef2d(iff,clef_stations(iff),o_rsut4co2%flag,o_rsut4co2%name, & 1665 "TOA Out SW in 4xCO2 atmosphere", "W/m2") 1666 CALL histdef2d(iff,clef_stations(iff),o_rlut4co2%flag,o_rlut4co2%name, & 1667 "TOA Out LW in 4xCO2 atmosphere", "W/m2") 1668 CALL histdef2d(iff,clef_stations(iff),o_rsutcs4co2%flag,o_rsutcs4co2%name, & 1669 "TOA Out CS SW in 4xCO2 atmosphere", "W/m2") 1670 CALL histdef2d(iff,clef_stations(iff),o_rlutcs4co2%flag,o_rlutcs4co2%name, & 1671 "TOA Out CS LW in 4xCO2 atmosphere", "W/m2") 1672 1673 CALL histdef3d(iff,clef_stations(iff),o_rsu4co2%flag,o_rsu4co2%name, & 1674 "Upwelling SW 4xCO2 atmosphere", "W/m2") 1675 CALL histdef3d(iff,clef_stations(iff),o_rlu4co2%flag,o_rlu4co2%name, & 1676 "Upwelling LW 4xCO2 atmosphere", "W/m2") 1677 CALL histdef3d(iff,clef_stations(iff),o_rsucs4co2%flag,o_rsucs4co2%name, & 1678 "Upwelling CS SW 4xCO2 atmosphere", "W/m2") 1679 CALL histdef3d(iff,clef_stations(iff),o_rlucs4co2%flag,o_rlucs4co2%name, & 1680 "Upwelling CS LW 4xCO2 atmosphere", "W/m2") 1681 1682 CALL histdef3d(iff,clef_stations(iff),o_rsd4co2%flag,o_rsd4co2%name, & 1683 "Downwelling SW 4xCO2 atmosphere", "W/m2") 1684 CALL histdef3d(iff,clef_stations(iff),o_rld4co2%flag,o_rld4co2%name, & 1685 "Downwelling LW 4xCO2 atmosphere", "W/m2") 1686 CALL histdef3d(iff,clef_stations(iff),o_rsdcs4co2%flag,o_rsdcs4co2%name, & 1687 "Downwelling CS SW 4xCO2 atmosphere", "W/m2") 1688 CALL histdef3d(iff,clef_stations(iff),o_rldcs4co2%flag,o_rldcs4co2%name, & 1689 "Downwelling CS LW 4xCO2 atmosphere", "W/m2") 1690 1691 endif 1692 1693 1694 IF (nqtot>=3) THEN 1695 DO iq=3,nqtot 1696 iiq=niadv(iq) 1697 o_trac(iq-2) = ctrl_out((/ 4, 5, 1, 1, 1, 10 /),tname(iiq)) 1698 CALL histdef3d (iff,clef_stations(iff), & 1699 o_trac(iq-2)%flag,o_trac(iq-2)%name,'Tracer '//ttext(iiq), "-" ) 1700 ENDDO 1701 ENDIF 1702 1703 CALL histend(nid_files(iff)) 1704 1705 ndex2d = 0 1706 ndex3d = 0 1707 1708 ENDIF ! clef_files 1709 1710 ENDDO ! iff 1711 print*,'Fin phys_output_mod.F90' 1712 end subroutine phys_output_open 1713 1714 SUBROUTINE histdef2d (iff,lpoint,flag_var,nomvar,titrevar,unitvar) 1715 1716 use ioipsl 1717 USE dimphy 1718 USE mod_phys_lmdz_para 1719 USE iophy 1720 1721 IMPLICIT NONE 1722 1723 include "dimensions.h" 1724 include "temps.h" 1725 include "indicesol.h" 1726 include "clesphys.h" 1727 1728 integer :: iff 1729 logical :: lpoint 1730 integer, dimension(nfiles) :: flag_var 1731 character(len=20) :: nomvar 1732 character(len=*) :: titrevar 1733 character(len=*) :: unitvar 1734 1735 real zstophym 1736 1737 if (type_ecri(iff)=='inst(X)'.OR.type_ecri(iff)=='once') then 1738 zstophym=zoutm(iff) 1739 else 1740 zstophym=zdtime 1741 endif 1742 1743 ! Appel a la lecture des noms et niveau d'ecriture des variables dans output.def 1744 call conf_physoutputs(nomvar,flag_var) 1745 1746 if(.NOT.lpoint) THEN 1742 1747 if ( flag_var(iff)<=lev_files(iff) ) then 1743 call histdef (nid_files(iff),nomvar,titrevar,unitvar, &1748 call histdef (nid_files(iff),nomvar,titrevar,unitvar, & 1744 1749 iim,jj_nb,nhorim(iff), 1,1,1, -99, 32, & 1745 1750 type_ecri(iff), zstophym,zoutm(iff)) 1746 endif 1747 1751 endif 1752 else 1748 1753 if ( flag_var(iff)<=lev_files(iff) ) then 1749 call histdef (nid_files(iff),nomvar,titrevar,unitvar, &1754 call histdef (nid_files(iff),nomvar,titrevar,unitvar, & 1750 1755 npstn,1,nhorim(iff), 1,1,1, -99, 32, & 1751 1756 type_ecri(iff), zstophym,zoutm(iff)) 1752 endif1753 endif1754 end subroutine histdef2d1755 1756 SUBROUTINE histdef3d (iff,lpoint,flag_var,nomvar,titrevar,unitvar)1757 1758 use ioipsl1759 USE dimphy1760 USE mod_phys_lmdz_para1761 USE iophy1762 1763 IMPLICIT NONE1764 1765 include "dimensions.h"1766 include "temps.h"1767 include "indicesol.h"1768 include "clesphys.h"1769 1770 integer :: iff1771 logical :: lpoint1772 integer, dimension(nfiles) :: flag_var1773 character(len=20) :: nomvar1774 character(len=*) :: titrevar1775 character(len=*) :: unitvar1776 1777 real zstophym1778 1779 ! Appel a la lecture des noms et niveau d'ecriture des variables dans output.def1780 call conf_physoutputs(nomvar,flag_var)1781 1782 if (type_ecri(iff)=='inst(X)'.OR.type_ecri(iff)=='once') then1783 zstophym=zoutm(iff)1784 else1785 zstophym=zdtime1786 1757 endif 1787 1788 if(.NOT.lpoint) THEN 1758 endif 1759 end subroutine histdef2d 1760 1761 SUBROUTINE histdef3d (iff,lpoint,flag_var,nomvar,titrevar,unitvar) 1762 1763 use ioipsl 1764 USE dimphy 1765 USE mod_phys_lmdz_para 1766 USE iophy 1767 1768 IMPLICIT NONE 1769 1770 include "dimensions.h" 1771 include "temps.h" 1772 include "indicesol.h" 1773 include "clesphys.h" 1774 1775 integer :: iff 1776 logical :: lpoint 1777 integer, dimension(nfiles) :: flag_var 1778 character(len=20) :: nomvar 1779 character(len=*) :: titrevar 1780 character(len=*) :: unitvar 1781 1782 real zstophym 1783 1784 ! Appel a la lecture des noms et niveau d'ecriture des variables dans output.def 1785 call conf_physoutputs(nomvar,flag_var) 1786 1787 if (type_ecri(iff)=='inst(X)'.OR.type_ecri(iff)=='once') then 1788 zstophym=zoutm(iff) 1789 else 1790 zstophym=zdtime 1791 endif 1792 1793 if(.NOT.lpoint) THEN 1789 1794 if ( flag_var(iff)<=lev_files(iff) ) then 1790 1795 call histdef (nid_files(iff), nomvar, titrevar, unitvar, & … … 1793 1798 zstophym, zoutm(iff)) 1794 1799 endif 1795 1800 else 1796 1801 if ( flag_var(iff)<=lev_files(iff) ) then 1797 1802 call histdef (nid_files(iff), nomvar, titrevar, unitvar, & … … 1800 1805 type_ecri(iff), zstophym,zoutm(iff)) 1801 1806 endif 1802 1803 1804 1805 1807 endif 1808 end subroutine histdef3d 1809 1810 SUBROUTINE conf_physoutputs(nam_var,flag_var) 1806 1811 !!! Lecture des noms et niveau de sortie des variables dans output.def 1807 ! en utilisant les routines getin de IOIPSL1808 1809 1810 1811 1812 1813 1814 1815 1816 1817 1818 1819 1820 1821 1822 1823 1824 1825 1826 1827 1828 1829 1830 1831 1832 1833 1834 1835 1836 1837 1838 1839 1840 !1841 1842 1843 1844 1845 1846 1847 1848 1849 1850 1851 1852 1853 1854 1855 1856 1857 1858 1859 1860 1861 1862 1863 1864 1865 1866 1812 ! en utilisant les routines getin de IOIPSL 1813 use ioipsl 1814 1815 IMPLICIT NONE 1816 1817 include 'iniprint.h' 1818 1819 character(len=20) :: nam_var 1820 integer, dimension(nfiles) :: flag_var 1821 1822 IF(prt_level>10) WRITE(lunout,*)'Avant getin: nam_var flag_var ',nam_var,flag_var(:) 1823 call getin('flag_'//nam_var,flag_var) 1824 call getin('name_'//nam_var,nam_var) 1825 IF(prt_level>10) WRITE(lunout,*)'Apres getin: nam_var flag_var ',nam_var,flag_var(:) 1826 1827 END SUBROUTINE conf_physoutputs 1828 1829 SUBROUTINE convers_timesteps(str,dtime,timestep) 1830 1831 use ioipsl 1832 USE phys_cal_mod 1833 1834 IMPLICIT NONE 1835 1836 character(len=20) :: str 1837 character(len=10) :: type 1838 integer :: ipos,il 1839 real :: ttt,xxx,timestep,dayseconde,dtime 1840 parameter (dayseconde=86400.) 1841 include "temps.h" 1842 include "comconst.h" 1843 1844 ipos=scan(str,'0123456789.',.true.) 1845 ! 1846 il=len_trim(str) 1847 print*,ipos,il 1848 read(str(1:ipos),*) ttt 1849 print*,ttt 1850 type=str(ipos+1:il) 1851 1852 1853 if ( il == ipos ) then 1854 type='day' 1855 endif 1856 1857 if ( type == 'day'.or.type == 'days'.or.type == 'jours'.or.type == 'jour' ) timestep = ttt * dayseconde 1858 if ( type == 'mounths'.or.type == 'mth'.or.type == 'mois' ) then 1859 print*,'annee_ref,day_ref mon_len',annee_ref,day_ref,ioget_mon_len(annee_ref,day_ref) 1860 timestep = ttt * dayseconde * mth_len 1861 endif 1862 if ( type == 'hours'.or.type == 'hr'.or.type == 'heurs') timestep = ttt * dayseconde / 24. 1863 if ( type == 'mn'.or.type == 'minutes' ) timestep = ttt * 60. 1864 if ( type == 's'.or.type == 'sec'.or.type == 'secondes' ) timestep = ttt 1865 if ( type == 'TS' ) timestep = ttt * dtime 1866 1867 print*,'type = ',type 1868 print*,'nb j/h/m = ',ttt 1869 print*,'timestep(s)=',timestep 1870 1871 END SUBROUTINE convers_timesteps 1867 1872 1868 1873 END MODULE phys_output_mod -
LMDZ5/trunk/libf/phylmd/phys_output_write.h
r1539 r1562 614 614 ENDIF 615 615 616 IF (o_plcl%flag(iff)<=lev_files(iff)) THEN 617 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 618 $o_plcl%name,itau_w,plcl) 619 ENDIF 620 621 IF (o_plfc%flag(iff)<=lev_files(iff)) THEN 622 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 623 $o_plfc%name,itau_w,plfc) 624 ENDIF 625 626 IF (o_wbeff%flag(iff)<=lev_files(iff)) THEN 627 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 628 $o_wbeff%name,itau_w,wbeff) 629 ENDIF 630 616 if (iflag_con /= 30) then 617 if (o_plcl%flag(iff)<=lev_files(iff)) THEN 618 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 619 $ o_plcl%name,itau_w,plcl) 620 ENDIF 621 622 IF (o_plfc%flag(iff)<=lev_files(iff)) THEN 623 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 624 $ o_plfc%name,itau_w,plfc) 625 ENDIF 626 627 IF (o_wbeff%flag(iff)<=lev_files(iff)) THEN 628 CALL histwrite_phy(nid_files(iff),clef_stations(iff), 629 $ o_wbeff%name,itau_w,wbeff) 630 ENDIF 631 end if 631 632 632 633 IF (o_prw%flag(iff)<=lev_files(iff)) THEN
Note: See TracChangeset
for help on using the changeset viewer.
