Changeset 768 for branches/dev_001_GM/NEMO/TOP_SRC/CFC/trcini_cfc.F90

Timestamp:

2007-12-16T15:46:18+01:00 (16 years ago)

Author:

gm

Message:

dev_001_GM - create 1 trclsm_ module by trc model (CFC, LOBSTER, PISCES..) + some bug corrections

File:

• branches/dev_001_GM/NEMO/TOP_SRC/CFC/trcini_cfc.F90 (modified) (7 diffs)

branches/dev_001_GM/NEMO/TOP_SRC/CFC/trcini_cfc.F90

@@ -12 +12 @@
    !! trc_ini_cfc      : CFC model initialisation
    !!----------------------------------------------------------------------
+   USE oce_trc         ! Ocean variables
    USE par_trc         ! TOP parameters
+   USE trc             ! TOP variables
    USE trccfc          ! CFC sms trends
 
@@ -28 +30 @@
    !!----------------------------------------------------------------------
    !! NEMO/TOP 2.0 , LOCEAN-IPSL (2007) 
-   !! $Id:$ 
+   !! $Id$ 
    !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
    !!----------------------------------------------------------------------
@@ -52 +54 @@
       ! Initialization of boundaries conditions
       ! --------------------------------------- 
-      pp_cfc(:,:,:) = 0.e0
       qtr   (:,:,:) = 0.e0
       xphem (:,:)   = 0.e0
-      DO jn = 1, jptra
+      DO jn = jp_cfc0, jp_cfc1
          DO jm = 1, jphem
             DO jl = 1, jpyear
@@ -71 +72 @@
             WRITE(numout,*) 'Initialization de qint ; No restart : qint equal zero '
          ENDIF
-         DO jn = 1, jptra
+         DO jn = jp_cfc0, jp_cfc1
             qint(:,:,jn) = 0.e0
          END DO
@@ -88 +89 @@
       REWIND(inum)
       
-
-      ! Skip over 1st six descriptor lines
-      !-----------------------------------
-      DO jm = 1, 6
+      DO jm = 1, 6        ! Skip over 1st six descriptor lines
          READ(inum,'(1x)')
       END DO
    
-   
-      !   Read file
-      ! ---------
-      DO jn = 31, 98
+      DO jn = 31, 98      !   Read file
          READ(inum,*) zyy, p_cfc(jn,1,jp11), p_cfc(jn,1,jp12), &
             &              p_cfc(jn,2,jp11), p_cfc(jn,2,jp12)
@@ -106 +101 @@
       END DO
 
-      p_cfc(32,1:2,jp11) = 5.e-4
+      p_cfc(32,1:2,jp11) = 5.e-4      ! modify the values of the first years
       p_cfc(33,1:2,jp11) = 8.e-4
       p_cfc(34,1:2,jp11) = 1.e-6
@@ -115 +110 @@
       p_cfc(39,1:2,jp11) = 1.e-2
       
-      
-      IF(lwp) THEN
+      IF(lwp) THEN        ! Control print
          WRITE(numout,*)
          WRITE(numout,*) ' Year   p11HN    p11HS    p12HN    p12HS '