Changeset 4192 for dynamico_lmdz/simple_physics

Timestamp:

Dec 20, 2019, 12:56:43 AM (5 years ago)

Author:

dubos

Message:

simple_physics : cleanup astronomy

Location:

dynamico_lmdz/simple_physics/phyparam

Files:

• param/iniphyparam.F (modified) (1 diff)
• param/orbite.F (deleted)
• param/phyparam.F (modified) (3 diffs)
• param/sw.F (modified) (1 diff)
• param/zerofi.F (deleted)
• physics/astronomy.F90 (modified) (2 diffs)

dynamico_lmdz/simple_physics/phyparam/param/iniphyparam.F

@@ -64 +64 @@
       INTEGER ig,ierr,offset
  
-      EXTERNAL inifrict,orbite
+      EXTERNAL inifrict
  
       print*,'INIPHYPARAM'

dynamico_lmdz/simple_physics/phyparam/param/phyparam.F

@@ -164 +164 @@
 
       EXTERNAL convadj
-      EXTERNAL orbite,mucorr
+      EXTERNAL mucorr
       EXTERNAL ismin,ismax
 
@@ -482 +482 @@
          ENDDO
 
-         CALL zerophys(ngrid*nlayer,zdum1)
-         CALL zerophys(ngrid*nlayer,zdum2)
-c        CALL zerophys(ngrid*nlayer,zdum3)
+         zdum1(:,:)=0.
+         zdum2(:,:)=0.
+
          do l=1,nlayer
             do ig=1,ngrid
@@ -546 +546 @@
             ENDDO
          ENDDO
-         CALL zerophys(ngrid*nlayer,zdufr)
-         CALL zerophys(ngrid*nlayer,zdvfr)
-         CALL zerophys(ngrid*nlayer,zdhfr)
+
+         zdufr(:,:)=0.
+         zdvfr(:,:)=0.
+         zdhfr(:,:)=0.
+
          CALL convadj(ngrid,nlayer,ptimestep,
      $                pplay,pplev,zpopsk,

dynamico_lmdz/simple_physics/phyparam/param/sw.F

@@ -204 +204 @@
 
       IF (ldiurn) THEN
-         CALL zerophys(ngrid,fsrfvis)
+         fsrfvis(:)=0.
          CALL monscatter(ncount,fsrfvis,index,zflux)
-         CALL zerophys(ngrid*nlayer,dtsw)
+         dtsw(:,:)=0.
          DO l=1,nlayer
             CALL monscatter(ncount,dtsw(1,l),index,zdtsw(1,l))

dynamico_lmdz/simple_physics/phyparam/physics/astronomy.F90

@@ -7 +7 @@
   REAL, PARAMETER :: unitastr=149.597927, & ! millions of km
        pi=2.*ASIN(1.)
+
 CONTAINS
  
@@ -166 +167 @@
   END SUBROUTINE iniorbit
 
+  PURE SUBROUTINE orbite(pls,pdist_sol,pdecli)
+    !=======================================================================
+    !
+    !   Objet:
+    !   ------
+    !
+    !   Distance from sun and declimation as a function of the solar
+    !   longitude Ls
+    !
+    !   Arguments:
+    !   ----------
+    !
+    !   Input:
+    !   ------
+    !   pls          Ls
+    !
+    !   Output:
+    !   -------
+    !   pdist_sol     Distance Sun-Planet in UA
+    !   pdecli        declinaison ( en radians )
+    !
+    !=======================================================================
+    !-----------------------------------------------------------------------
+    !   Declarations:
+    !   -------------
+    
+    ! arguments:
+    ! ----------
+
+    REAL, INTENT(IN) :: pls
+    REAL, INTENT(OUT) :: pdist_sol,pdecli
+    
+    !-----------------------------------------------------------------------
+    
+    ! Distance Sun-Planet
+    
+    pdist_sol=p_elips/(1.+e_elips*cos(pls+timeperi))
+    
+    ! Solar declination
+    
+    pdecli= asin (sin(pls)*sin(obliquit*pi/180.))
+    
+  END SUBROUTINE orbite
+
 END MODULE astronomy