Changeset 1133 for trunk/LMDZ.GENERIC

Timestamp:

Dec 27, 2013, 1:53:08 AM (11 years ago)

Author:

aslmd

Message:

LMDZ.GENERIC. LMDZ.UNIVERSAL. A first version of ring shadowing subroutine. Authors: M. Sylvestre, M. Capderou, S. Guerlet, A. Spiga.

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• callcorrk.F90 (modified) (1 diff)
• callkeys.h (modified) (2 diffs)
• inifis.F (modified) (1 diff)
• physiq.F90 (modified) (5 diffs)
• radcommon_h.F90 (modified) (1 diff)
• rings.F90 (added)

trunk/LMDZ.GENERIC/libf/phystd/callcorrk.F90

@@ -697 +697 @@
             else
                do nw=1,L_NSPECTV
-                  stel_fract(nw)= stel(nw) * fract(ig)
+                  stel_fract(nw)= stel(nw) * fract(ig)*(1.-eclipse(ig))
                end do
             endif

trunk/LMDZ.GENERIC/libf/phystd/callkeys.h

@@ -30 +30 @@
      &   , co2supsat                                                    &
      &   , cloudlvl                                                     &
-     &   , pceil                                                    
+     &   , pceil                                                        &
+     &   , rings_shadow                                         
 
       logical callrad,corrk,calldifv,UseTurbDiff                        &
@@ -36 +37 @@
      &   , season,diurnal,tlocked,lwrite                                &
      &   , callstats,calleofdump                                        &
-     &   , callgasvis,continuum,H2Ocont_simple,graybody
+     &   , callgasvis,continuum,H2Ocont_simple,graybody                 &
+     &   , rings_shadow
 
       logical enertest

trunk/LMDZ.GENERIC/libf/phystd/inifis.F

@@ -153 +153 @@
          write(*,*) "tlocked = ",tlocked
 
+         write(*,*) "Saturn ring shadowing ?"
+         rings_shadow = .false.
+         call getin("rings_shadow", rings_shadow)
+         write(*,*) "rings_shadow = ", rings_shadow
+         
 ! Test of incompatibility:
 ! if tlocked, then diurnal should be false

trunk/LMDZ.GENERIC/libf/phystd/physiq.F90

@@ -11 +11 @@
       use watercommon_h, only : RLVTT, Psat_water,epsi
       use gases_h
-      use radcommon_h, only: sigma
+      use radcommon_h, only: sigma, eclipse
       use radii_mod, only: h2o_reffrad, co2_reffrad, h2o_cloudrad
       use aerosol_mod
@@ -112 +112 @@
 !           Loops converted to F90 matrix format: J. Leconte (2012)
 !           No more ngridmx/nqmx, F90 commons and adaptation to parallel: A. Spiga (2012)
-!     
 !==================================================================
 
@@ -460 +459 @@
          ALLOCATE(fract(ngrid))
 
+         !! this is defined in radcommon_h
+         ALLOCATE(eclipse(ngrid))
+
 !        variables set to 0
 !        ~~~~~~~~~~~~~~~~~~
@@ -700 +702 @@
                ! WARNING: this function appears not to work in 1D
 
+           endif
+
+           !! Eclipse incoming sunlight (e.g. Saturn ring shadowing)
+           if (rings_shadow) then
+             if (diurnal) then
+               call rings(ngrid, declin, ptime, rad, 0., eclipse)
+               call writediagfi(ngrid,"shad","rings"," ", 2, 1.-eclipse)
+               write(*,*) "Ring shadow activated."
+             else
+               write(*,*) "Ring shadow + diurnal=F not supported yet."
+             endif
+           else
+               eclipse(:) = 0.0
            endif
 
@@ -1891 +1906 @@
         IF ( ALLOCATED(sinlon)) DEALLOCATE(sinlon)
         IF ( ALLOCATED(coslon)) DEALLOCATE(coslon)
-
       endif
 

trunk/LMDZ.GENERIC/libf/phystd/radcommon_h.F90

@@ -135 +135 @@
       save Cmk
 
+      ! extinction of incoming sunlight (Saturn's rings, eclipses, etc...)
+      REAL, DIMENSION(:), ALLOCATABLE :: eclipse
 
       end module radcommon_h