Changeset 3195 for trunk/LMDZ.PLUTO/libf/phypluto/aerosol_mod.F90

Timestamp:

Jan 31, 2024, 4:36:51 PM (10 months ago)

Author:

afalco

Message:

Pluto PCM:
Imported condense n2 from pluto.old.
Aerosol data from Pluto.old not yet working.
AF

File:

• trunk/LMDZ.PLUTO/libf/phypluto/aerosol_mod.F90 (modified) (2 diffs)

trunk/LMDZ.PLUTO/libf/phypluto/aerosol_mod.F90

@@ -8 +8 @@
 !                 corresponding aerosol was not activated in callphys.def
 !                 -- otherwise a value is set via iniaerosol
-      integer, save, protected :: iaero_n2 = 0 
-      integer, save, protected :: iaero_dust = 0
-      integer, save, protected :: iaero_h2so4 = 0
+      integer, save :: iaero_haze = 0
+      integer, save :: i_haze = 0
       logical, save, protected :: noaero = .false.
-!$OMP THREADPRIVATE(iaero_n2,iaero_h2o,iaero_dust,iaero_h2so4,noaero)
+!$OMP THREADPRIVATE(iaero_haze,i_haze,noaero)
 
 ! two-layer simple aerosol model
       integer, save, protected :: iaero_back2lay = 0
- ! NH3 cloud
-      integer, save, protected :: iaero_nh3 = 0
 ! N-layer aerosol model (replaces the 2-layer and hard-coded clouds)
       integer,dimension(:), allocatable, save, protected :: iaero_nlay
-! Auroral aerosols
-      integer, save, protected :: iaero_aurora = 0
-!$OMP THREADPRIVATE(iaero_back2lay,iaero_nh3,iaero_nlay,iaero_aurora)
+!$OMP THREADPRIVATE(iaero_back2lay,iaero_nlay)
 
 ! Generic aerosols
@@ -33 +28 @@
 contains
 
-  SUBROUTINE iniaerosol
+  !==================================================================
+   subroutine haze_prof(ngrid,nlayer,zzlay,pplay,pt,reffrad,profmmr)
+!==================================================================
+!     Purpose
+!     -------
+!     Get fixed haze properties
+!     profile of haze (from txt file) and fixed radius profile
+!
+!==================================================================
+      use radinc_h, only: naerkind
+      use datafile_mod
+      use tracer_h
+      use comcstfi_mod, only: r, pi
 
-  use mod_phys_lmdz_para, only : is_master
-  use radinc_h, only: naerkind
-  use tracer_h, only: n_rgcs, nqtot, is_rgcs
-  use callkeys_mod, only: aeron2, dusttau, &
-                          aeroback2lay, aeronh3, nlayaero, aeronlay, &
-                          aeroaurora, aerogeneric
+!-----------------------------------------------------------------------
+!     Arguments
+      Implicit none
 
-  IMPLICIT NONE
-!=======================================================================
-!   subject:
-!   --------
-!   Initialization related to aerosols 
-!   (N2 aerosols, dust, water, chemical species, ice...)   
-!
-!   author: Laura Kerber, S. Guerlet
-!   ------
-!        
-!=======================================================================
+      integer,intent(in) :: ngrid
+      integer,intent(in) :: nlayer
+      real,intent(in) :: zzlay(ngrid,nlayer)
+      real,intent(in) :: pplay(ngrid,nlayer)
+      real,intent(in) :: pt(ngrid,nlayer)
+      real, intent(in) :: reffrad(ngrid,nlayer,naerkind)      ! haze particles radii (m)
 
-  integer :: i, ia, iq
+      real, intent(out) :: profmmr(ngrid,nlayer)              ! mmr haze kg/kg
 
-  ! Special case, dyn. allocation for n-layer depending on callphys.def
-  IF(.NOT.ALLOCATED(iaero_nlay)) ALLOCATE(iaero_nlay(nlayaero))
-  iaero_nlay(:) = 0
-  ! Do the same for iaero_generic and i_rgcs_ice
-  IF (.not. allocated(iaero_generic)) allocate(iaero_generic(aerogeneric))
-  if (.not. allocated(i_rgcs_ice)) allocate(i_rgcs_ice(aerogeneric)) 
+!     Local variables
+      integer :: iaer,l,ig,ifine
 
-  ! Init of i_rgcs_ice
-  i_rgcs_ice(:) =0
-  ia = 1
-  do iq=1,nqtot
-    if (is_rgcs(iq) .eq. 1) then 
-        i_rgcs_ice(ia)=iq
-        ia = ia+1
-     endif
-  enddo
+      LOGICAL firstcall
+      SAVE firstcall
+      DATA firstcall/.true./
 
