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

Timestamp:

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

Author:

smasson

Message:

supress useless variables in phycst, see ticket:602

File:

• trunk/NEMO/OPA_SRC/SBC/sbcana.F90 (modified) (6 diffs)

trunk/NEMO/OPA_SRC/SBC/sbcana.F90

@@ -146 +146 @@
       REAL(wp) ::   zcdrag = 1.5e-3       ! drag coefficient
       REAL(wp) ::   ztx, zty, zmod, zcoef ! temporary variables
+      REAL(wp) ::   zyydd                 ! number of days in one year
       !!---------------------------------------------------------------------
-         
+      zyydd = REAL(nyear_len(1),wp)
+
       ! ---------------------------- !
       !  heat and freshwater fluxes  !
@@ -160 +162 @@
       ! current day (in hours) since january the 1st of the current year
       ztime = REAL( kt ) * rdt / (rmmss * rhhmm)   &       !  total incrementation (in hours)
-         &      - (nyear  - 1) * rjjhh * raajj             !  minus years since beginning of experiment (in hours)
+         &      - (nyear  - 1) * rjjhh * zyydd             !  minus years since beginning of experiment (in hours)
 
       ztimemax1 = ((5.*30.)+21.)* 24.                      ! 21th june     at 24h in hours
-      ztimemin1 = ztimemax1 + rjjhh * raajj / 2            ! 21th december        in hours
+      ztimemin1 = ztimemax1 + rjjhh * zyydd / 2            ! 21th december        in hours
       ztimemax2 = ((6.*30.)+21.)* 24.                      ! 21th july     at 24h in hours
-      ztimemin2 = ztimemax2 - rjjhh * raajj / 2            ! 21th january         in hours
-      !                                                    ! NB: rjjhh * raajj / 4 = one seasonal cycle in hours
+      ztimemin2 = ztimemax2 - rjjhh * zyydd / 2            ! 21th january         in hours
+      !                                                    ! NB: rjjhh * zyydd / 4 = one seasonal cycle in hours
 
       ! amplitudes
@@ -243 +245 @@
       ! day (in hours) since january the 1st
       ztime = FLOAT( kt ) * rdt / (rmmss * rhhmm)  &  ! incrementation in hour
-         &     - (nyear - 1) * rjjhh * raajj          !  - nber of hours the precedent years
+         &     - (nyear - 1) * rjjhh * zyydd          !  - nber of hours the precedent years
       ztimemax = ((5.*30.)+21.)* 24.               ! 21th june     in hours
-      ztimemin = ztimemax + rjjhh * raajj / 2      ! 21th december in hours
-      !                                            ! NB: rjjhh * raajj / 4 = 1 seasonal cycle in hours
+      ztimemin = ztimemax + rjjhh * zyydd / 2      ! 21th december in hours
+      !                                            ! NB: rjjhh * zyydd / 4 = 1 seasonal cycle in hours
 
       ! mean intensity at 0.105 ; srqt(2) because projected with 45deg angle
@@ -287 +289 @@
          WRITE(numout,*)'           nmonth     = ', nmonth
          WRITE(numout,*)'           nday       = ', nday
-         WRITE(numout,*)'           nday_year  = ',nday_year
+         WRITE(numout,*)'           nday_year  = ', nday_year
          WRITE(numout,*)'           ztime      = ', ztime
+         WRITE(numout,*)'           ztimemax   = ', ztimemax
+         WRITE(numout,*)'           ztimemin   = ', ztimemin
          WRITE(numout,*)'           ztimemax1  = ', ztimemax1
          WRITE(numout,*)'           ztimemin1  = ', ztimemin1
@@ -297 +301 @@
          WRITE(numout,*)'           zday0      = ', zday0
          WRITE(numout,*)'           zday_year0 = ', zday_year0
-         WRITE(numout,*)'           raajj      = ', raajj
+         WRITE(numout,*)'           zyydd      = ', zyydd
          WRITE(numout,*)'           zemp_S     = ', zemp_S
          WRITE(numout,*)'           zemp_N     = ', zemp_N
@@ -306 +310 @@
          WRITE(numout,*)'           ztrp       = ', ztrp
          WRITE(numout,*)'           zconv      = ', zconv
-
-         WRITE(numout,*)'           ndastp     = ',ndastp
-         WRITE(numout,*)'           adatrj     = ',adatrj
-         WRITE(numout,*)'           ztime      = ',ztime
-
-         WRITE(numout,*)'           ztimemax   = ',ztimemax
-         WRITE(numout,*)'           ztimemin   = ',ztimemin
-         WRITE(numout,*)'           zyear0     = ', zyear0
-         WRITE(numout,*)'           zmonth0    = ', zmonth0
-         WRITE(numout,*)'           zday0      = ', zday0
-         WRITE(numout,*)'           zday_year0 = ',zday_year0
-         WRITE(numout,*)'           raajj  = ', raajj
+         WRITE(numout,*)'           ndastp     = ', ndastp
+         WRITE(numout,*)'           adatrj     = ', adatrj
+
       ENDIF
       !