Changeset 1255 for trunk/NEMO/TOP_SRC/CFC/trclsm_cfc.F90

Timestamp:

2009-01-13T11:20:17+01:00 (16 years ago)

Author:

cetlod

Message:

minor modifications in all top models, see ticket:299

File:

• trunk/NEMO/TOP_SRC/CFC/trclsm_cfc.F90 (modified) (3 diffs)

trunk/NEMO/TOP_SRC/CFC/trclsm_cfc.F90

@@ -42 +42 @@
       !!----------------------------------------------------------------------
       CHARACTER (len=32) ::   clname = 'namelist_cfc'
-      INTEGER ::   numnat
+      INTEGER ::   numnatc
+#if defined key_trc_diaadd
+      ! definition of additional diagnostic as a structure
+      INTEGER :: jl, jn
+      TYPE DIAG
+         CHARACTER(len = 20)  :: snamedia   !: short name
+         CHARACTER(len = 80 ) :: lnamedia   !: long name
+         CHARACTER(len = 20 ) :: unitdia    !: unit
+      END TYPE DIAG
+
+      TYPE(DIAG) , DIMENSION(jp_cfc_2d) :: cfcdia2d
+#endif
       !!
       NAMELIST/namcfcdate/ ndate_beg, nyear_res
+#if defined key_trc_diaadd
+      NAMELIST/namcfcdia/nwritedia, cfcdia2d     ! additional diagnostics
+#endif
       !!-------------------------------------------------------------------
 
@@ -51 +65 @@
 
       !                             ! Open namelist file
-      CALL ctlopn( numnat, clname, 'OLD', 'FORMATTED', 'SEQUENTIAL',   &
-         &           1, numout, .FALSE., 1 )
+      CALL ctlopn( numnatc, clname, 'OLD', 'FORMATTED', 'SEQUENTIAL',  1, numout, .FALSE., 1 )
          
-      READ( numnat , namcfcdate )     ! read namelist
+      READ( numnatc , namcfcdate )     ! read namelist
 
       IF(lwp) THEN                  ! control print
@@ -66 +79 @@
       IF(lwp) WRITE(numout,*) '    initial year (aa)                       nyear_beg = ', nyear_beg
       !
+#if defined key_trc_diaadd
+
+      ! Namelist namcfcdia
+      ! -------------------
+      nwritedia = 10                   ! default values
+
+      DO jl = 1, jp_cfc_2d
+         jn = jp_cfc0_2d + jl - 1 
+         WRITE(ctrc2d(jn),'("2D_",I1)') jn                      ! short name
+         WRITE(ctrc2l(jn),'("2D DIAGNOSTIC NUMBER ",I2)') jn    ! long name
+         ctrc2u(jn) = ' '                                       ! units
+      END DO
+
+      REWIND( numnatc )               ! read natrtd
+      READ  ( numnatc, namcfcdia )
+
+      DO jl = 1, jp_cfc_2d
+         jn = jp_cfc0_2d + jl - 1
+         ctrc2d(jn) = cfcdia2d(jl)%snamedia
+         ctrc2l(jn) = cfcdia2d(jl)%lnamedia
+         ctrc2u(jn) = cfcdia2d(jl)%unitdia
+      END DO
+
+
+      IF(lwp) THEN                   ! control print
+         WRITE(numout,*)
+         WRITE(numout,*) ' Namelist : natadd'
+         WRITE(numout,*) '    frequency of outputs for additional arrays nwritedia = ', nwritedia
+         DO jl = 1, jp_cfc_2d
+            jn = jp_cfc0_2d + jl - 1
+            WRITE(numout,*) '   2d output field No : ',jn
+            WRITE(numout,*) '   short name         : ', TRIM(ctrc2d(jn))
+            WRITE(numout,*) '   long name          : ', TRIM(ctrc2l(jn))
+            WRITE(numout,*) '   unit               : ', TRIM(ctrc2u(jn))
+            WRITE(numout,*) ' '
+         END DO
+      ENDIF
+#endif
+
    END SUBROUTINE trc_lsm_cfc