-  iaero_generic(:)=0
-  ia=0
-  if (aeron2) then
-     ia=ia+1
-     iaero_n2=ia
-  endif
-  if (is_master) write(*,*) '--- N2 aerosol = ', iaero_n2
+      !!read altitudes and haze mmrs
+      integer Nfine
+      !parameter(Nfine=21)
+      parameter(Nfine=701)
+      character(len=100) :: file_path
+      real,save :: levdat(Nfine),densdat(Nfine)
 
-  if (dusttau.gt.0) then
-     ia=ia+1
-     iaero_dust=ia
-  endif
-  if (is_master) write(*,*) '--- Dust aerosol = ', iaero_dust
+!---------------- INPUT ------------------------------------------------
 
-     
-  if (aeroback2lay) then
-     ia=ia+1
-     iaero_back2lay=ia
-  endif
-  if (is_master) write(*,*) '--- Two-layer aerosol = ', iaero_back2lay
+      !! Read data
+      IF (firstcall) then
+        firstcall=.false.
+        file_path=trim(datadir)//'/haze_prop/hazemmr.txt'
+        open(224,file=file_path,form='formatted')
+        do ifine=1,Nfine
+           read(224,*) levdat(ifine), densdat(ifine)
+        enddo
+        close(224)
+        print*, 'TB22 read Haze MMR profile'
+      ENDIF
 
-  if (aeronh3) then
-     ia=ia+1
-     iaero_nh3=ia
-  endif
-  if (is_master) write(*,*) '--- NH3 Cloud = ', iaero_nh3
+      !! Interpolate on the model vertical grid
+      do ig=1,ngrid
+        CALL interp_line(levdat,densdat,Nfine,zzlay(ig,:)/1000.,profmmr(ig,:),nlayer)
+      enddo
 
-  if (aeronlay) then
-     do i=1,nlayaero
-       ia=ia+1
-       iaero_nlay(i)=ia
-     enddo
-  endif
-  if (is_master) write(*,*) '--- N-layer aerosol = ', iaero_nlay
+      !! Get profile Mass mixing ratio from number density:    part.cm-3 --> m-3 --> m3 m-3
+      !                                --> kg m-3 --> kg/kg
+      do iaer=1,naerkind
+            if(iaer.eq.iaero_haze.and.1.eq.2) then !TB22 activate/deactivate mmr or part density
+              !print*, 'Haze profile is fixed'
+              do ig=1,ngrid
+                 do l=1,nlayer
+                    !from number density in cm-3
+                    profmmr(ig,l)=profmmr(ig,l)*1.e6*4./3.*pi*reffrad(ig,l,iaer)**3*rho_q(i_haze)/(pplay(ig,l)/(r*pt(ig,l)))
+                 enddo
+              enddo
+            endif
+      enddo
+   end subroutine haze_prof
+!==================================================================
 
-  if (aeroaurora) then
-     ia=ia+1
-     iaero_aurora=ia
-  endif
-  if (is_master) write(*,*) '--- Auroral aerosols = ', iaero_aurora
-
-  if (aerogeneric .ne. 0) then 
-     do i=1,aerogeneric
-        ia = ia+1
-        iaero_generic(i) = ia
-     enddo
-  endif 
-      
-  if (is_master) then
-    write(*,*)'--- Radiative Generic Condensable Species = ',iaero_generic
-
-    write(*,*) '=== Number of aerosols= ', ia
-  endif ! of is_master
-
-! For the zero aerosol case, we currently make a dummy n2 aerosol which is zero everywhere.
-! (See aeropacity.F90 for how this works). A better solution would be to turn off the 
-! aerosol machinery in the no aerosol case, but this would be complicated. LK
-
-  if (ia.eq.0) then  !For the zero aerosol case. 
-     ia = 0
-     noaero = .true.
-  endif
-
-  if (.not.noaero .and. ia.ne.naerkind) then
-    if (is_master) then
-      print*, 'Aerosols counted not equal to naerkind'
-      print*, 'set correct value for nearkind in callphys.def'
-      print*, 'which should be ',ia
-      print*, 'according to current options in callphys.def'
-      print*, 'or change/correct incompatible options there'
-      print*, 'Abort in iniaerosol'
-    endif
-    call abort_physic("iniaerosl",'wrong number of aerosols',1)
-  endif ! of if (ia.ne.naerkind)
-
-  END SUBROUTINE iniaerosol
 
 end module aerosol_mod