Changeset 167 for trunk/NEMO/OPA_SRC/SBC/ocesbc.F90

Timestamp:

2004-11-05T14:39:09+01:00 (20 years ago)

Author:

opalod

Message:

CT : BUGFIX111 : remove all special characters and lines too long for the GYRE configuration part

File:

• trunk/NEMO/OPA_SRC/SBC/ocesbc.F90 (modified) (1 diff)

trunk/NEMO/OPA_SRC/SBC/ocesbc.F90

@@ -655 +655 @@
           ztimemax = zdaymax * 24
           ! 21th december day in hours
-          ztimemin = ztimemax + rjjhh * raajj / 2         ! rjjhh * raajj / 4 = 1 seasonal cycle in hours
+          ! rjjhh * raajj / 4 = 1 seasonal cycle in hours
+          ztimemin = ztimemax + rjjhh * raajj / 2         
           ! amplitudes
           zemp_S   = 0.7       ! intensity of COS in the South
           zemp_N   = 0.8       ! intensity of COS in the North
           zemp_sais= 0.1
-          zTstar   = 28.3      ! intemsity from 28.3 à -5°
-          zcos_sais = COS( (ztime - ztimemax) / (ztimemin - ztimemax) * rpi )    ! 1/2 period between 21th June and 21th December
+          zTstar   = 28.3      ! intemsity from 28.3 a -5 deg
+          ! 1/2 period between 21th June and 21th December
+          zcos_sais = COS( (ztime - ztimemax) / (ztimemin - ztimemax) * rpi )    
           ztrp= - 40.          ! retroaction term (W/m2/K)
           zconv = 3.16e-5      ! convert 1m/yr->3.16e-5mm/s
           DO jj = 1, jpj
-            DO ji = 1, jpi
-              t_star (ji,jj) = zTstar * ( 1 + 1. / 50. * zcos_sais )   &                  ! domain from 15° to 50°
-                 &             * COS( rpi * (gphit(ji,jj) - 5.)        &                  ! between 27 and 28 °C at 15N, -3 
-                 &                        / (53.5 * ( 1 + 11 / 53.5 * zcos_sais ) * 2.) ) ! and 13°C at 50N 53.5 + or - 11 
-                 !                                                                        ! = 1/4 period :
-                 !                                                                        ! 64.5 in summer, 42.5 in winter
-              qt  (ji,jj) = ztrp * ( tb(ji,jj,1) - t_star(ji,jj) )
-              IF( gphit(ji,jj) >= 14.845 .AND. 37.2 >= gphit(ji,jj)) THEN
-                 emp  (ji,jj) =   zemp_S * zconv   &
-                    &         * SIN( rpi / 2 * (gphit(ji,jj) - 37.2) / (24.6 - 37.2) )  & ! zero at 37.8°, max at 24.6°
-                    &         * ( 1 - zemp_sais / zemp_S * zcos_sais)
-                 emps (ji,jj) = emp  (ji,jj)
-              ELSE
-                 emp (ji,jj) =  - zemp_N * zconv   &
-                    &         * SIN( rpi / 2 * (gphit(ji,jj) - 37.2) / (46.8 - 37.2) )  & ! zero at 37.8°, max at 46.8°
-                    &         * ( 1 - zemp_sais / zemp_N * zcos_sais )
-                  emps (ji,jj) = emp  (ji,jj)
-              ENDIF
-              qsr (ji,jj) =  230 * COS( 3.1415 * ( gphit(ji,jj) - 23.5 * zcos_sais ) / ( 0.9 * 180 ) )   ! 23.5° : tropics
-            END DO
+             DO ji = 1, jpi
+                ! domain from 15 deg to 50 deg between 27 and 28  degC at 15N, -3 
+                ! and 13 degC at 50N 53.5 + or - 11 = 1/4 period :
+                ! 64.5 in summer, 42.5 in winter
+                t_star (ji,jj) = zTstar * ( 1 + 1. / 50. * zcos_sais )   &                  
+                   &             * COS( rpi * (gphit(ji,jj) - 5.)        &                  
+                   &                        / (53.5 * ( 1 + 11 / 53.5 * zcos_sais ) * 2.) ) 
+                qt  (ji,jj) = ztrp * ( tb(ji,jj,1) - t_star(ji,jj) )
+                IF( gphit(ji,jj) >= 14.845 .AND. 37.2 >= gphit(ji,jj)) THEN
+                   ! zero at 37.8 deg, max at 24.6 deg
+                   emp  (ji,jj) =   zemp_S * zconv   &
+                      &         * SIN( rpi / 2 * (gphit(ji,jj) - 37.2) / (24.6 - 37.2) )  & 
+                      &         * ( 1 - zemp_sais / zemp_S * zcos_sais)
+                   emps (ji,jj) = emp  (ji,jj)
+                ELSE
+                   ! zero at 37.8 deg, max at 46.8 deg
+                   emp (ji,jj) =  - zemp_N * zconv   &
+                      &         * SIN( rpi / 2 * (gphit(ji,jj) - 37.2) / (46.8 - 37.2) )  & 
+                      &         * ( 1 - zemp_sais / zemp_N * zcos_sais )
+                    emps (ji,jj) = emp  (ji,jj)
+                ENDIF
+                ! 23.5 deg : tropics
+                qsr (ji,jj) =  230 * COS( 3.1415 * ( gphit(ji,jj) - 23.5 * zcos_sais ) / ( 0.9 * 180 ) )   
+             END DO
           END DO
           ! compute the emp flux such as its integration on the whole domain and at each time be zero