Changeset 3232 for trunk/LMDZ.PLUTO

Timestamp:

Feb 21, 2024, 5:45:11 PM (9 months ago)

Author:

afalco

Message:

Pluto PCM:
1D functional
AF

Location:

trunk/LMDZ.PLUTO/libf/phypluto

Files:

• condense_n2.F90 (modified) (2 diffs)
• dyn1d/calcenergy_kcm.F90 (deleted)
• dyn1d/kcmprof_fn.F90 (modified) (3 diffs)
• dyn1d/rcm1d.F (modified) (9 diffs)
• physiq_mod.F90 (modified) (1 diff)

trunk/LMDZ.PLUTO/libf/phypluto/condense_n2.F90

@@ -495 +495 @@
           end do
 
+!             Value transfert at the surface interface when condensation sublimation:
+
+          if (zmflux(1).lt.0) then
+!               Surface condensation
+            zum(1)= zu(1)
+            zvm(1)= zv(1)
+            ztm(1) = ztc(1)
+          else
+!               Surface sublimation:
+            ztm(1) = ztsrf(ig) + pdtsrfc(ig)*ptimestep
+            zum(1) = 0
+            zvm(1) = 0
+          end if
+          do iq=1,nq
+              zqm(1,iq)=0. ! most tracer do not condense !
+          enddo
+!             Special case if the tracer is n2 gas
+          if (igcm_n2.ne.0) zqm(1,igcm_n2)=1.
+
            ztm(2:klev+1)=0.
            zum(2:klev+1)=0.
@@ -528 +547 @@
            end if
 
-!             Value transfert at the surface interface when condensation sublimation:
-
-          if (zmflux(1).lt.0) then
-!               Surface condensation
-            zum(1)= zu(1)
-            zvm(1)= zv(1)
-            ztm(1) = ztc(1)
-          else
-!               Surface sublimation:
-            ztm(1) = ztsrf(ig) + pdtsrfc(ig)*ptimestep
-            zum(1) = 0
-            zvm(1) = 0
-          end if
-          do iq=1,nq
-             zqm(1,iq)=0. ! most tracer do not condense !
-          enddo
-!             Special case if the tracer is n2 gas
-          if (igcm_n2.ne.0) zqm(1,igcm_n2)=1.
 
           !!! Source haze: 0.02 pourcent when n2 sublimes

trunk/LMDZ.PLUTO/libf/phypluto/dyn1d/kcmprof_fn.F90

@@ -2 +2 @@
 
 use params_h
-use watercommon_h, only : mH2O
 use gases_h
 use comcstfi_mod, only: mugaz, cpp, g
-use callkeys_mod, only: co2cond
+use callkeys_mod, only: n2cond
 implicit none
 
@@ -97 +96 @@
      m_v   = 1.0
      tcrit = 1000.0
-  elseif(trim(gnom(vgas)).eq.'H2O')then
-     m_v   = dble(mH2O/1000.)
-     tcrit = 6.47d2
+!   elseif(trim(gnom(vgas)).eq.'H2O')then
+   !   m_v   = dble(mH2O/1000.) ! AF24: no water in Pluton
+   !   tcrit = 6.47d2
   elseif(trim(gnom(vgas)).eq.'NH3')then
      m_v   = 17.031/1000.
@@ -368 +367 @@
   enddo
 
-!    CO2 condensation 'haircut' of temperature profile if necessary
-  if(co2cond)then
-     print*,'CO2 condensation haircut - assumes CO2-dominated atmosphere!'
-     do ilay=2,nlayer
-        if(P_rcm(ilay).lt.518000.)then
-           TCO2cond = (-3167.8)/(log(.01*P_rcm(ilay))-23.23) ! Fanale's formula
-        else
-           TCO2cond = 684.2-92.3*log(P_rcm(ilay))+4.32*log(P_rcm(ilay))**2 
-           ! liquid-vapour transition (based on CRC handbook 2003 data)
-        endif
-
-        print*,'p=',P_rcm(ilay),', T=',T_rcm(ilay),' Tcond=',TCO2cond
-        if(T_rcm(ilay).lt.TCO2cond)then
-           T_rcm(ilay)=TCO2cond
-        endif
-     enddo
-  endif
-
   return
 end subroutine kcmprof_fn

trunk/LMDZ.PLUTO/libf/phypluto/dyn1d/rcm1d.F

@@ -89 +89 @@
       REAL,ALLOCATABLE :: qsurf(:)    ! tracer surface budget (e.g. kg.m-2)
       REAL,ALLOCATABLE :: tsoil(:)    ! subsurface soil temperature (K)
-!      REAL co2ice               ! co2ice layer (kg.m-2) !not used anymore
-      integer :: i_co2_ice=0     ! tracer index of co2 ice
+!      REAL n2ice               ! n2ice layer (kg.m-2) !not used anymore
+      integer :: i_n2_ice=0     ! tracer index of n2 ice
       integer :: i_h2o_ice=0     ! tracer index of h2o ice
       integer :: i_h2o_vap=0     ! tracer index of h2o vapor
