Changeset 1732 for trunk/NEMO/OPA_SRC/SBC/sbcblk_clio.F90

Timestamp:

2009-11-16T16:15:24+01:00 (14 years ago)

Author:

smasson

Message:

supress useless variables in phycst, see ticket:602

File:

• trunk/NEMO/OPA_SRC/SBC/sbcblk_clio.F90 (modified) (7 diffs)

trunk/NEMO/OPA_SRC/SBC/sbcblk_clio.F90

@@ -63 +63 @@
 #endif
 
-   REAL(wp) ::   yearday     !: number of days per year   
-   REAL(wp) ::   rdtbs2      !: number of days per year   
+   REAL(wp) ::   rdtbs2      !:   
    
    REAL(wp), DIMENSION(19)  ::  budyko            ! BUDYKO's coefficient (cloudiness effect on LW radiation)
@@ -725 +724 @@
             END DO
          END DO
-         IF    ( nleapy == 1 ) THEN   ;   yearday = 366.e0
-         ELSEIF( nleapy == 0 ) THEN   ;   yearday = 365.e0
-         ELSEIF( nleapy == 30) THEN   ;   yearday = 360.e0
-         ENDIF
          lbulk_init = .FALSE.
       ENDIF
@@ -764 +759 @@
       !  correction factor added for computation of shortwave flux to take into account the variation of
       !  the distance between the sun and the earth during the year (Oberhuber 1988)
-      zdist    = zxday * 2. * rpi / yearday
+      zdist    = zxday * 2. * rpi / REAL(nyear_len(1), wp)
       zdaycor  = 1.0 + 0.0013 * SIN( zdist ) + 0.0342 * COS( zdist )
 
@@ -779 +774 @@
             zlsset (ji,jj) = - zlsrise(ji,jj)
             !  dividing the solar day into jp24 segments of length zdlha
-            zdlha  (ji,jj) = ( zlsrise(ji,jj) - zlsset(ji,jj) ) / REAL( jp24 )
+            zdlha  (ji,jj) = ( zlsrise(ji,jj) - zlsset(ji,jj) ) / REAL( jp24, wp )
          END DO
       END DO
@@ -894 +889 @@
       !  correction factor added for computation of shortwave flux to take into account the variation of
       !  the distance between the sun and the earth during the year (Oberhuber 1988)
-      zdist    = zxday * 2. * rpi / yearday
+      zdist    = zxday * 2. * rpi / REAL(nyear_len(1), wp)
       zdaycor  = 1.0 + 0.0013 * SIN( zdist ) + 0.0342 * COS( zdist )
 
@@ -909 +904 @@
             zlsset (ji,jj) = - zlsrise(ji,jj)
             !  dividing the solar day into jp24 segments of length zdlha
-            zdlha  (ji,jj) = ( zlsrise(ji,jj) - zlsset(ji,jj) ) / REAL( jp24 )
+            zdlha  (ji,jj) = ( zlsrise(ji,jj) - zlsset(ji,jj) ) / REAL( jp24, wp )
          END DO
       END DO
@@ -991 +986 @@
       !!---------------------------------------------------------------------
             
-      IF    ( ky == 1 )  THEN   ;   zday = REAL( kday ) - 0.5
-      ELSEIF( ky == 3 )  THEN   ;   zday = REAL( kday ) - 1.
-      ELSE                      ;   zday = REAL( kday )
+      IF    ( ky == 1 )  THEN   ;   zday = REAL( kday, wp ) - 0.5
+      ELSEIF( ky == 3 )  THEN   ;   zday = REAL( kday, wp ) - 1.
+      ELSE                      ;   zday = REAL( kday, wp )
       ENDIF
       
-      zp = rpi * ( 2.0 * zday - 367.0 ) / yearday
+      zp = rpi * ( 2.0 * zday - 367.0 ) / REAL(nyear_len(1), wp)
       
       pdecl  = a0                                                                      &