Changeset 353

Timestamp:

Nov 7, 2011, 1:05:27 PM (13 years ago)

Author:

acolaitis

Message:

Emissivity of surface when co2 ice is present has been changed to 1 (for lwxb.F and lwxd.F) for gaz exclusively. (It remains the decreased value for lwdiff.F).

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• callradite.F (modified) (3 diffs)
• lwmain.F (modified) (5 diffs)
• physiq.F (modified) (2 diffs)

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

@@ -2 +2 @@
      $     emis,mu0,pplev,pplay,pt,tsurf,fract,dist_sol,igout,
      $     dtlw,dtsw,fluxsurf_lw,fluxsurf_sw,fluxtop_lw,fluxtop_sw,
-     &     tauref,tau,aerosol,ccn,rdust,rice,nuice)
+     &     tauref,tau,aerosol,ccn,rdust,rice,nuice,co2ice)
 
        IMPLICIT NONE
@@ -178 +178 @@
       REAL rice(ngridmx,nlayermx)   ! Ice geometric mean radius (m)
       REAL nuice(ngridmx,nlayermx)  ! Estimated effective variance
+      REAL co2ice(ngridmx)           ! co2 ice surface layer (kg.m-2)
+
 c
 c    Local variables :
@@ -491 +493 @@
      .        ,fluxsurf_lw(ig0+1),fluxtop_lw(ig0+1)
      .        ,znetrad
-     &        ,zQIRsQREF3d,zomegaIR3d,zgIR3d)
+     &        ,zQIRsQREF3d,zomegaIR3d,zgIR3d
+     &        ,co2ice(ig0+1))
 
 c       Shortwave ("sw") radiative transfer (= solar radiation)

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

@@ -4 +4 @@
      .                   ,fluxground,fluxtop
      .                   ,netrad
-     &                   ,QIRsQREF3d,omegaIR3d,gIR3d)
+     &                   ,QIRsQREF3d,omegaIR3d,gIR3d
+     &                   ,co2ice)
 
 c----------------------------------------------------------------------
@@ -76 +77 @@
       real newcoolrate(ndlon,nflev) ! cooling rate (K/s) / with implicite scheme
 
+      REAL co2ice(ndlo2)           ! co2 ice surface layer (kg.m-2)
+      REAL emis_gaz(ndlo2)         ! emissivity for gaz computations
+
       integer jk,jkk,ja,jl
 
@@ -103 +107 @@
       endif
 
+      DO jl=1 , kdlon
+         IF(co2ice(jl) .GT. 20.e-3) THEN
+             emis_gaz(jl)=1. 
+         ELSE
+             emis_gaz(jl)=emis(jl)
+         ENDIF
+      ENDDO
+
 c----------------------------------------------------------------------
 c         1.0   planck function
@@ -129 +141 @@
 
 c     print*, 'CALL of DISTANTS'
-      call lwxd ( ig0, kdlon, kflev, emis
+      call lwxd ( ig0, kdlon, kflev, emis_gaz
      .          , aer_t, co2_u, co2_up)
 
@@ -148 +160 @@
 
 c     print*, 'CALL of BOUNDARIES'
-      call lwxb ( ig0, kdlon, kflev, emis
+      call lwxb ( ig0, kdlon, kflev, emis_gaz
      .          , aer_t, co2_u, co2_up)
 

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

@@ -548 +548 @@
      $     emis,mu0,pplev,pplay,pt,tsurf,fract,dist_sol,igout,
      $     zdtlw,zdtsw,fluxsurf_lw,fluxsurf_sw,fluxtop_lw,fluxtop_sw,
-     &     tauref,tau,aerosol,ccn,rdust,rice,nuice)
+     &     tauref,tau,aerosol,ccn,rdust,rice,nuice,co2ice)
 
 c          Outputs for basic check (middle of domain)
@@ -1692 +1692 @@
 
 
-         z_out=1. !0.
+         z_out=0.
          if (calltherm .and. (z_out .gt. 0.)) then
          call surflayer_interpol(ngrid,nlayer,z0,g,zzlay,zu,zv,wmax_th