Ignore:
Timestamp:
Dec 13, 2023, 4:09:00 PM (13 months ago)
Author:
csegonne
Message:

MARS PCM
Cleaning of conduction.F, euvheat.F90, moldiff.F and molvis.F, some commented lines referring to a local calculation of layers/levels altitudes have been removed.

Location:
trunk/LMDZ.MARS/libf/aeronomars
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • TabularUnified trunk/LMDZ.MARS/libf/aeronomars/conduction.F

    r2615 r3158  
    8282
    8383        zt(1)=pt(ig,1)+pdt(ig,1)*ptimestep
    84 c        zlay(1)=-log(pplay(ig,1)/pplev(ig,1))*Rnew(ig,1)*zt(1)/g
    85 c        zlev(1)=0.0
    8684        zlay(1)=zzlay(ig,1)
    8785        zlev(1)=zzlev(ig,1)
     
    9088
    9189          zt(i)=pt(ig,i)+pdt(ig,i)*ptimestep
    92 c          tmean=zt(i)
    93 c          if(zt(i).ne.zt(i-1))
    94 c     &    tmean=(zt(i)-zt(i-1))/log(zt(i)/zt(i-1))
    95 c          zlay(i)= zlay(i-1)
    96 c     &          -log(pplay(ig,i)/pplay(ig,i-1))*Rnew(ig,i-1)*tmean/g
    97 c          zlev(i)= zlev(i-1)
    98 c     &         -log(pplev(ig,i)/pplev(ig,i-1))*Rnew(ig,i-1)*tmean/g
    99         zlay(i)=zzlay(ig,i)
    100         zlev(i)=zzlev(ig,i)
     90          zlay(i)=zzlay(ig,i)
     91          zlev(i)=zzlev(ig,i)
    10192        enddo
    102        
    103 c        zlev(nlayer+1)= zlev(nlayer)
    104 c     &         -log(max(pplev(ig,nlayer+1),1.e-30)/pplev(ig,nlayer))
    105 c     &           *Rnew(ig,nlayer)*tmean/g
    106 c        if(pplev(ig,nlayer+1).eq.0.)
    107 c     &     zlev(nlayer+1)=zlev(nlayer)+(zlay(nlayer)-zlay(nlayer-1))
    108      
     93           
    10994        zlev(nlayer+1)= zlev(nlayer)+10000.
    11095
  • TabularUnified trunk/LMDZ.MARS/libf/aeronomars/euvheat.F90

    r2615 r3158  
    411411         enddo
    412412
    413 !        zlocal(1)=-log(pplay(ig,1)/pplev(ig,1))
    414 !     &            *Rnew(ig,1)*zt(ig,1)/g
    415413         zlocal(1)=zzlay(ig,1)
    416414         zlocal(1)=zlocal(1)/1000.
  • TabularUnified trunk/LMDZ.MARS/libf/aeronomars/moldiff.F

    r2615 r3158  
    308308c
    309309
    310 c      zlocal(1)=-log(pplay(ig,1)/pplev(ig,1))* Rnew(ig,1)*tt(1)/g
    311310      zlocal(1)=zzlay(ig,1)
    312311
     
    325324        enddo
    326325
    327 c        tmean=tt(l)
    328 c        if(tt(l).ne.tt(l-1))
    329 c     &        tmean=(tt(l)-tt(l-1))/log(tt(l)/tt(l-1))
    330 c        zlocal(l)= zlocal(l-1)
    331 c     &         -log(pplay(ig,l)/pplay(ig,l-1))*rnew(ig,l)*tmean/g
    332326      zlocal(l)=zzlay(ig,l)
    333327      enddo
    334328
    335 c      zlocal(nz)= zlocal(nz-1)
    336 c     &         -log(pplay(ig,nz)/pplay(ig,nz-1))*rnew(ig,nz)*tmean/g
    337329      zlocal(nz)=zzlay(ig,nz)
    338330       
  • TabularUnified trunk/LMDZ.MARS/libf/aeronomars/molvis.F

    r2615 r3158  
    9595        zt(1)=pt(ig,1)+(pdteuv(ig,1)+pdtconduc(ig,1))*ptimestep
    9696        zvel(1)=pvel(ig,1)
    97 c        zlay(1)=-log(pplay(ig,1)/pplev(ig,1))*Rnew(ig,1)*zt(1)/g
    98 c        zlev(1)=0.0
    99 
    10097        zlay(1)=zzlay(ig,1)
    10198        zlev(1)=zzlev(ig,1)
     
    104101          zt(l)=pt(ig,l)+(pdteuv(ig,l)+pdtconduc(ig,l))*ptimestep
    105102          zvel(l)=pvel(ig,l)
    106 c          tmean=zt(l)
    107 c          if(zt(l).ne.zt(l-1)) tmean=(zt(l)-zt(l-1))/log(zt(l)/zt(l-1))
    108 c          zlay(l)= zlay(l-1)
    109 c     &          -log(pplay(ig,l)/pplay(ig,l-1))*Rnew(ig,l-1)*tmean/g
    110 c          zlev(l)= zlev(l-1)
    111 c     &         -log(pplev(ig,l)/pplev(ig,l-1))*Rnew(ig,l-1)*tmean/g
    112         zlay(l)=zzlay(ig,l)
    113         zlev(l)=zzlev(ig,l)
     103          zlay(l)=zzlay(ig,l)
     104          zlev(l)=zzlev(ig,l)
    114105        enddo
    115 
    116 c          zlev(nz+1)= zlev(nz)
    117 c     &         -log(max(pplev(ig,nz+1),1.e-30)/pplev(ig,nz))
    118 c     &          *Rnew(ig,nz)*tmean/g
    119 c          if(pplev(ig,nz+1).eq.0.)
    120 c     &       zlev(nz+1)=zlev(nz)+(zlay(nz)-zlay(nz-1))
    121106
    122107          zlev(nz+1)= zlev(nz)+10000.
Note: See TracChangeset for help on using the changeset viewer.