Changeset 433 for trunk/NEMO/TOP_SRC/trcdta.F90

Timestamp:

2006-04-10T17:40:29+02:00 (18 years ago)

Author:

opalod

Message:

nemo_v1_update_044 : CT : update the passive tracers TOP component and the standard GYRE configuration

File:

• trunk/NEMO/TOP_SRC/trcdta.F90 (modified) (9 diffs)

trunk/NEMO/TOP_SRC/trcdta.F90

@@ -1 +1 @@
 MODULE trcdta
    !!======================================================================
-   !!                     ***  MODULE  dtatem  ***
-   !! Ocean data  :  read passive tracer data from monthly atlas data
+   !!                     ***  MODULE  trcdta  ***
+   !! Ocean data  :  reads passive tracer data 
    !!=====================================================================
    !!  TOP 1.0,  LOCEAN-IPSL (2005) 
@@ -29 +29 @@
    !! * Shared module variables
    REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk,jptra) ::   &  !:
-      trdta             !: temperature data at given time-step
+      trdta             !: tracer data at given time-step
 
    !! * Module variables
    REAL(wp), DIMENSION(jpi,jpj,jpk,jptra,2) ::   &
-      tracdta            ! temperature data at two consecutive times
+      tracdta            ! tracer data at two consecutive times
    INTEGER , DIMENSION(jptra) :: &
       nlectr  ,   &    !!: switch for reading once
-      ntrc1   ,   &    !!: ????
-      ntrc2            !!: ????
+      ntrc1   ,   &    !!: number of first month when reading 12 monthly value
+      ntrc2            !!: number of second month when reading 12 monthly value
 
    !! * Substitutions
@@ -99 +99 @@
                !! 3D tracer data
                IF(lwp)WRITE(numout,*)
-               IF(lwp)WRITE(numout,*) ' trcdta: reading tracer' 
-               IF(lwp)WRITE(numout,*) ' data file ', jn
+               IF(lwp)WRITE(numout,*) ' dta_trc: reading tracer' 
+               IF(lwp)WRITE(numout,*) ' data file ', jn, ctrcnm(jn)
                IF(lwp)WRITE(numout,*)
                nlectr(jn) = 0
@@ -112 +112 @@
             ipi = jpiglo
             ipj = jpjglo
+            ipk = jpk
 
             ! First call kt=nit000
@@ -120 +121 @@
                IF(lwp) THEN
                   WRITE(numout,*)
-                  WRITE(numout,*) ' Tracer monthly fields' 
+                  WRITE(numout,*) ' Tracer data  fields' 
                   WRITE(numout,*) ' ~~~~~~~~~~~~~~~~~~~~~'
                   WRITE(numout,*) ' NetCDF FORMAT'
@@ -127 +128 @@
 
                ! open file 
-
+#if defined key_trc_pisces
                clname(jn) = 'LEVITUS_'//ctrcnm(jn)
+#else
+               itime=1
+               clname(jn) = ctrcnm(jn)
+#endif
                CALL flinopen(TRIM(clname(jn)),mig(1),nlci,mjg(1),nlcj,    &
                   .FALSE.,ipi,ipj,ipk,zlon,zlat,zlev,itime,    &
                   istep,zdate0,rdt,numtr(jn)               )
 
+#if defined key_trc_pisces
                ! title, dimensions and tests
                IF( itime /= jpmois ) THEN
@@ -140 +146 @@
                      WRITE(numout,*) ' itime ',itime,' jpmois ',jpmois
                   ENDIF
-                  STOP 'trc_dta'
-               ENDIF
+                  STOP 'dta_trc'
+               ENDIF
+#endif
 
                IF( ipi /= jpidta .OR. ipj /= jpjdta .OR. ipk /= jpk ) THEN
@@ -151 +158 @@
                      WRITE(numout,*) ' ipk ',ipk,' jpk ',jpk
                   ENDIF
-                  STOP 'trc_dta'
-               ENDIF
-               IF(lwp)WRITE(numout,*) itime,istep,zdate0,rdt,numtr(jn)
+                  STOP 'dta_trc'
+               ENDIF
+               IF(lwp)WRITE(numout,*) itime,istep(1),zdate0,rdt,numtr(jn)
                trdta(:,:,:,jn) = 0.
 
             ENDIF
 
-
+#if defined key_trc_pisces
             ! Read montly file
             IF( ( kt == nit000 .AND. nlectr(jn) == 0)   & 
@@ -244 +251 @@
             IF( jn == jpsil) trdta(:,:,:,jn) = trdta(:,:,:,jn) * 1.E-6
             IF( jn == jppo4) trdta(:,:,:,jn) = trdta(:,:,:,jn) * 122.E-6
-
-
-         ENDIF
-
-      END DO
+#else
+            ! Read init file only
+            IF( kt == nit000  ) THEN
+               CALL flinget( numtr(jn),ctrcnm(jn),jpidta,jpjdta,jpk,    &
+               1,1,1,mig(1),nlci,mjg(1),nlcj,  &
+               trdta(1:nlci,1:nlcj,1:jpk,jn) )
+               trdta(:,:,:,jn)=trdta(:,:,:,jn)*tmask(:,:,:)
+            ENDIF 
+#endif
+
+        ENDIF
+
+       END DO
 
    END SUBROUTINE dta_trc