Changeset 890 for trunk/LMDZ.MARS/libf

Timestamp:

Feb 27, 2013, 3:51:29 PM (12 years ago)

Author:

emillour

Message:

Mars GCM:
minor bug correction in newcondens.F and some code tidying (in co2snow.F also).
EM

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• co2snow.F (modified) (2 diffs)
• newcondens.F (modified) (7 diffs)

trunk/LMDZ.MARS/libf/phymars/co2snow.F

@@ -20 +20 @@
 #include "callkeys.h"
 
-c     input
-c     -----
+c     Arguments
+c     ---------
 
-      INTEGER ngrid,nlayer
-      REAL ptimestep , emisref(ngrid) ! grd or ice  emissivity without snow
-      logical condsub(ngrid)
+      INTEGER,INTENT(IN) :: ngrid ! number of atmospheric columns
+      INTEGER,INTENT(IN) :: nlayer ! number of atmospheric layers
+      REAL,INTENT(IN) :: ptimestep ! timestep of the physics (s)
+      REAL,INTENT(IN) :: emisref(ngrid) ! grd or ice  emissivity without snow
+      logical,intent(in) :: condsub(ngrid) ! true if there is CO2 condensation
+                                           ! or sublimation in the column
+      REAL,INTENT(IN) :: pplev(ngrid,nlayer+1) ! interlayer pressure (Pa)
+      REAL,INTENT(IN) :: pcondicea(ngrid,nlayer) ! CO2 condensation rate (kg/m2/s)
+      REAL,INTENT(IN) :: pcondices(ngrid) ! CO2 condensation rate (kg/m2/s)
+                                          ! on the surface
+      REAL,INTENT(IN) :: pfallice(ngrid,nlayer+1) ! falling CO2 ice (kg/m2/s)
 
-      REAL pplev(ngrid,nlayer+1)
-      REAL pcondicea(ngrid,nlayer), pcondices(ngrid)
-      REAL pfallice(ngrid,nlayer+1)
-
-c     output
-c     ------
-
-      REAL pemisurf(ngrid)
+      REAL,INTENT(OUT) :: pemisurf(ngrid) ! surface emissivity
 
 c     local
@@ -45 +46 @@
 c     saved
 c     -----
-      REAL Kscat(2), scaveng
-      LOGICAL firstcall
-      SAVE firstcall
-      save Kscat , scaveng
-      DATA firstcall/.true./
+      REAL,save :: Kscat(2), scaveng
+      LOGICAL,SAVE :: firstcall=.true.
 
 c --------------

trunk/LMDZ.MARS/libf/phymars/newcondens.F

@@ -66 +66 @@
 c    Arguments :
 c    ---------
