Changeset 2359 for LMDZ5/trunk/libf/phylmd/orbite.F90

Timestamp:

Sep 4, 2015, 6:48:31 PM (9 years ago)

Author:

oboucher

Message:

Recoded diurnal cycle in physiq.F90
iflag_diurnal_cycle is the new flag
=0 for no diurnal cycle
=1 for old diurnal cycle
=2 for new diurnal cycle with updates
in between two calls to radiation according
to how rmu0 has changed over the previous
timestep of physiq

File:

• LMDZ5/trunk/libf/phylmd/orbite.F90 (modified) (5 diffs)

LMDZ5/trunk/libf/phylmd/orbite.F90

@@ -102 +102 @@
 END SUBROUTINE angle
 ! ====================================================================
-SUBROUTINE zenang(longi, gmtime, pdtrad, lat, long, pmu0, frac)
+SUBROUTINE zenang(longi, gmtime, pdtrad1, pdtrad2, lat, long, pmu0, frac)
   USE dimphy
   IMPLICIT NONE
@@ -114 +114 @@
   ! fournit des moyennes de pmu0 et non des valeurs
   ! instantanees, du coup frac prend toutes les valeurs
-  ! entre 0 et 1.
+  ! entre 0 et 1. La routine integre entre gmtime+pdtrad1 et 
+  ! gmtime+pdtrad2 avec pdtrad1 et pdtrad2 exprimes en secondes.
   ! Date   : premiere version le 13 decembre 1994
   ! revu pour  GCM  le 30 septembre 1996
+  ! revu le 3 septembre 2015 pour les bornes de l'integrale
   ! ===============================================================
   ! longi : la longitude vraie de la terre dans son plan
   ! solaire a partir de l'equinoxe de printemps (degre)
   ! gmtime : temps universel en fraction de jour
-  ! pdtrad : pas de temps du rayonnement (secondes)
+  ! pdtrad1 : borne inferieure du pas de temps du rayonnement (secondes)
+  ! pdtrad2 : borne inferieure du pas de temps du rayonnement (secondes)
+  ! pdtrad2-pdtrad1 correspond a pdtrad, le pas de temps du rayonnement (secondes)
   ! lat------INPUT : latitude en degres
   ! long-----INPUT : longitude en degres
-  ! pmu0-----OUTPUT: angle zenithal moyen entre gmtime et gmtime+pdtrad
-  ! frac-----OUTPUT: ensoleillement moyen entre gmtime et gmtime+pdtrad
-  ! ================================================================
-  include "YOMCST.h"
-  ! ================================================================
-  REAL, INTENT (IN) :: longi, gmtime, pdtrad
+  ! pmu0-----OUTPUT: angle zenithal moyen entre gmtime+pdtrad1 et gmtime+pdtrad2
+  ! frac-----OUTPUT: ensoleillement moyen entre gmtime+pdtrad1 et gmtime+pdtrad2
+  ! ================================================================
+  include "YOMCST.h"
+  ! ================================================================
+  REAL, INTENT (IN) :: longi, gmtime, pdtrad1, pdtrad2
   REAL lat(klon), long(klon), pmu0(klon), frac(klon)
   ! ================================================================
@@ -143 +147 @@
   REAL lat_sun ! declinaison en radian
   REAL lon_sun ! longitude solaire en radian
-  REAL latr ! latitude du pt de grille en radian
+  REAL latr    ! latitude du pt de grille en radian
   ! ================================================================
 
@@ -153 +157 @@
   lat_sun = asin(sin(lon_sun)*sin(incl))
 
-  gmtime1 = gmtime*86400.
-  gmtime2 = gmtime*86400. + pdtrad
+  gmtime1 = gmtime*86400. + pdtrad1
+  gmtime2 = gmtime*86400. + pdtrad2
 
   DO i = 1, klon
@@ -160 +164 @@
     latr = lat(i)*pi_local/180.
 
-    ! --pose probleme quand lat=+/-90 degres
-
-    ! omega = -TAN(latr)*TAN(lat_sun)
-    ! omega = ACOS(omega)
-    ! IF (latr.GE.(pi_local/2.+lat_sun)
-    ! .    .OR. latr.LE.(-pi_local/2.+lat_sun)) THEN
-    ! omega = 0.0       ! nuit polaire
-    ! ENDIF
-    ! IF (latr.GE.(pi_local/2.-lat_sun)
-    ! .          .OR. latr.LE.(-pi_local/2.-lat_sun)) THEN
-    ! omega = pi_local  ! journee polaire
-    ! ENDIF
-
-    ! --remplace par cela (le cas par defaut est different)
-
     omega = 0.0 !--nuit polaire
+
     IF (latr>=(pi_local/2.-lat_sun) .OR. latr<=(-pi_local/2.-lat_sun)) THEN
       omega = pi_local ! journee polaire
     END IF
+
     IF (latr<(pi_local/2.+lat_sun) .AND. latr>(-pi_local/2.+lat_sun) .AND. &
         latr<(pi_local/2.-lat_sun) .AND. latr>(-pi_local/2.-lat_sun)) THEN