@@ -235 +235 @@
       emisice(2)=0.95            ! Emissivite calotte sud  
 
-      iceradius(1) = 100.e-6     ! mean scat radius of CO2 snow (north)
-      iceradius(2) = 100.e-6     ! mean scat radius of CO2 snow (south)
+      iceradius(1) = 100.e-6     ! mean scat radius of n2 snow (north)
+      iceradius(2) = 100.e-6     ! mean scat radius of n2 snow (south)
       dtemisice(1) = 2.          ! time scale for snow metamorphism (north)
       dtemisice(2) = 2.          ! time scale for snow metamorphism (south
@@ -313 +313 @@
             endif
           enddo !of do iq=1,nq
-! check for co2_ice / h2o_ice tracers:
-         i_co2_ice=0
+! check for n2_ice / h2o_ice tracers:
+         i_n2_ice=0
          i_h2o_ice=0
          i_h2o_vap=0
          do iq=1,nq
-           if (tname(iq)=="co2_ice") then
-             i_co2_ice=iq
+           if (tname(iq)=="n2") then
+             i_n2_ice=iq
            elseif (tname(iq)=="h2o_ice") then
              i_h2o_ice=iq
@@ -390 +390 @@
       call getin("rad",rad)
       ! Planetary  radius is needed to compute shadow of the rings
-      IF (rad.eq.-99999. .and. rings_shadow .eqv. .true.) THEN
+      IF (rad.eq.-99999.) THEN
           PRINT *,"STOP. I NEED rad IN RCM1D.DEF."
           STOP
@@ -623 +623 @@
             write(*,*)"  tracer:",trim(txt)
               
-            ! CO2
-            if (txt.eq."co2_ice") then
+            ! n2
+            if (txt.eq."n2_ice") then
                q(:,iq)=0.   ! kg/kg of atmosphere
                qsurf(iq)=0. ! kg/m2 at the surface               
-               ! Look for a "profile_co2_ice" input file
-               open(91,file='profile_co2_ice',status='old',
+               ! Look for a "profile_n2_ice" input file
+               open(91,file='profile_n2_ice',status='old',
      &         form='formatted',iostat=ierr)
                if (ierr.eq.0) then
@@ -636 +636 @@
                   enddo
                else
-                  write(*,*) "No profile_co2_ice file!"
+                  write(*,*) "No profile_n2_ice file!"
                endif
                close(91)
-            endif ! of if (txt.eq."co2")
+            endif ! of if (txt.eq."n2")
           
             ! WATER VAPOUR
@@ -739 +739 @@
       ENDDO
 
-c  emissivity / surface co2 ice ( + h2o ice??)
+c  emissivity / surface n2 ice ( + h2o ice??)
 c  -------------------------------------------
       emis(1)=emissiv ! default value for emissivity
@@ -745 +745 @@
       call getin("emis",emis(1))
       write(*,*) " emis = ",emis(1)
-      emissiv=emis(1) ! we do this so that condense_co2 sets things to the right
+      emissiv=emis(1) ! we do this so that condense_n2 sets things to the right
                    ! value if there is no snow
 
-      if(i_co2_ice.gt.0)then
-         qsurf(i_co2_ice)=0 ! default value for co2ice
-         print*,'Initial CO2 ice on the surface (kg.m-2)'
-         call getin("co2ice",qsurf(i_co2_ice))
-         write(*,*) " co2ice = ",qsurf(i_co2_ice)
-         IF (qsurf(i_co2_ice).ge.1.E+0) THEN
-            ! if we have some CO2 ice on the surface, change emissivity
+      if(i_n2_ice.gt.0)then
+         qsurf(i_n2_ice)=0 ! default value for n2ice
+         print*,'Initial n2 ice on the surface (kg.m-2)'
+         call getin("n2ice",qsurf(i_n2_ice))
+         write(*,*) " n2ice = ",qsurf(i_n2_ice)
+         IF (qsurf(i_n2_ice).ge.1.E+0) THEN
+            ! if we have some n2 ice on the surface, change emissivity
             if (latitude(1).ge.0) then ! northern hemisphere
               emis(1)=emisice(1)
@@ -911 +911 @@
       call physdem1("startfi.nc",nsoilmx,1,llm,nq,
      &                dtphys,time,
-     &                tsurf,tsoil,emis,q2,qsurf,
-     &                cloudfrac,totcloudfrac,hice, 
-     &                rnat,pctsrf_sic,tslab,tsea_ice,sea_ice)
+     &                tsurf,tsoil,emis,q2,qsurf)
 
 c=======================================================================

trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90

@@ -551 +551 @@
              tsoil(1:ngrid,isoil)=tsurf(1:ngrid)
            enddo
-           if (is_master) write(*,*) "Physiq: initializing day_ini to pdat !"
+           if (is_master) write(*,*) "Physiq: initializing day_ini to pday !"
            day_ini=pday
          endif