Changeset 1283 for trunk/NEMO/TOP_SRC/TRP

Timestamp:

2009-02-03T14:15:21+01:00 (15 years ago)

Author:

cetlod

Message:

Initialization of biological trends variable name, see ticket:312

Location:

trunk/NEMO/TOP_SRC/TRP

Files:

• trdmld_trc.F90 (modified) (3 diffs)
• trdmld_trc_rst.F90 (modified) (2 diffs)

trunk/NEMO/TOP_SRC/TRP/trdmld_trc.F90

@@ -1338 +1338 @@
       nmoymltrdbio   = 0
       tmltrd_sum_bio     (:,:,:) = 0.e0     ;   tmltrd_csum_ln_bio (:,:,:) = 0.e0
+      DO jl = 1, jp_lobster_trd
+          ctrd_bio(jl,1) = ctrbil(jl)   ! long name
+          ctrd_bio(jl,2) = ctrbio(jl)   ! short name
+       ENDDO
 #endif
 
@@ -1458 +1462 @@
 
 #if defined key_lobster
-
-          ctrd_bio(1,:) = "NO3PHY"
-          ctrd_bio(2,:) = "NH4PHY"
-          ctrd_bio(3,:) = "PHYNH4"
-          ctrd_bio(4,:) = "PHYDOM"
-          ctrd_bio(5,:) = "PHYZOO"
-          ctrd_bio(6,:) = "PHYDET"
-          ctrd_bio(7,:) = "DETZOO"
-          ctrd_bio(8,:) = "DETSED"
-          ctrd_bio(9,:) = "ZOODET"
-          ctrd_bio(10,:) = "ZOOBOD"
-          ctrd_bio(11,:) = "ZOONH4"
-          ctrd_bio(12,:) = "ZOODOM"
-          ctrd_bio(13,:) = "NH4NO3"
-          ctrd_bio(14,:) = "DOMNH4"
-          ctrd_bio(15,:) = "DETNH4"
-          ctrd_bio(16,:) = "DETDOM"
-          ctrd_bio(17,:) = "SEDNO3"
-
-
           !-- Create a NetCDF file and enter the define mode
           CALL dia_nam( clhstnam, ntrd_trc, 'trdbio' )
           CALL histbeg( clhstnam, jpi, glamt, jpj, gphit,                                            &
              &             1, jpi, 1, jpj, 0, zjulian, rdt, nh_tb, nidtrdbio, domain_id=nidom )
-
 #endif
 
@@ -1528 +1511 @@
 
 #if defined key_lobster
-      DO jl = 1, jpdiabio
+      DO jl = 1, jp_lobster_trd
          CALL histdef(nidtrdbio, TRIM("ML_"//ctrd_bio(jl,2)), TRIM(clmxl//" ML_"//ctrd_bio(jl,1))   ,            &
              &    cltrcu, jpi, jpj, nh_tb, 1  , 1, 1  , -99 , 32, clop, zsto, zout ) ! IOIPSL: time mean

trunk/NEMO/TOP_SRC/TRP/trdmld_trc_rst.F90

@@ -51 +51 @@
            WRITE(clkt,'(i8.8)') nitrst
          ENDIF
-         clname = TRIM(cexper)//"_"//TRIM(ADJUSTL(clkt))//"_restart_mld_trc"
+         clname = TRIM(cexper)//"_"//TRIM(ADJUSTL(clkt))//"_"//TRIM(cn_trdrst_trc_out)
          IF(lwp) WRITE(numout,*) '             open ocean restart_mld_trc NetCDF  '//clname
          CALL iom_open( clname, nummldw_trc, ldwrt = .TRUE., kiolib = jprstlib )
@@ -140 +140 @@
       ENDIF
       
-      CALL iom_open( 'restart_mld_trc', inum, kiolib = jprstlib ) 
+      CALL iom_open( cn_trdrst_trc_in, inum, kiolib = jprstlib ) 
       
       IF( ln_trdmld_trc_instant ) THEN