Changeset 3357 for trunk


Ignore:
Timestamp:
Jun 4, 2024, 5:56:15 PM (7 months ago)
Author:
afalco
Message:

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

Location:
trunk/LMDZ.PLUTO.old/libf/phypluto
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.PLUTO.old/libf/phypluto/aerosol_mod.F90

    r3353 r3357  
    6262      IF (firstcall) then
    6363        firstcall=.false.
    64       !   file_path=trim(datadir)//'/haze_prop/hazemmr.txt'
    65         file_path=trim(datadir)//'/haze_prop/dens_haze.txt'
     64        if (hazemmr_file/='None')then
     65            file_name = hazemmr_file
     66            print*, 'Read Haze MMR file: ',hazemmr_file
     67        else if(hazedens_file/='None')then
     68            file_name = hazedens_file
     69            print*, 'Read Haze density: ',hazedens_file
     70        else
     71            STOP "No filename given for haze profile. Either set hazemmr_file or hazedens_file"
     72        endif
     73
     74        file_path=trim(datadir)//'/haze_prop/'//file_name
    6675        open(224,file=file_path,form='formatted')
    6776        do ifine=1,Nfine
     
    8089      !                                --> kg m-3 --> kg/kg
    8190      do iaer=1,naerkind
    82             if(iaer.eq.iaero_haze.and.1.eq.1) then !TB22 activate/deactivate mmr or part density   
     91            if(iaer.eq.iaero_haze.and.hazedens_file/='None') then !AF24 activate/deactivate mmr or part density
    8392              !print*, 'Haze profile is fixed'
    8493              do ig=1,ngrid
  • trunk/LMDZ.PLUTO.old/libf/phypluto/datafile_mod.F90

    r3353 r3357  
    1818      ! Data haze properties
    1919      character(len=300),save :: hazeprop
     20      character(len=300),save :: hazerad_file
     21      character(len=300),save :: hazemmr_file
     22      character(len=300),save :: hazedens_file
    2023
    2124      end module datafile_mod
  • trunk/LMDZ.PLUTO.old/libf/phypluto/inifis.F

    r3175 r3357  
    132132         call getin("hazeprop",hazeprop)
    133133         write(*,*) " hazeprop = ",trim(hazeprop)
     134         write(*,*) "Haze radii datafile"
     135         hazerad_file="hazerad.txt"  ! default file
     136         call getin("hazerad_file",hazerad_file)
     137         write(*,*) trim(rname)//" hazerad_file = ",trim(hazerad_file)
     138         write(*,*) "Haze mmr datafile"
     139         hazemmr_file="None"  ! default file
     140         call getin("hazemmr_file",hazemmr_file)
     141         write(*,*) trim(rname)//" hazemmr_file = ",trim(hazemmr_file)
     142         write(*,*) "Haze dens datafile"
     143         hazedens_file="None"  ! default file
     144         call getin("hazedens_file",hazedens_file)
     145         write(*,*) trim(rname)//" hazedens_file = ",trim(hazedens_file)
    134146
    135147!***************************************************************
  • trunk/LMDZ.PLUTO.old/libf/phypluto/radii_mod.F90

    r3353 r3357  
    9191      IF (firstcall) then
    9292        firstcall=.false.
    93         file_path=trim(datadir)//'/haze_prop/rad_haze.txt'
     93         file_path=trim(datadir)//'/haze_prop/'//hazerad_file
    9494      !   file_path=trim(datadir)//'/haze_prop/hazerad.txt'
    9595        open(223,file=file_path,form='formatted')
Note: See TracChangeset for help on using the changeset viewer.