Changeset 4191 for dynamico_lmdz/simple_physics/phyparam/param/zenang.F

Timestamp:

Dec 20, 2019, 12:41:10 AM (5 years ago)

Author:

dubos

Message:

simple_physics : cleanup astronomy

File:

• dynamico_lmdz/simple_physics/phyparam/param/zenang.F (modified) (3 diffs)

dynamico_lmdz/simple_physics/phyparam/param/zenang.F

@@ -1 +1 @@
       SUBROUTINE zenang(klon,longi,gmtime,pdtrad,lat,long,
      s                  pmu0,frac)
+      USE astronomy
       IMPLICIT none
 c=============================================================
@@ -24 +25 @@
 c frac-----OUTPUT: ensoleillement moyen entre gmtime et gmtime+pdtrad
 c================================================================
-#include "comorbit.h"
       integer klon
 c================================================================
@@ -155 +155 @@
       END
 c===================================================================
-#ifdef PASLA
-      SUBROUTINE zenith (klon,longi, gmtime, lat, long,
-     s                   pmu0, fract)
-      IMPLICIT none
-c
-c Auteur(s): Z.X. Li (LMD/ENS)
-c
-c Objet: calculer le cosinus de l'angle zenithal du soleil en
-c        connaissant la declinaison du soleil, la latitude et la
-c        longitude du point sur la terre, et le temps universel
-c
-c Arguments d'entree:
-c     longi  : declinaison du soleil (en degres)
-c     gmtime : temps universel en second qui varie entre 0 et 86400
-c     lat    : latitude en degres
-c     long   : longitude en degres
-c Arguments de sortie:
-c     pmu0   : cosinus de l'angle zenithal
-c
-c====================================================================
-#include "YOMCST.h"
-c====================================================================
-      INTETER klon
-      REAL longi, gmtime
-      REAL lat(klon), long(klon), pmu0(klon), fract(klon)
-c=====================================================================
-      INTEGER n
-      REAL zpi, zpir, omega, zgmtime
-      REAL incl, lat_sun, lon_sun
-c----------------------------------------------------------------------
-      zpi = 4.0*ATAN(1.0)
-      zpir = zpi / 180.0
-      zgmtime=gmtime*86400.
-c
-      incl=R_incl * zpir
-c
-      lon_sun = longi * zpir
-      lat_sun = ASIN (SIN(lon_sun)*SIN(incl) )
-c
-c--initialisation a la nuit
-c
-      DO n =1, klon
-        pmu0(n)=0.
-        fract(n)=0.0
-      ENDDO
-c
-c 1 degre en longitude = 240 secondes en temps
-c
-      DO n = 1, klon
-         omega = zgmtime + long(n)*86400.0/360.0
-         omega = omega / 86400.0 * 2.0 * zpi
-         omega = MOD(omega + 2.0 * zpi, 2.0 * zpi)
-         omega = omega - zpi
-         pmu0(n) = sin(lat(n)*zpir) * sin(lat_sun)
-     .           + cos(lat(n)*zpir) * cos(lat_sun)
-     .           * cos(omega)
-         pmu0(n) = MAX (pmu0(n), 0.0)
-         IF (pmu0(n).GT.1.E-6) fract(n)=1.0
-      ENDDO
-c
-      RETURN
-      END
-#endif