Changeset 3572 for trunk/LMDZ.PLUTO/libf/phypluto/aeropacity.F90
- Timestamp:
- Jan 13, 2025, 2:52:33 PM (3 weeks ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.PLUTO/libf/phypluto/aeropacity.F90
r3275 r3572 6 6 7 7 Subroutine aeropacity(ngrid,nlayer,nq,pplay,pplev,pt, pq, & 8 aerosol,reffrad,nueffrad, QREFvis3d,QREFir3d,tau_col, & 9 cloudfrac,totcloudfrac,clearsky) 8 aerosol,reffrad,nueffrad, QREFvis3d,QREFir3d,tau_col) 10 9 11 10 use radinc_h, only : L_TAUMAX,naerkind 12 use aerosol_mod, only: iaero_haze, i_haze, iaero_generic, & 13 noaero 11 use aerosol_mod, only: iaero_haze, i_haze 14 12 USE tracer_h, only: noms,rho_n2,rho_ice,rho_q,mmol 15 13 use comcstfi_mod, only: g, pi, mugaz, avocado 16 14 use geometry_mod, only: latitude 17 use callkeys_mod, only: aerofixn2,kastprof,cloudlvl, & 18 pres_bottom_tropo,pres_top_tropo,obs_tau_col_tropo, & 19 pres_bottom_strato,pres_top_strato,obs_tau_col_strato, & 20 nlayaero, aeronlay_tauref, aeronlay_choice, & 21 aeronlay_pbot, aeronlay_ptop, aeronlay_sclhght, & 22 aerogeneric 23 use generic_tracer_index_mod, only: generic_tracer_index 15 use callkeys_mod, only: kastprof 24 16 implicit none 25 17 … … 70 62 REAL,INTENT(IN) :: QREFir3d(ngrid,nlayer,naerkind) 71 63 REAL,INTENT(OUT):: tau_col(ngrid) !column integrated visible optical depth 72 ! BENJAMIN MODIFS73 real,intent(in) :: cloudfrac(ngrid,nlayer) ! cloud fraction74 real,intent(out) :: totcloudfrac(ngrid) ! total cloud fraction75 logical,intent(in) :: clearsky76 64 77 65 real aerosol0, obs_tau_col_aurora, pm … … 96 84 integer igen_ice,igen_gas ! to store the index of generic tracer 97 85 logical dummy_bool ! dummy boolean just in case we need one 98 ! integer i_rgcs_ice(aerogeneric)99 86 ! for venus clouds 100 87 real :: p_bot,p_top,h_bot,h_top,mode_dens,h_lay … … 112 99 endif 113 100 enddo 114 115 if (noaero) then116 print*, "No active aerosols found in aeropacity"117 else118 print*, "If you would like to use aerosols, make sure any old"119 print*, "start files are updated in newstart using the option"120 print*, "q=0"121 write(*,*) "Active aerosols found in aeropacity:"122 endif123 101 124 102 firstcall=.false.
Note: See TracChangeset
for help on using the changeset viewer.