Changeset 2823 for LMDZ5/trunk/libf
- Timestamp:
- Mar 14, 2017, 11:10:55 PM (7 years ago)
- Location:
- LMDZ5/trunk/libf/phylmd
- Files:
-
- 8 edited
-
aero_mod.F90 (modified) (2 diffs)
-
phys_local_var_mod.F90 (modified) (4 diffs)
-
phys_output_ctrlout_mod.F90 (modified) (1 diff)
-
phys_output_mod.F90 (modified) (2 diffs)
-
phys_output_write_mod.F90 (modified) (4 diffs)
-
readaerosol.F90 (modified) (4 diffs)
-
readaerosol_optic.F90 (modified) (8 diffs)
-
rrtm/readaerosol_optic_rrtm.F90 (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/phylmd/aero_mod.F90
r2550 r2823 46 46 ! 2/ Total number of aerosols for which an aerosol mass is provided 47 47 48 INTEGER, PARAMETER :: naero_spc = 1 048 INTEGER, PARAMETER :: naero_spc = 13 49 49 50 50 ! Corresponding names for the aerosols … … 59 59 "CIDUSTM", & 60 60 "AIBCM ", & 61 "AIPOMM " /) 61 "AIPOMM " ,& 62 "ASNO3M ", & 63 "CSNO3M ", & 64 "CINO3M " /) 62 65 63 66 ! 3/ Number of aerosol groups -
LMDZ5/trunk/libf/phylmd/phys_local_var_mod.F90
r2806 r2823 209 209 REAL, SAVE, ALLOCATABLE :: loaddust(:) 210 210 !$OMP THREADPRIVATE(loaddust) 211 REAL, SAVE, ALLOCATABLE :: loadno3(:) 212 !$OMP THREADPRIVATE(loadno3) 211 213 REAL, SAVE, ALLOCATABLE :: load_tmp1(:) 212 214 !$OMP THREADPRIVATE(load_tmp1) … … 215 217 REAL, SAVE, ALLOCATABLE :: load_tmp3(:) 216 218 !$OMP THREADPRIVATE(load_tmp3) 217 REAL, SAVE, ALLOCATABLE :: load_tmp4(:)218 !$OMP THREADPRIVATE(load_tmp4)219 REAL, SAVE, ALLOCATABLE :: load_tmp5(:)220 !$OMP THREADPRIVATE(load_tmp5)221 REAL, SAVE, ALLOCATABLE :: load_tmp6(:)222 !$OMP THREADPRIVATE(load_tmp6)223 REAL, SAVE, ALLOCATABLE :: load_tmp7(:)224 !$OMP THREADPRIVATE(load_tmp7)225 219 226 220 !IM ajout variables CFMIP2/CMIP5 … … 605 599 allocate(loadss(klon)) 606 600 allocate(loaddust(klon)) 601 allocate(loadno3(klon)) 607 602 allocate(load_tmp1(klon)) 608 603 allocate(load_tmp2(klon)) 609 604 allocate(load_tmp3(klon)) 610 allocate(load_tmp4(klon))611 allocate(load_tmp5(klon))612 allocate(load_tmp6(klon))613 allocate(load_tmp7(klon))614 605 615 606 !IM ajout variables CFMIP2/CMIP5 … … 866 857 deallocate(loadss) 867 858 deallocate(loaddust) 859 deallocate(loadno3) 868 860 deallocate(load_tmp1) 869 861 deallocate(load_tmp2) 870 862 deallocate(load_tmp3) 871 deallocate(load_tmp4)872 deallocate(load_tmp5)873 deallocate(load_tmp6)874 deallocate(load_tmp7)875 863 deallocate(du_gwd_hines,dv_gwd_hines,d_t_hin) 876 864 deallocate(d_q_ch4) -
LMDZ5/trunk/libf/phylmd/phys_output_ctrlout_mod.F90
r2806 r2823 1136 1136 TYPE(ctrl_out), SAVE :: o_loaddust = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11, 11/), & 1137 1137 'loaddust', 'Column Load of Dust ', 'kg/m2', (/ ('', i=1, 10) /)) 1138 TYPE(ctrl_out), SAVE :: o_loadno3 = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11, 11/), & 1139 'loadno3', 'Column Load of Nitrate ', 'kg/m2', (/ ('', i=1, 10) /)) 1138 1140 TYPE(ctrl_out), SAVE :: o_swtoaas_nat = ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11, 11/), & 1139 1141 'swtoaas_nat', 'Natural aerosol radiative forcing all-sky at TOA', 'W/m2', (/ ('', i=1, 10) /)) -
LMDZ5/trunk/libf/phylmd/phys_output_mod.F90
r2752 r2823 5 5 USE indice_sol_mod 6 6 USE phys_output_var_mod 7 USE aero_mod, only : naero_spc,name_aero8 7 USE phys_output_write_mod, ONLY : phys_output_write 9 8 REAL, DIMENSION(nfiles),SAVE :: ecrit_files … … 40 39 USE phys_cal_mod, only : hour, calend 41 40 USE mod_phys_lmdz_para 42 USE aero_mod, only : naero_spc,name_aero43 41 !Martin 44 42 USE surface_data, ONLY : ok_snow -
LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90
r2816 r2823 105 105 o_concoa, o_concbc, o_concss, o_concdust, & 106 106 o_loadso4, o_loadoa, o_loadbc, o_loadss, & 107 o_loaddust, o_ tausumaero, o_tausumaero_lw, &107 o_loaddust, o_loadno3, o_tausumaero, o_tausumaero_lw, & 108 108 o_topswad, o_topswad0, o_solswad, o_solswad0, & 109 109 o_toplwad, o_toplwad0, o_sollwad, o_sollwad0, & … … 256 256 sconcoa, sconcbc, sconcss, sconcdust, concso4, concno3, & 257 257 concoa, concbc, concss, concdust, loadso4, & 258 loadoa, loadbc, loadss, loaddust, tausum_aero, &258 loadoa, loadbc, loadss, loaddust, loadno3, tausum_aero, & 259 259 topswad_aero, topswad0_aero, solswad_aero, & 260 260 solswad0_aero, topsw_aero, solsw_aero, & … … 325 325 USE geometry_mod, ONLY: cell_area 326 326 USE surface_data, ONLY: type_ocean, version_ocean, ok_veget, ok_snow 327 ! USE aero_mod, ONLY: naero_spc328 327 USE aero_mod, ONLY: naero_tot, id_STRAT_phy 329 328 USE ioipsl, ONLY: histend, histsync … … 1192 1191 CALL histwrite_phy(o_loadss, loadss) 1193 1192 CALL histwrite_phy(o_loaddust, loaddust) 1193 CALL histwrite_phy(o_loadno3, loadno3) 1194 1194 !--STRAT AER 1195 1195 ENDIF -
LMDZ5/trunk/libf/phylmd/readaerosol.F90
r2346 r2823 222 222 REAL, ALLOCATABLE, DIMENSION(:) :: varktmp 223 223 224 REAL, DIMENSION(nbp_lon,nbp_lat,12) 224 REAL, DIMENSION(nbp_lon,nbp_lat,12) :: psurf_glo2D ! Surface pression for 12 months on dynamics global grid 225 225 REAL, DIMENSION(klon_glo,12) :: psurf_glo1D ! -"- on physical global grid 226 REAL, DIMENSION(nbp_lon,nbp_lat,12) 226 REAL, DIMENSION(nbp_lon,nbp_lat,12) :: load_glo2D ! Load for 12 months on dynamics global grid 227 227 REAL, DIMENSION(klon_glo,12) :: load_glo1D ! -"- on physical global grid 228 REAL, DIMENSION(nbp_lon,nbp_lat) 229 REAL, DIMENSION(nbp_lon) 230 REAL, DIMENSION(nbp_lat) 228 REAL, DIMENSION(nbp_lon,nbp_lat) :: vartmp 229 REAL, DIMENSION(nbp_lon) :: lon_src ! longitudes in file 230 REAL, DIMENSION(nbp_lat) :: lat_src, lat_src_inv ! latitudes in file 231 231 LOGICAL :: new_file ! true if new file format detected 232 232 LOGICAL :: invert_lat ! true if the field has to be inverted for latitudes … … 245 245 246 246 WRITE(lunout,*) 'reading variable ',TRIM(varname),' in file ', TRIM(fname) 247 CALL check_err( nf90_open(TRIM(fname), NF90_NOWRITE, ncid), "pb open "//trim( varname) )247 CALL check_err( nf90_open(TRIM(fname), NF90_NOWRITE, ncid), "pb open "//trim(fname) ) 248 248 249 249 ! Test for equal longitudes and latitudes in file and model … … 338 338 !**************************************************************************************** 339 339 ! Get variable id 340 CALL check_err( nf90_inq_varid(ncid, TRIM(varname), varid),"pb inq var "//TRIM(varname) ) 341 342 ! Get the variable 343 CALL check_err( nf90_get_var(ncid, varid, varyear(:,:,:,:)),"pb get var "//TRIM(varname) ) 340 !CALL check_err( nf90_inq_varid(ncid, TRIM(varname), varid),"pb inq var "//TRIM(varname) ) 341 print *,'readaerosol ', TRIM(varname) 342 IF ( nf90_inq_varid(ncid, TRIM(varname), varid) /= NF90_NOERR ) THEN 343 ! Variable is not there 344 WRITE(lunout,*) 'Attention '//TRIM(varname)//' is not in aerosol input file' 345 varyear(:,:,:,:)=0.0 346 ELSE 347 ! Get the variable 348 CALL check_err( nf90_get_var(ncid, varid, varyear(:,:,:,:)),"pb get var "//TRIM(varname) ) 349 ENDIF 344 350 345 351 ! ++) Read surface pression, 12 month in one variable … … 353 359 !**************************************************************************************** 354 360 ! Get variable id 355 CALL check_err( nf90_inq_varid(ncid, "load_"//TRIM(varname), varid) ,"pb inq var load_"//TRIM(varname)) 356 ! Get the variable 357 CALL check_err( nf90_get_var(ncid, varid, load_glo2D),"pb get var load_"//TRIM(varname) ) 361 !CALL check_err( nf90_inq_varid(ncid, "load_"//TRIM(varname), varid) ,"pb inq var load_"//TRIM(varname)) 362 IF ( nf90_inq_varid(ncid, "load_"//TRIM(varname), varid) /= NF90_NOERR ) THEN 363 WRITE(lunout,*) 'Attention load_'//TRIM(varname)//' is not in aerosol input file' 364 load_glo2D(:,:,:)=0.0 365 ELSE 366 ! Get the variable 367 CALL check_err( nf90_get_var(ncid, varid, load_glo2D),"pb get var load_"//TRIM(varname) ) 368 ENDIF 358 369 359 370 ! ++) Read ap -
LMDZ5/trunk/libf/phylmd/readaerosol_optic.F90
r2642 r2823 16 16 USE phys_local_var_mod, only: sconcso4,sconcno3,sconcoa,sconcbc,sconcss,sconcdust, & 17 17 concso4,concno3,concoa,concbc,concss,concdust,loadso4,loadoa,loadbc,loadss,loaddust, & 18 load_tmp1,load_tmp2,load_tmp3 ,load_tmp4,load_tmp5,load_tmp6,load_tmp718 load_tmp1,load_tmp2,load_tmp3 19 19 IMPLICIT NONE 20 20 … … 93 93 94 94 ! Get bc aerosol distribution 95 CALL readaerosol_interp(id_ASBCM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, bcsol, bcsol_pi, load_tmp1 96 CALL readaerosol_interp(id_AIBCM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, bcins, bcins_pi, load_tmp2 95 CALL readaerosol_interp(id_ASBCM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, bcsol, bcsol_pi, load_tmp1) 96 CALL readaerosol_interp(id_AIBCM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, bcins, bcins_pi, load_tmp2) 97 97 loadbc(:)=load_tmp1(:)+load_tmp2(:) 98 98 ELSE … … 107 107 flag_aerosol .EQ. 6 ) THEN 108 108 109 CALL readaerosol_interp(id_ASPOMM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, pomsol, pomsol_pi, load_tmp 3)110 CALL readaerosol_interp(id_AIPOMM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, pomins, pomins_pi, load_tmp 4)111 loadoa(:)=load_tmp 3(:)+load_tmp4(:)109 CALL readaerosol_interp(id_ASPOMM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, pomsol, pomsol_pi, load_tmp1) 110 CALL readaerosol_interp(id_AIPOMM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, pomins, pomins_pi, load_tmp2) 111 loadoa(:)=load_tmp1(:)+load_tmp2(:) 112 112 ELSE 113 113 pomsol(:,:) = 0. ; pomsol_pi(:,:) = 0. … … 121 121 flag_aerosol .EQ. 6 ) THEN 122 122 123 CALL readaerosol_interp(id_SSSSM_phy ,itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, sssupco, sssupco_pi, load_tmp5)124 CALL readaerosol_interp(id_CSSSM_phy ,itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, sscoarse,sscoarse_pi, load_tmp 6)125 CALL readaerosol_interp(id_ASSSM_phy ,itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, ssacu, ssacu_pi, load_tmp7)126 loadss(:)=load_tmp5(:)+load_tmp6(:)+load_tmp7(:)127 ELSE 128 sscoarse(:,:) = 0. ; sscoarse_pi(:,:) = 0.129 ssacu(:,:) = 0. ; ssacu_pi(:,:) = 0.130 sssupco(:,:) = 0. ; sssupco_pi = 0.131 loadss=0.123 CALL readaerosol_interp(id_SSSSM_phy ,itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, sssupco, sssupco_pi, load_tmp1) 124 CALL readaerosol_interp(id_CSSSM_phy ,itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, sscoarse,sscoarse_pi, load_tmp2) 125 CALL readaerosol_interp(id_ASSSM_phy ,itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, ssacu, ssacu_pi, load_tmp3) 126 loadss(:)=load_tmp1(:)+load_tmp2(:)+load_tmp3(:) 127 ELSE 128 sscoarse(:,:) = 0. ; sscoarse_pi(:,:) = 0. 129 ssacu(:,:) = 0. ; ssacu_pi(:,:) = 0. 130 sssupco(:,:) = 0. ; sssupco_pi = 0. 131 loadss=0. 132 132 ENDIF 133 133 … … 137 137 138 138 CALL readaerosol_interp(id_CIDUSTM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, cidust, cidust_pi, loaddust) 139 140 139 ELSE 141 140 cidust(:,:) = 0. ; cidust_pi(:,:) = 0. … … 174 173 m_allaer_pi(:,:,id_CSNO3M_phy) = 0.0 175 174 m_allaer_pi(:,:,id_CINO3M_phy) = 0.0 176 177 175 ! 178 176 ! Calculate the total mass of all soluble aersosols … … 214 212 215 213 END IF 216 217 214 218 215 ! Diagnostics calculation for CMIP5 protocol … … 230 227 concdust(:,:)=m_allaer(:,:,id_CIDUSTM_phy)*1.e-9 231 228 232 233 229 END SUBROUTINE readaerosol_optic -
LMDZ5/trunk/libf/phylmd/rrtm/readaerosol_optic_rrtm.F90
r2738 r2823 17 17 USE phys_local_var_mod, only: sconcso4,sconcno3,sconcoa,sconcbc,sconcss,sconcdust, & 18 18 concso4,concno3,concoa,concbc,concss,concdust,loadso4,loadoa,loadbc,loadss,loaddust, & 19 load _tmp1,load_tmp2,load_tmp3,load_tmp4,load_tmp5,load_tmp6,load_tmp719 loadno3, load_tmp1,load_tmp2,load_tmp3 20 20 21 21 USE infotrac_phy … … 195 195 IF ( flag_aerosol .EQ. 3 .OR. flag_aerosol .EQ. 6 ) THEN 196 196 197 CALL readaerosol_interp(id_ASPOMM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, pomsol, pomsol_pi, load_tmp 3)198 CALL readaerosol_interp(id_AIPOMM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, pomins, pomins_pi, load_tmp 4)199 loadoa(:)=load_tmp 3(:)+load_tmp4(:)197 CALL readaerosol_interp(id_ASPOMM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, pomsol, pomsol_pi, load_tmp1) 198 CALL readaerosol_interp(id_AIPOMM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, pomins, pomins_pi, load_tmp2) 199 loadoa(:)=load_tmp1(:)+load_tmp2(:) 200 200 ELSE 201 201 pomsol(:,:) = 0. ; pomsol_pi(:,:) = 0. … … 208 208 209 209 CALL readaerosol_interp(id_SSSSM_phy ,itap, pdtphys,rjourvrai, & 210 debut, pplay, paprs, t_seri, sssupco, sssupco_pi, load_tmp 5)210 debut, pplay, paprs, t_seri, sssupco, sssupco_pi, load_tmp1) 211 211 CALL readaerosol_interp(id_CSSSM_phy ,itap, pdtphys,rjourvrai, & 212 debut, pplay, paprs, t_seri, sscoarse,sscoarse_pi, load_tmp 6)213 CALL readaerosol_interp(id_ASSSM_phy ,itap, pdtphys, 214 debut, pplay, paprs, t_seri, ssacu, ssacu_pi, load_tmp 7)215 loadss(:)=load_tmp 5(:)+load_tmp6(:)+load_tmp7(:)212 debut, pplay, paprs, t_seri, sscoarse,sscoarse_pi, load_tmp2) 213 CALL readaerosol_interp(id_ASSSM_phy ,itap, pdtphys,rjourvrai, & 214 debut, pplay, paprs, t_seri, ssacu, ssacu_pi, load_tmp3) 215 loadss(:)=load_tmp1(:)+load_tmp2(:)+load_tmp3(:) 216 216 ELSE 217 217 sscoarse(:,:) = 0. ; sscoarse_pi(:,:) = 0. … … 231 231 ENDIF 232 232 ! 233 ! Read and interpolate cidustm 234 IF (flag_aerosol .EQ. 6) THEN 235 236 CALL readaerosol_interp(id_ASNO3M_phy, itap, pdtphys, rjourvrai, & 237 debut, pplay, paprs, t_seri, nitracc, nitracc_pi, load_tmp1) 238 CALL readaerosol_interp(id_CSNO3M_phy, itap, pdtphys, rjourvrai, & 239 debut, pplay, paprs, t_seri, nitrcoarse, nitrcoarse_pi, load_tmp2) 240 CALL readaerosol_interp(id_CINO3M_phy, itap, pdtphys, rjourvrai, & 241 debut, pplay, paprs, t_seri, nitrinscoarse, nitrinscoarse_pi, load_tmp3) 242 loadss(:)=load_tmp1(:)+load_tmp2(:)+load_tmp3(:) 243 244 ELSE 245 nitracc(:,:) = 0.0 ; nitracc_pi(:,:) = 0.0 246 nitrcoarse(:,:) = 0.0 ; nitrcoarse_pi(:,:) = 0.0 247 nitrinscoarse(:,:) = 0.0 ; nitrinscoarse_pi(:,:)= 0.0 248 loadno3(:)=0.0 249 ENDIF 250 ! 251 ! CSSO4M is set to 0 as not reliable 233 252 sulfcoarse(:,:) = 0.0 ! CSSO4M (=SO4) + CSMSAM (=MSA) 234 253 sulfcoarse_pi(:,:) = 0.0 ! CSSO4M (=SO4) + CSMSAM (=MSA) pre-ind 235 !236 !--placeholder for offline nitrate237 !238 nitracc(:,:) = 0.0239 nitracc_pi(:,:) = 0.0240 nitrcoarse(:,:) = 0.0241 nitrcoarse_pi(:,:) = 0.0242 nitrinscoarse(:,:) = 0.0243 nitrinscoarse_pi(:,:)= 0.0244 254 245 255 ENDIF !--not aerosol_couple
Note: See TracChangeset
for help on using the changeset viewer.
