Changeset 2715 for trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbcana.F90

Timestamp:

2011-03-30T17:58:35+02:00 (13 years ago)

Author:

rblod

Message:

First attempt to put dynamic allocation on the trunk

File:

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

trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbcana.F90

@@ -27 +27 @@
    PUBLIC   sbc_gyre   ! routine called in sbcmod module
 
-   !                               !!* Namelist namsbc_ana *
-   INTEGER  ::   nn_tau000 = 1      ! nb of time-step during which the surface stress
-   !                                ! increase from 0 to its nominal value 
-   REAL(wp) ::   rn_utau0  = 0.e0   ! constant wind stress value in i-direction
-   REAL(wp) ::   rn_vtau0  = 0.e0   ! constant wind stress value in j-direction
-   REAL(wp) ::   rn_qns0   = 0.e0   ! non solar heat flux
-   REAL(wp) ::   rn_qsr0   = 0.e0   !     solar heat flux
-   REAL(wp) ::   rn_emp0   = 0.e0   ! net freshwater flux
+   !                                !!* Namelist namsbc_ana *
+   INTEGER  ::   nn_tau000 = 1       ! nb of time-step during which the surface stress
+   !                                 ! increase from 0 to its nominal value 
+   REAL(wp) ::   rn_utau0  = 0._wp   ! constant wind stress value in i-direction
+   REAL(wp) ::   rn_vtau0  = 0._wp   ! constant wind stress value in j-direction
+   REAL(wp) ::   rn_qns0   = 0._wp   ! non solar heat flux
+   REAL(wp) ::   rn_qsr0   = 0._wp   !     solar heat flux
+   REAL(wp) ::   rn_emp0   = 0._wp   ! net freshwater flux
    
    !! * Substitutions
@@ -42 +42 @@
    !! NEMO/OPA 3.3 , NEMO Consortium (2010)
    !! $Id$
-   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
-   !!----------------------------------------------------------------------
-
+   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
+   !!----------------------------------------------------------------------
 CONTAINS
 
@@ -64 +63 @@
       !!----------------------------------------------------------------------
       INTEGER, INTENT(in) ::   kt       ! ocean time step
-      !!
-      REAL(wp)            ::   zfacto                ! local scalar
-      REAL(wp)            ::   zrhoa  = 1.22         ! Air density kg/m3
-      REAL(wp)            ::   zcdrag = 1.5e-3       ! drag coefficient
-      REAL(wp)            ::   ztx, zty, zmod, zcoef ! temporary variables
+      !
+      REAL(wp) ::   zfacto                ! local scalar
+      REAL(wp) ::   zrhoa  = 1.22_wp      ! Air density kg/m3
+      REAL(wp) ::   zcdrag = 1.5e-3_wp    ! drag coefficient
+      REAL(wp) ::   ztx, zty, zmod, zcoef ! temporary variables
       !!
       NAMELIST/namsbc_ana/ nn_tau000, rn_utau0, rn_vtau0, rn_qns0, rn_qsr0, rn_emp0
@@ -75 +74 @@
       IF( kt == nit000 ) THEN
          !
-         REWIND ( numnam )                   ! Read Namelist namsbc : surface fluxes
-         READ   ( numnam, namsbc_ana )
+         REWIND( numnam )                    ! Read Namelist namsbc : surface fluxes
+         READ  ( numnam, namsbc_ana )
          !
          IF(lwp) WRITE(numout,*)' '
@@ -88 +87 @@
          IF(lwp) WRITE(numout,*)'              net heat flux          rn_emp0   = ', rn_emp0  , ' Kg/m2/s'
          !
-         nn_tau000 = MAX( nn_tau000, 1 )   ! must be >= 1
+         nn_tau000 = MAX( nn_tau000, 1 )     ! must be >= 1
          !
       ENDIF
 
-      qns   (:,:) = rn_qns0
-      qsr   (:,:) = rn_qsr0
-      emp   (:,:) = rn_emp0
-      emps  (:,:) = rn_emp0
+      qns (:,:) = rn_qns0
+      qsr (:,:) = rn_qsr0
+      emp (:,:) = rn_emp0
+      emps(:,:) = rn_emp0
    
       ! Increase the surface stress to its nominal value during the first nn_tau000 time-steps
@@ -305 +304 @@
          WRITE(numout,*)'           ndastp     = ', ndastp
          WRITE(numout,*)'           adatrj     = ', adatrj
-
       ENDIF
       !