Changeset 3373 for trunk/LMDZ.PLUTO.old


Ignore:
Timestamp:
Jun 13, 2024, 2:52:10 PM (5 months ago)
Author:
afalco
Message:

Pluto.old PCM:
Corrected some runtime issues with gfortran.
AF

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

Legend:

Unmodified
Added
Removed

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

    r3357 r3373  
    5555      parameter(Nfine=701)
    5656      character(len=100) :: file_path
     57      character(len=100) :: file_name
    5758      real,save :: levdat(Nfine),densdat(Nfine)
    5859

  • trunk/LMDZ.PLUTO.old/libf/phypluto/callcorrk.F

    r3275 r3373  
    317317!     Prepare NON LTE correction in Pluto atmosphere
    318318      IF (nlte) then
    319     !     CALL nlte_ch4(ngrid,nlayer,nq,pplay,pplev,pt,vmrch4,
    320     ! &             eps_nlte_sw23,eps_nlte_sw33,eps_nlte_lw)
     319        CALL nlte_ch4(ngrid,nlayer,nq,pplay,pplev,pt,vmrch4,
     320    &             eps_nlte_sw23,eps_nlte_sw33,eps_nlte_lw)
    321321      ENDIF
    322322c     Net atmospheric radiative cooling rate from C2H2 (K.s-1):

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

    r3357 r3373  
    135135         hazerad_file="hazerad.txt"  ! default file
    136136         call getin("hazerad_file",hazerad_file)
    137          write(*,*) trim(rname)//" hazerad_file = ",trim(hazerad_file)
     137         write(*,*) " hazerad_file = ",trim(hazerad_file)
    138138         write(*,*) "Haze mmr datafile"
    139139         hazemmr_file="None"  ! default file
    140140         call getin("hazemmr_file",hazemmr_file)
    141          write(*,*) trim(rname)//" hazemmr_file = ",trim(hazemmr_file)
     141         write(*,*) " hazemmr_file = ",trim(hazemmr_file)
    142142         write(*,*) "Haze dens datafile"
    143143         hazedens_file="None"  ! default file
    144144         call getin("hazedens_file",hazedens_file)
    145          write(*,*) trim(rname)//" hazedens_file = ",trim(hazedens_file)
     145         write(*,*) " hazedens_file = ",trim(hazedens_file)
    146146
    147147!***************************************************************

  • trunk/LMDZ.PLUTO.old/libf/phypluto/optci.F90

    r3175 r3373  
    6666
    6767!     temporary variables for multiple aerosol calculation
    68       real*8 atemp, btemp
     68      real*8 atemp
     69      real*8 btemp(L_NLAYRAD,L_NSPECTI)
    6970
    7071!     variables for k in units m^-1
    … …  
    186187
    187188      DO NW=1,L_NSPECTI
    188          DO K=2,L_LEVELS+1
     189         DO K=2,L_LEVELS
    189190            do iaer=1,naerkind
    190191               TAUAEROLK(K,NW,IAER) = TAUAERO(K,IAER)*QSIAER(K,NW,IAER)
    … …  
    201202
    202203            atemp = 0.
    203             btemp = 0.
     204            btemp(L,NW) = 0.
     205            do iaer=1,naerkind
     206               atemp = atemp +                                     &
     207                     GIAER(K,NW,IAER)   * TAUAEROLK(K,NW,IAER) +    &
     208                     GIAER(K+1,NW,IAER) * TAUAEROLK(K+1,NW,IAER)
     209               btemp(L,NW) = btemp(L,NW) + TAUAEROLK(K,NW,IAER) + TAUAEROLK(K+1,NW,IAER)
     210!     *                    + 1.e-10
     211            end do
     212           
    204213            if(DTAUI(L,NW,NG) .GT. 1.0E-9) then
    205                do iaer=1,naerkind
    206                   atemp = atemp +                                     &
    207                       GIAER(K,NW,IAER)   * TAUAEROLK(K,NW,IAER) +    &
    208                       GIAER(K+1,NW,IAER) * TAUAEROLK(K+1,NW,IAER)
    209                   btemp = btemp + TAUAEROLK(K,NW,IAER) + TAUAEROLK(K+1,NW,IAER)
    210 !     *                    + 1.e-10
    211                end do
    212                WBARI(L,nw,ng) = btemp  / DTAUI(L,NW,NG)
     214               WBARI(L,nw,ng) = btemp(L,NW)  / DTAUI(L,NW,NG)
    213215            else
    214216               WBARI(L,nw,ng) = 0.0D0
    … …  
    216218            endif
    217219
    218             if(btemp .GT. 0.0) then
    219                cosbi(L,NW,NG) = atemp/btemp
     220            if(btemp(L,NW) .GT. 0.0) then
     221               cosbi(L,NW,NG) = atemp/btemp(L,NW)
    220222            else
    221223               cosbi(L,NW,NG) = 0.0D0
    … …  
    223225
    224226         END DO ! L vertical loop
     227         ! Last level
    225228         
    226      ! Last level
    227      
    228      L              = L_NLAYRAD
    229      K              = 2*L+1
    230      DTAUI(L,nw,ng) = DTAUKI(K,NW,NG) ! + 1.e-50
    231 
    232      atemp = 0.
    233      if(DTAUI(L,NW,NG) .GT. 1.0D-9) then
    234         do iaer=1,naerkind
    235            atemp = atemp + GIAER(K,NW,IAER)   * TAUAEROLK(K,NW,IAER)
    236         end do
    237         WBARI(L,nw,ng) = btemp  / DTAUI(L,NW,NG)
    238      else
    239         WBARI(L,nw,ng) = 0.0D0
    240         DTAUI(L,NW,NG) = 1.0D-9
    241      endif
    242 
    243      if(btemp .GT. 0.0d0) then
    244         cosbi(L,NW,NG) = atemp/btemp
    245      else
    246         cosbi(L,NW,NG) = 0.0D0
    247      end if
    248 
    249 !     Now the other Gauss points, if needed.
     229         L              = L_NLAYRAD
     230         K              = 2*L+1
     231         DTAUI(L,nw,ng) = DTAUKI(K,NW,NG) ! + 1.e-50
     232         btemp(L,NW) = 0
     233         do iaer=1,naerkind
     234               btemp(L,NW) = btemp(L,NW) + TAUAEROLK(K,NW,IAER)
     235         enddo
     236         
     237         atemp = 0.
     238         if(DTAUI(L,NW,NG) .GT. 1.0D-9) then
     239            do iaer=1,naerkind
     240               atemp = atemp + GIAER(K,NW,IAER)   * TAUAEROLK(K,NW,IAER)
     241            end do
     242            WBARI(L,nw,ng) = btemp(L,NW)  / DTAUI(L,NW,NG)
     243         else
     244            WBARI(L,nw,ng) = 0.0D0
     245            DTAUI(L,NW,NG) = 1.0D-9
     246         endif
     247
     248         if(btemp(L,NW) .GT. 0.0d0) then
     249            cosbi(L,NW,NG) = atemp/btemp(L,NW)
     250         else
     251            cosbi(L,NW,NG) = 0.0D0
     252         end if
     253
     254         !     Now the other Gauss points, if needed.
    250255
    251256         DO NG=1,L_NGAUSS-1
    … …  
    256261                  DTAUI(L,nw,ng) = DTAUKI(K,NW,NG)+DTAUKI(K+1,NW,NG)! + 1.e-50
    257262
    258                   btemp = 0.
    259263                  if(DTAUI(L,NW,NG) .GT. 1.0E-9) then
    260 
    261                      do iaer=1,naerkind
    262                         btemp = btemp + TAUAEROLK(K,NW,IAER) + TAUAEROLK(K+1,NW,IAER)
    263                      end do
    264                      WBARI(L,nw,ng) = btemp  / DTAUI(L,NW,NG)
    265 
     264                     WBARI(L,nw,ng) = btemp(L,NW)  / DTAUI(L,NW,NG)
    266265                  else
    267266                     WBARI(L,nw,ng) = 0.0D0

  • trunk/LMDZ.PLUTO.old/libf/phypluto/physiq.F

    r3329 r3373  
    11631163          zdqphot_prec(:,:)=0.
    11641164          zdqphot_ch4(:,:)=0.
     1165          zdqhaze(:,:,:)=0.
    11651166          ! Forcing to a fixed haze profile if haze_proffix
    11661167          if (haze_proffix.and.i_haze.gt.0.) then
Note: See TracChangeset for help on using the changeset viewer.