-      INTEGER ngrid, nlayer, nq
-
-      REAL ptimestep 
-      REAL pplay(ngrid,nlayer),pplev(ngrid,nlayer+1)
-      REAL pcapcal(ngrid)
-      REAL pt(ngrid,nlayer)
-      REAL ptsrf(ngrid)
-      REAL pphi(ngrid,nlayer)
-      REAL pdt(ngrid,nlayer),pdtsrf(ngrid),pdtc(ngrid,nlayer)
-      REAL pdtsrfc(ngrid),pdpsrf(ngrid)
-      REAL piceco2(ngrid),psolaralb(ngrid,2),pemisurf(ngrid)
-
-      REAL pu(ngrid,nlayer) ,  pv(ngrid,nlayer)
-      REAL pdu(ngrid,nlayer) , pdv(ngrid,nlayer)
-      REAL pduc(ngrid,nlayer) , pdvc(ngrid,nlayer)
-      REAL pq(ngridmx,nlayer,nq),pdq(ngrid,nlayer,nq)
-      REAL pdqc(ngrid,nlayer,nq)
-      REAL fluxsurf_sw(ngrid,2) ! added to calculate flux dependent albedo 
-      real zls ! solar longitude (rad)
+      INTEGER,INTENT(IN) :: ngrid  ! number of atmospheric columns
+      INTEGER,INTENT(IN) :: nlayer ! number of vertical layers
+      INTEGER,INTENT(IN) :: nq  ! number of tracers
+
+      REAL,INTENT(IN) :: ptimestep ! physics timestep (s)
+      REAL,INTENT(IN) :: pcapcal(ngrid)
+      REAL,INTENT(IN) :: pplay(ngrid,nlayer) !mid-layer pressure (Pa)
+      REAL,INTENT(IN) :: pplev(ngrid,nlayer+1) ! inter-layer pressure (Pa)
+      REAL,INTENT(IN) :: ptsrf(ngrid) ! surface temperature (K)
+      REAL,INTENT(IN) :: pt(ngrid,nlayer) ! atmospheric temperature (K)
+      REAL,INTENT(IN) :: pphi(ngrid,nlayer) ! geopotential (m2.s-2)
+      REAL,INTENT(IN) :: pdt(ngrid,nlayer) ! tendency on temperature from
+                                           ! previous physical processes (K/s)
+      REAL,INTENT(IN) :: pdu(ngrid,nlayer) ! tendency on zonal wind (m/s2)
+                                           ! from previous physical processes
+      REAL,INTENT(IN) :: pdv(ngrid,nlayer) ! tendency on meridional wind (m/s2)
+                                           ! from previous physical processes
+      REAL,INTENT(IN) :: pdtsrf(ngrid) ! tendency on surface temperature from
+                                       ! previous physical processes (K/s)
+      REAL,INTENT(IN) :: pu(ngrid,nlayer) ! zonal wind (m/s)
+      REAL,INTENT(IN) :: pv(ngrid,nlayer) ! meridional wind (m/s)
+      REAL,INTENT(IN) :: pq(ngrid,nlayer,nq) ! tracers (../kg_air)
+      REAL,INTENT(IN) :: pdq(ngrid,nlayer,nq) ! tendency on tracers from
+                                              ! previous physical processes
+      REAL,INTENT(INOUT) :: piceco2(ngrid) ! CO2 ice on the surface (kg.m-2)
+      REAL,INTENT(INOUT) :: psolaralb(ngrid,2) ! albedo of the surface
+      REAL,INTENT(INOUT) :: pemisurf(ngrid) ! emissivity of the surface
+      
+      ! tendencies due to CO2 condensation/sublimation:
+      REAL,INTENT(OUT) :: pdtc(ngrid,nlayer) ! tendency on temperature (K/s)
+      REAL,INTENT(OUT) :: pdtsrfc(ngrid) ! tendency on surface temperature (K/s)
+      REAL,INTENT(OUT) :: pdpsrf(ngrid) ! tendency on surface pressure (Pa/s)
+      REAL,INTENT(OUT) :: pduc(ngrid,nlayer) ! tendency on zonal wind (m.s-2)
+      REAL,INTENT(OUT) :: pdvc(ngrid,nlayer) ! tendency on meridional wind (m.s-2)
+      REAL,INTENT(OUT) :: pdqc(ngrid,nlayer,nq) ! tendency on tracers
+      
+      ! added to calculate flux dependent albedo:
+      REAL,intent(in) :: fluxsurf_sw(ngrid,2)
+      real,intent(in) :: zls ! solar longitude (rad)
+
 c
 c    Local variables :
@@ -100 +121 @@
 c   the incident solar flux      
       PARAMETER(fluxdependent=.false.)
-      REAL slopy,alpha,constA,constB,constT,albediceF_new(ngridmx)
-      REAL zt(ngridmx,nlayermx)
+      REAL slopy,alpha,constA,constB,constT,albediceF_new(ngrid)
+      REAL zt(ngrid,nlayer)
       REAL zcpi
-      REAL ztcond (ngridmx,nlayermx+1)
-      REAL ztcondsol(ngridmx) 
-      REAL zdiceco2(ngridmx)
-      REAL zcondicea(ngridmx,nlayermx)
-      REAL zcondices(ngridmx)
-      REAL zfallice(ngridmx,nlayermx+1) , zfallheat
-      REAL zmflux(nlayermx+1)
-      REAL zu(nlayermx),zv(nlayermx)
-      REAL zq(nlayermx,nqmx),zq1(nlayermx)
-      REAL ztsrf(ngridmx)
-      REAL ztc(nlayermx), ztm(nlayermx+1) 
-      REAL zum(nlayermx+1) , zvm(nlayermx+1)
-      REAL zqm(nlayermx+1,nqmx),zqm1(nlayermx+1)
-      REAL masse(nlayermx),w(nlayermx+1)
-      REAL Sm(nlayermx),Smq(nlayermx,nqmx),mixmas,qmix
-      LOGICAL condsub(ngridmx)
+      REAL ztcond (ngrid,nlayer+1)
+      REAL ztcondsol(ngrid) 
+      REAL zdiceco2(ngrid)
+      REAL zcondicea(ngrid,nlayer)
+      REAL zcondices(ngrid)
+      REAL zfallice(ngrid,nlayer+1) , zfallheat
+      REAL zmflux(nlayer+1)
+      REAL zu(nlayer),zv(nlayer)
+      REAL zq(nlayer,nq),zq1(nlayer)
+      REAL ztsrf(ngrid)
+      REAL ztc(nlayer), ztm(nlayer+1) 
+      REAL zum(nlayer+1) , zvm(nlayer+1)
+      REAL zqm(nlayer+1,nq),zqm1(nlayer+1)
+      REAL masse(nlayer),w(nlayer+1)
+      REAL Sm(nlayer),Smq(nlayer,nq),mixmas,qmix
+      LOGICAL condsub(ngrid)
 
 c variable speciale diagnostique
