Changeset 775 for branches/dev_001_GM/NEMO/TOP_SRC/PISCES_SMS/p4zday.F90

Timestamp:

2007-12-19T14:45:15+01:00 (16 years ago)

Author:

gm

Message:

dev_001_GM - PISCES in F90 : encapsulation of all p4z...F files in module F90 + doctor norme for local variables - compilation OK

File:

• branches/dev_001_GM/NEMO/TOP_SRC/PISCES_SMS/p4zday.F90 (moved) (moved from branches/dev_001_GM/NEMO/TOP_SRC/PISCES_SMS/p4zday.F) (1 diff)

branches/dev_001_GM/NEMO/TOP_SRC/PISCES_SMS/p4zday.F90

@@ -1 @@
-CCC$Header$
-CCC  TOP 1.0 , LOCEAN-IPSL (2005)
-C This software is governed by CeCILL licence see modipsl/doc/NEMO_CeCILL.txt
-C ---------------------------------------------------------------------------
-CDIR$ LIST
-      SUBROUTINE p4zday
-#if defined key_top && defined key_pisces
-CCC---------------------------------------------------------------------
-CCC
-CCC          ROUTINE p4zday : PISCES MODEL
-CCC          *****************************
-CCC
-CCC  PURPOSE :
-CCC  ---------
-CCC        PISCES : compute the day length depending on latitude
-CCC                 and the day
-CCC
-CC   INPUT :
-CC   -----
-CC      argument
-CC              ktask           : task identificator
-CC      common
-CC              all the common defined in opa
-CC
-CC
-CC   OUTPUT :                   : no
-CC   ------
-CC
-CC   EXTERNAL :
-CC   --------
-CC            None
-CC
-CC   MODIFICATIONS:
-CC   --------------
-CC      original  : E. Maier-Reimer (GBC 1993)
-CC      additions : C. Le Quere (1999)
-CC      modifications : O. Aumont (2004)
-CC----------------------------------------------------------------------
-CC parameters and commons
-CC ======================
-      USE oce_trc
-      USE trp_trc
-      USE sms
-      IMPLICIT NONE
-CC----------------------------------------------------------------------
-CC local declarations
-CC ==================
-      INTEGER ji, jj, iyy
-      REAL rum, delta, codel, phi, argu
-C
-C Get year
-C --------
-C
-      iyy = ndastp/10000
+MODULE p4zday
+   !!======================================================================
+   !!                         ***  MODULE p4zday  ***
+   !! TOP :   PISCES compute the day length depending on latitude and the day
+   !!======================================================================
+   !! History :    -   !  1993     (E. Maier-Reimer) Original code GBC 1993
+   !!              -   !  1999     (C. Le Quere)
+   !!             1.0  !  2004     (O. Aumont) Original code
+   !!             2.0  !  2007-12  (C. Ethe, G. Madec)  F90
+   !!----------------------------------------------------------------------
+#if defined key_pisces
+   !!----------------------------------------------------------------------
+   !!   'key_pisces'                                       PISCES bio-model
+   !!----------------------------------------------------------------------
+   !!   p4z_day       :   compute the day length depending on latitude and the day
+   !!----------------------------------------------------------------------
+   USE oce_trc         !
+   USE trp_trc         ! 
+   USE sms             ! 
 
-       IF(lwp) write(numout,*)
-       IF(lwp) write(numout,*) 'p4zday - Julian day ', nday_year
-       IF(lwp) write(numout,*)
+   IMPLICIT NONE
+   PRIVATE
+
+   PUBLIC   p4z_day    ! called in p4zprod.F90
+
+   !!* Substitution
+#  include "domzgr_substitute.h90"
+   !!----------------------------------------------------------------------
+   !! NEMO/TOP 2.0 , LOCEAN-IPSL (2007) 
+   !! $Header:$ 
+   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
+   !!----------------------------------------------------------------------
+
+CONTAINS
+
+   SUBROUTINE p4z_day
+      !!---------------------------------------------------------------------
+      !!                     ***  ROUTINE p4z_day  ***
+      !!
+      !! ** Purpose :   compute the day length depending on latitude and the day
+      !!
+      !! ** Method  : - ???
+      !!---------------------------------------------------------------------
+      INTEGER  ::   ji, jj
+      INTEGER  ::   iyy
+      REAL(wp) ::   zrum, zdelta, zcodel, zphi, zargu
+      !!---------------------------------------------------------------------
+
+      ! Get year
+      ! --------
+
+      iyy = ndastp / 10000
+
+      IF(lwp) write(numout,*)
+      IF(lwp) write(numout,*) 'p4zday : - Julian day ', nday_year
+      IF(lwp) write(numout,*) '~~~~~~'
 
 
-      IF (nleapy.EQ.1 .AND. MOD(iyy,4).EQ.0) THEN
-          rum = FLOAT(nday_year-80)/366.
+      IF( nleapy == 1 .AND. MOD( iyy, 4 ) == 0 ) THEN
+         zrum = FLOAT( nday_year - 80 ) / 366.
       ELSE
-          rum = FLOAT(nday_year-80)/365.
+         zrum = FLOAT( nday_year - 80 ) / 365.
       ENDIF
 
-C
-      delta = SIN(rum*rpi*2.)*sin(rpi*23.5/180.)
-      codel = asin(delta)
+      zdelta = SIN( zrum * rpi * 2. ) * sin( rpi * 23.5 / 180. )
+      zcodel = ASIN( zdelta )
 
-      DO jj = 1,jpj
-        DO ji = 1,jpi
-          phi = gphit(ji,jj)*rpi/180.
-          argu = tan(codel)*tan(phi)
-          strn(ji,jj) = 0.
-          argu=min(1.,argu)
-          argu=max(-1.,argu)
-          strn(ji,jj)=24.-2.*acos(argu)*180./rpi/15.
-          strn(ji,jj)=max(strn(ji,jj),0.)
-        END DO
+      DO jj = 1, jpj
+         DO ji = 1, jpi
+            zphi = gphit(ji,jj) * rpi / 180.
+            zargu = TAN( zcodel ) * TAN( zphi )
+            strn(ji,jj) = 0.e0
+            zargu = MIN(  1., zargu )
+            zargu = MAX( -1., zargu )
+            strn(ji,jj) = 24.- 2.* ACOS( zargu ) * 180./ rpi / 15.
+            strn(ji,jj) = MAX( strn(ji,jj), 0.e0 )
+         END DO
       END DO
-C
-#endif
-      RETURN
-      END
+      !
+   END SUBROUTINE p4z_day
+
+#else
+   !!======================================================================
+   !!  Dummy module :                                   No PISCES bio-model
+   !!======================================================================
+CONTAINS
+   SUBROUTINE p4z_day                    ! Empty routine
+   END SUBROUTINE p4z_day
+#endif 
+
+   !!======================================================================
+END MODULE  p4zday