Changeset 3357 for trunk/LMDZ.PLUTO.old

Timestamp:

Jun 4, 2024, 5:56:15 PM (6 months ago)

Author:

afalco

Message:

Pluto.old PCM:
Read haze mmr, rad and dens files.
AF

Location:

trunk/LMDZ.PLUTO.old/libf/phypluto

Files:

• aerosol_mod.F90 (modified) (2 diffs)
• datafile_mod.F90 (modified) (1 diff)
• inifis.F (modified) (1 diff)
• radii_mod.F90 (modified) (1 diff)

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

@@ -62 +62 @@
       IF (firstcall) then
         firstcall=.false.
-      !   file_path=trim(datadir)//'/haze_prop/hazemmr.txt'
-        file_path=trim(datadir)//'/haze_prop/dens_haze.txt'
+        if (hazemmr_file/='None')then
+            file_name = hazemmr_file
+            print*, 'Read Haze MMR file: ',hazemmr_file
+        else if(hazedens_file/='None')then
+            file_name = hazedens_file
+            print*, 'Read Haze density: ',hazedens_file
+        else
+            STOP "No filename given for haze profile. Either set hazemmr_file or hazedens_file"
+        endif
+
+        file_path=trim(datadir)//'/haze_prop/'//file_name
         open(224,file=file_path,form='formatted')
         do ifine=1,Nfine
@@ -80 +89 @@
       !                                --> kg m-3 --> kg/kg
       do iaer=1,naerkind
-            if(iaer.eq.iaero_haze.and.1.eq.1) then !TB22 activate/deactivate mmr or part density   
+            if(iaer.eq.iaero_haze.and.hazedens_file/='None') then !AF24 activate/deactivate mmr or part density
               !print*, 'Haze profile is fixed'
               do ig=1,ngrid

trunk/LMDZ.PLUTO.old/libf/phypluto/datafile_mod.F90

@@ -18 +18 @@
       ! Data haze properties
       character(len=300),save :: hazeprop 
+      character(len=300),save :: hazerad_file 
+      character(len=300),save :: hazemmr_file 
+      character(len=300),save :: hazedens_file 
 
       end module datafile_mod

trunk/LMDZ.PLUTO.old/libf/phypluto/inifis.F

@@ -132 +132 @@
          call getin("hazeprop",hazeprop) 
          write(*,*) " hazeprop = ",trim(hazeprop)
+         write(*,*) "Haze radii datafile"
+         hazerad_file="hazerad.txt"  ! default file
+         call getin("hazerad_file",hazerad_file)
+         write(*,*) trim(rname)//" hazerad_file = ",trim(hazerad_file)
+         write(*,*) "Haze mmr datafile"
+         hazemmr_file="None"  ! default file
+         call getin("hazemmr_file",hazemmr_file)
+         write(*,*) trim(rname)//" hazemmr_file = ",trim(hazemmr_file)
+         write(*,*) "Haze dens datafile"
+         hazedens_file="None"  ! default file
+         call getin("hazedens_file",hazedens_file)
+         write(*,*) trim(rname)//" hazedens_file = ",trim(hazedens_file)
 
 !***************************************************************

trunk/LMDZ.PLUTO.old/libf/phypluto/radii_mod.F90

@@ -91 +91 @@
       IF (firstcall) then
         firstcall=.false.
-        file_path=trim(datadir)//'/haze_prop/rad_haze.txt'
+         file_path=trim(datadir)//'/haze_prop/'//hazerad_file
       !   file_path=trim(datadir)//'/haze_prop/hazerad.txt'
         open(223,file=file_path,form='formatted')