-      real tconda1(ngridmx,nlayermx)
-      real tconda2(ngridmx,nlayermx)
-c     REAL zdiceco2a(ngridmx) ! for diagnostic only
-      real zdtsig (ngridmx,nlayermx)
-      real zdt (ngridmx,nlayermx)
-      real vmr_co2(ngridmx,nlayermx) ! co2 volume mixing ratio
+      real tconda1(ngrid,nlayer)
+      real tconda2(ngrid,nlayer)
+c     REAL zdiceco2a(ngrid) ! for diagnostic only
+      real zdtsig (ngrid,nlayer)
+      real zdt (ngrid,nlayer)
+      real vmr_co2(ngrid,nlayer) ! co2 volume mixing ratio
 ! improved_ztcond flag: If set to .true. (AND running with a 'co2' tracer)
 ! then condensation temperature is computed using partial pressure of CO2
@@ -137 +158 @@
 
 c   local saved variables
-      integer ico2
-      real qco2min,qco2,mmean
-      save ico2,qco2min
-      REAL emisref(ngridmx)
-      REAL latcond,tcond1mb
-      REAL acond,bcond,ccond,cpice
-!      REAL albediceF(ngridmx)
-      SAVE emisref
-      SAVE latcond,acond,bcond,ccond,cpice
-!      SAVE albediceF
-      real m_co2, m_noco2, A , B
-      save A, B, m_co2, m_noco2
-
-      LOGICAL firstcall !,firstcall2
-      SAVE firstcall !,firstcall2
-!      REAL SSUM
-!      EXTERNAL SSUM
-
-      DATA latcond,tcond1mb/5.9e5,136.27/
-      DATA cpice /1000./
-      DATA firstcall/.true./
-!      DATA firstcall2/.true./
+      integer,save :: ico2 ! index of CO2 tracer
+      real,save :: qco2min,qco2,mmean
+      real,allocatable,save :: emisref(:)
+      real,parameter :: latcond=5.9e5 ! (J/kg) Latent heat of solid CO2 ice
+      real,parameter :: tcond1mb=136.27 ! condensation temperature (K) at 1 mbar
+      real,parameter :: cpice=1000. ! (J.kg-1.K-1) specific heat of CO2 ice
+      REAL,SAVE :: acond,bcond,ccond
+!      REAL,SAVE :: albediceF(ngridmx)
+      real,save :: m_co2, m_noco2, A , B
+
+      LOGICAL,SAVE :: firstcall = .true. !,firstcall2=.true.
 
       integer flag
@@ -168 +178 @@
 c
       IF (firstcall) THEN
+         
+         allocate(emisref(ngrid))
+         
          bcond=1./tcond1mb
          ccond=cpp/(g*latcond)
@@ -198 +211 @@
            qco2min =0.75  
          end if
-      ENDIF
+      ENDIF ! of IF (firstcall)
       zcpi=1./cpp
 c
@@ -348 +361 @@
                  zfallheat=zfallice(ig,1)*
      &           (pphi(ig,1)- phisfi(ig) +
-     &           cpice*(ztcond(ig,1)-ztcondsol(ig)))/(latcond*ptimestep)
+     &           cpice*(ztcond(ig,1)-ztcondsol(ig)))/latcond
             ELSE
                  zfallheat=0.
@@ -733 +746 @@
 !! Specific stuff to bound co2 tracer values ....
       if (bound_qco2.and.(ico2.ne.0)) then
-        do ig=1,ngridmx
-          do l=1,nlayermx
+        do ig=1,ngrid
+          do l=1,nlayer
             zqco2=pq(ig,l,ico2)
      &            +(pdq(ig,l,ico2)+pdqc(ig,l,ico2))*ptimestep