Changeset 3158

Timestamp:

Dec 13, 2023, 4:09:00 PM (12 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

Files:

• changelog.txt (modified) (1 diff)
• libf/aeronomars/conduction.F (modified) (2 diffs)
• libf/aeronomars/euvheat.F90 (modified) (1 diff)
• libf/aeronomars/moldiff.F (modified) (2 diffs)
• libf/aeronomars/molvis.F (modified) (2 diffs)

trunk/LMDZ.MARS/changelog.txt

@@ -4404 +4404 @@
 * zzlay and zzlev are updated at the end of physics after call co2condens with updated pressure and temperature.
 * call concentrations, when photochem or callthermos is true, has been moved before the first calculation of zzlay and zzlev to be able to use varying reduced gas constant rnew.
+
+== 13/12/2023 == CS
+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.

trunk/LMDZ.MARS/libf/aeronomars/conduction.F

@@ -82 +82 @@
 
         zt(1)=pt(ig,1)+pdt(ig,1)*ptimestep
-c        zlay(1)=-log(pplay(ig,1)/pplev(ig,1))*Rnew(ig,1)*zt(1)/g
-c        zlev(1)=0.0
         zlay(1)=zzlay(ig,1)
         zlev(1)=zzlev(ig,1)
@@ -90 +88 @@
 
           zt(i)=pt(ig,i)+pdt(ig,i)*ptimestep 
-c          tmean=zt(i)
-c          if(zt(i).ne.zt(i-1))
-c     &    tmean=(zt(i)-zt(i-1))/log(zt(i)/zt(i-1))
-c          zlay(i)= zlay(i-1)
-c     &          -log(pplay(ig,i)/pplay(ig,i-1))*Rnew(ig,i-1)*tmean/g
-c          zlev(i)= zlev(i-1)
-c     &         -log(pplev(ig,i)/pplev(ig,i-1))*Rnew(ig,i-1)*tmean/g
-        zlay(i)=zzlay(ig,i)
-        zlev(i)=zzlev(ig,i)
+          zlay(i)=zzlay(ig,i)
+          zlev(i)=zzlev(ig,i)
         enddo
-        
-c        zlev(nlayer+1)= zlev(nlayer)
-c     &         -log(max(pplev(ig,nlayer+1),1.e-30)/pplev(ig,nlayer))
-c     &           *Rnew(ig,nlayer)*tmean/g
-c        if(pplev(ig,nlayer+1).eq.0.) 
-c     &     zlev(nlayer+1)=zlev(nlayer)+(zlay(nlayer)-zlay(nlayer-1))
-      
+            
         zlev(nlayer+1)= zlev(nlayer)+10000.
 

trunk/LMDZ.MARS/libf/aeronomars/euvheat.F90

@@ -411 +411 @@
          enddo
 
-!        zlocal(1)=-log(pplay(ig,1)/pplev(ig,1))
-!     &            *Rnew(ig,1)*zt(ig,1)/g
          zlocal(1)=zzlay(ig,1)
          zlocal(1)=zlocal(1)/1000.

trunk/LMDZ.MARS/libf/aeronomars/moldiff.F

@@ -308 +308 @@
 c
 
-c      zlocal(1)=-log(pplay(ig,1)/pplev(ig,1))* Rnew(ig,1)*tt(1)/g
       zlocal(1)=zzlay(ig,1)
 
@@ -325 +324 @@
         enddo
 
-c        tmean=tt(l)
-c        if(tt(l).ne.tt(l-1))
-c     &        tmean=(tt(l)-tt(l-1))/log(tt(l)/tt(l-1))
-c        zlocal(l)= zlocal(l-1)
-c     &         -log(pplay(ig,l)/pplay(ig,l-1))*rnew(ig,l)*tmean/g
       zlocal(l)=zzlay(ig,l)
       enddo 
 
-c      zlocal(nz)= zlocal(nz-1)
-c     &         -log(pplay(ig,nz)/pplay(ig,nz-1))*rnew(ig,nz)*tmean/g
       zlocal(nz)=zzlay(ig,nz)
         

trunk/LMDZ.MARS/libf/aeronomars/molvis.F

@@ -95 +95 @@
         zt(1)=pt(ig,1)+(pdteuv(ig,1)+pdtconduc(ig,1))*ptimestep
         zvel(1)=pvel(ig,1)
-c        zlay(1)=-log(pplay(ig,1)/pplev(ig,1))*Rnew(ig,1)*zt(1)/g
-c        zlev(1)=0.0
-
         zlay(1)=zzlay(ig,1)
         zlev(1)=zzlev(ig,1)
@@ -104 +101 @@
           zt(l)=pt(ig,l)+(pdteuv(ig,l)+pdtconduc(ig,l))*ptimestep
           zvel(l)=pvel(ig,l)
-c          tmean=zt(l)
-c          if(zt(l).ne.zt(l-1)) tmean=(zt(l)-zt(l-1))/log(zt(l)/zt(l-1))
-c          zlay(l)= zlay(l-1)
-c     &          -log(pplay(ig,l)/pplay(ig,l-1))*Rnew(ig,l-1)*tmean/g
-c          zlev(l)= zlev(l-1)
-c     &         -log(pplev(ig,l)/pplev(ig,l-1))*Rnew(ig,l-1)*tmean/g
-        zlay(l)=zzlay(ig,l)
-        zlev(l)=zzlev(ig,l)
+          zlay(l)=zzlay(ig,l)
+          zlev(l)=zzlev(ig,l)
         enddo
-
-c          zlev(nz+1)= zlev(nz)
-c     &         -log(max(pplev(ig,nz+1),1.e-30)/pplev(ig,nz))
-c     &          *Rnew(ig,nz)*tmean/g
-c          if(pplev(ig,nz+1).eq.0.)
-c     &       zlev(nz+1)=zlev(nz)+(zlay(nz)-zlay(nz-1))
 
           zlev(nz+1)= zlev(nz)+10000.