- Timestamp:
- 2010-08-13T09:58:59+02:00 (14 years ago)
- Location:
- trunk/NEMO/TOP_SRC
- Files:
-
- 9 edited
-
C14b/par_c14b.F90 (modified) (3 diffs)
-
CFC/par_cfc.F90 (modified) (3 diffs)
-
CFC/trcctl_cfc.F90 (modified) (2 diffs)
-
CFC/trcini_cfc.F90 (modified) (5 diffs)
-
CFC/trcsms_cfc.F90 (modified) (2 diffs)
-
LOBSTER/par_lobster.F90 (modified) (1 diff)
-
MY_TRC/par_my_trc.F90 (modified) (3 diffs)
-
PISCES/par_pisces.F90 (modified) (4 diffs)
-
par_trc.F90 (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/TOP_SRC/C14b/par_c14b.F90
r1252 r2047 27 27 28 28 IMPLICIT NONE 29 PUBLIC30 29 31 INTEGER, P UBLIC, PARAMETER :: jp_lp= jp_lobster + jp_pisces + jp_cfc !: cum. number of pass. tracers32 INTEGER, P UBLIC, PARAMETER :: jp_lp_2d = jp_lobster_2d + jp_pisces_2d + jp_cfc_2d !:33 INTEGER, P UBLIC, PARAMETER :: jp_lp_3d = jp_lobster_3d + jp_pisces_3d + jp_cfc_3d !:34 INTEGER, P UBLIC, PARAMETER :: jp_lp_trd = jp_lobster_trd + jp_pisces_trd + jp_cfc_trd !:30 INTEGER, PARAMETER :: jp_lb = jp_lobster + jp_pisces + jp_cfc !: cum. number of pass. tracers 31 INTEGER, PARAMETER :: jp_lb_2d = jp_lobster_2d + jp_pisces_2d + jp_cfc_2d !: 32 INTEGER, PARAMETER :: jp_lb_3d = jp_lobster_3d + jp_pisces_3d + jp_cfc_3d !: 33 INTEGER, PARAMETER :: jp_lb_trd = jp_lobster_trd + jp_pisces_trd + jp_cfc_trd !: 35 34 36 35 #if defined key_c14b … … 43 42 INTEGER, PUBLIC, PARAMETER :: jp_c14b_3d = 1 !: additional 3d output arrays ('key_trc_diaadd') 44 43 INTEGER, PUBLIC, PARAMETER :: jp_c14b_trd = 0 !: number of sms trends for C14 45 46 ! assign an index in trc arrays for each CFC prognostic variables 47 INTEGER, PUBLIC, PARAMETER :: jpc14 = jp_lp + 1 !: C14 bomb 44 INTEGER, PUBLIC, PARAMETER :: jpc14 = jp_lb + 1 !: assign an index in trc arrays for C14 bomb 48 45 #else 49 46 !!--------------------------------------------------------------------- … … 58 55 59 56 ! Starting/ending C14 do-loop indices (N.B. no C14 : jp_c14b0 > jp_c14b1 the do-loop are never done) 60 INTEGER, PUBLIC, PARAMETER :: jp_c14b0 = jp_l p+ 1 !: First index of C14 tracer61 INTEGER, PUBLIC, PARAMETER :: jp_c14b1 = jp_l p+ jp_c14b !: Last index of C14 tracer62 INTEGER, PUBLIC, PARAMETER :: jp_c14b0_2d = jp_l p_2d + 1 !: First index of C14 tracer63 INTEGER, PUBLIC, PARAMETER :: jp_c14b1_2d = jp_l p_2d + jp_c14b_2d !: Last index of C14 tracer64 INTEGER, PUBLIC, PARAMETER :: jp_c14b0_3d = jp_l p_3d + 1 !: First index of C14 tracer65 INTEGER, PUBLIC, PARAMETER :: jp_c14b1_3d = jp_l p_3d + jp_c14b_3d !: Last index of C14 tracer66 INTEGER, PUBLIC, PARAMETER :: jp_c14b0_trd = jp_l p_trd + 1 !: First index of C14 tracer67 INTEGER, PUBLIC, PARAMETER :: jp_c14b1_trd = jp_l p_trd + jp_c14b_trd !: Last index of C14 tracer57 INTEGER, PUBLIC, PARAMETER :: jp_c14b0 = jp_lb + 1 !: First index of C14 tracer 58 INTEGER, PUBLIC, PARAMETER :: jp_c14b1 = jp_lb + jp_c14b !: Last index of C14 tracer 59 INTEGER, PUBLIC, PARAMETER :: jp_c14b0_2d = jp_lb_2d + 1 !: First index of C14 tracer 60 INTEGER, PUBLIC, PARAMETER :: jp_c14b1_2d = jp_lb_2d + jp_c14b_2d !: Last index of C14 tracer 61 INTEGER, PUBLIC, PARAMETER :: jp_c14b0_3d = jp_lb_3d + 1 !: First index of C14 tracer 62 INTEGER, PUBLIC, PARAMETER :: jp_c14b1_3d = jp_lb_3d + jp_c14b_3d !: Last index of C14 tracer 63 INTEGER, PUBLIC, PARAMETER :: jp_c14b0_trd = jp_lb_trd + 1 !: First index of C14 tracer 64 INTEGER, PUBLIC, PARAMETER :: jp_c14b1_trd = jp_lb_trd + jp_c14b_trd !: Last index of C14 tracer 68 65 69 66 !!====================================================================== -
trunk/NEMO/TOP_SRC/CFC/par_cfc.F90
r1253 r2047 21 21 22 22 IMPLICIT NONE 23 PUBLIC24 23 25 INTEGER, P UBLIC, PARAMETER :: jp_lp= jp_lobster + jp_pisces !: cumulative number of passive tracers26 INTEGER, P UBLIC, PARAMETER :: jp_lp_2d = jp_lobster_2d + jp_pisces_2d !:27 INTEGER, P UBLIC, PARAMETER :: jp_lp_3d = jp_lobster_3d + jp_pisces_3d !:28 INTEGER, P UBLIC, PARAMETER :: jp_lp_trd = jp_lobster_trd + jp_pisces_trd !:24 INTEGER, PARAMETER :: jp_lc = jp_lobster + jp_pisces !: cumulative number of passive tracers 25 INTEGER, PARAMETER :: jp_lc_2d = jp_lobster_2d + jp_pisces_2d !: 26 INTEGER, PARAMETER :: jp_lc_3d = jp_lobster_3d + jp_pisces_3d !: 27 INTEGER, PARAMETER :: jp_lc_trd = jp_lobster_trd + jp_pisces_trd !: 29 28 30 29 #if defined key_cfc … … 39 38 40 39 ! assign an index in trc arrays for each CFC prognostic variables 41 INTEGER, PUBLIC, PARAMETER :: jpc11 = jp_l p+ 1 !: CFC-1142 INTEGER, PUBLIC, PARAMETER :: jpc12 = jp_l p+ 2 !: CFC-1240 INTEGER, PUBLIC, PARAMETER :: jpc11 = jp_lc + 1 !: CFC-11 41 INTEGER, PUBLIC, PARAMETER :: jpc12 = jp_lc + 2 !: CFC-12 43 42 #else 44 43 !!--------------------------------------------------------------------- … … 53 52 54 53 ! Starting/ending CFC do-loop indices (N.B. no CFC : jp_cfc0 > jp_cfc1 the do-loop are never done) 55 INTEGER, PUBLIC, PARAMETER :: jp_cfc0 = jp_l p+ 1 !: First index of CFC tracers56 INTEGER, PUBLIC, PARAMETER :: jp_cfc1 = jp_l p+ jp_cfc !: Last index of CFC tracers57 INTEGER, PUBLIC, PARAMETER :: jp_cfc0_2d = jp_l p_2d + 1 !: First index of CFC tracers58 INTEGER, PUBLIC, PARAMETER :: jp_cfc1_2d = jp_l p_2d + jp_cfc_2d !: Last index of CFC tracers59 INTEGER, PUBLIC, PARAMETER :: jp_cfc0_3d = jp_l p_3d + 1 !: First index of CFC tracers60 INTEGER, PUBLIC, PARAMETER :: jp_cfc1_3d = jp_l p_3d + jp_cfc_3d !: Last index of CFC tracers61 INTEGER, PUBLIC, PARAMETER :: jp_cfc0_trd = jp_l p_trd + 1 !: First index of CFC tracers62 INTEGER, PUBLIC, PARAMETER :: jp_cfc1_trd = jp_l p_trd + jp_cfc_trd !: Last index of CFC tracers54 INTEGER, PUBLIC, PARAMETER :: jp_cfc0 = jp_lc + 1 !: First index of CFC tracers 55 INTEGER, PUBLIC, PARAMETER :: jp_cfc1 = jp_lc + jp_cfc !: Last index of CFC tracers 56 INTEGER, PUBLIC, PARAMETER :: jp_cfc0_2d = jp_lc_2d + 1 !: First index of CFC tracers 57 INTEGER, PUBLIC, PARAMETER :: jp_cfc1_2d = jp_lc_2d + jp_cfc_2d !: Last index of CFC tracers 58 INTEGER, PUBLIC, PARAMETER :: jp_cfc0_3d = jp_lc_3d + 1 !: First index of CFC tracers 59 INTEGER, PUBLIC, PARAMETER :: jp_cfc1_3d = jp_lc_3d + jp_cfc_3d !: Last index of CFC tracers 60 INTEGER, PUBLIC, PARAMETER :: jp_cfc0_trd = jp_lc_trd + 1 !: First index of CFC tracers 61 INTEGER, PUBLIC, PARAMETER :: jp_cfc1_trd = jp_lc_trd + jp_cfc_trd !: Last index of CFC tracers 63 62 64 63 !!====================================================================== -
trunk/NEMO/TOP_SRC/CFC/trcctl_cfc.F90
r1801 r2047 53 53 ! Check tracer names 54 54 ! ------------------ 55 IF ( ctrcnm(jpc11) /= 'CFC11' .OR. ctrcnm(jpc12) /= 'CFC12' ) THEN 56 ctrcnm(jpc11) = 'CFC11' 57 ctrcnl(jpc11) = 'Chlorofuorocarbone 11 concentration' 58 ctrcnm(jpc12) = 'CFC12' 59 ctrcnl(jpc12) = 'Chlorofuorocarbone 12 concentration' 55 ctrcnm(jpc11) = 'CFC11' 56 ctrcnl(jpc11) = 'Chlorofuorocarbone 11 concentration' 57 58 IF ( jp_cfc == 2 ) THEN 59 ctrcnm(jpc12) = 'CFC12' 60 ctrcnl(jpc12) = 'Chlorofuorocarbone 12 concentration' 60 61 ENDIF 61 62 62 63 IF(lwp) THEN 63 WRITE (numout,*) ' ===>>>> : w a r n i n g <<<<===' 64 WRITE (numout,*) ' ===>>>> : w a r n i n g ' 65 WRITE (numout,*) ' ======= ============= ' 64 66 WRITE (numout,*) ' we force tracer names' 65 67 DO jl = 1, jp_cfc … … 68 70 END DO 69 71 WRITE(numout,*) ' ' 70 ENDIF 72 ENDIF 71 73 72 74 -
trunk/NEMO/TOP_SRC/CFC/trcini_cfc.F90
r1581 r2047 22 22 23 23 PUBLIC trc_ini_cfc ! called by trcini.F90 module 24 25 CHARACTER (len=34) :: clname = 'cfc1112.atm' ! ??? 24 26 25 27 INTEGER :: inum ! unit number … … 55 57 ! --------------------------------------- 56 58 xphem (:,:) = 0.e0 57 DO jl = 1, jp_cfc 58 jn = jp_cfc0 + jl - 1 59 DO jm = 1, jphem 60 DO js = 1, jpyear 61 p_cfc(js,jm,jn) = 0.0 62 END DO 63 END DO 64 END DO 65 59 p_cfc(:,:,:) = 0.e0 66 60 67 61 ! Initialization of qint in case of no restart … … 73 67 WRITE(numout,*) 'Initialization de qint ; No restart : qint equal zero ' 74 68 ENDIF 69 qint_cfc(:,:,:) = 0.e0 75 70 DO jl = 1, jp_cfc 76 71 jn = jp_cfc0 + jl - 1 77 72 trn (:,:,:,jn) = 0.e0 78 qint_cfc(:,: ,jn) = 0.e079 73 END DO 80 74 ENDIF … … 88 82 IF(lwp) WRITE(numout,*) 'read of formatted file cfc1112atm' 89 83 90 CALL ctl_opn( inum, 'cfc1112.atm', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. )84 CALL ctl_opn( inum, clname, 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. ) 91 85 REWIND(inum) 92 86 … … 95 89 END DO 96 90 97 DO jn = 31, 98 ! Read file 98 READ(inum,*) zyy, p_cfc(jn,1,1), p_cfc(jn,1,2), p_cfc(jn,2,1), p_cfc(jn,2,2) 99 WRITE(numout,'(f7.2, 4f8.2)' ) & 91 ! file starts in 1931 do jn represent the year in the century.jhh 92 ! Read file till the end 93 jn = 31 94 DO WHILE ( 1 /= 2 ) 95 READ(inum,*,END=100) zyy, p_cfc(jn,1,1), p_cfc(jn,1,2), p_cfc(jn,2,1), p_cfc(jn,2,2) 96 IF ( lwp) THEN 97 WRITE(numout,'(f7.2, 4f8.2)' ) & 100 98 & zyy, p_cfc(jn,1,1), p_cfc(jn,1,2), p_cfc(jn,2,1), p_cfc(jn,2,2) 99 ENDIF 100 jn = jn + 1 101 101 END DO 102 100 npyear = jn - 1 103 IF ( lwp) WRITE(numout,*) ' ', npyear ,' years read' 102 104 103 105 p_cfc(32,1:2,1) = 5.e-4 ! modify the values of the first years -
trunk/NEMO/TOP_SRC/CFC/trcsms_cfc.F90
r1459 r2047 27 27 PUBLIC trc_sms_cfc ! called in ??? 28 28 29 INTEGER , PUBLIC, PARAMETER :: jpyear = 1 00 ! temporal parameter29 INTEGER , PUBLIC, PARAMETER :: jpyear = 150 ! temporal parameter 30 30 INTEGER , PUBLIC, PARAMETER :: jphem = 2 ! parameter for the 2 hemispheres 31 31 INTEGER , PUBLIC :: ndate_beg ! initial calendar date (aammjj) for CFC 32 32 INTEGER , PUBLIC :: nyear_res ! restoring time constant (year) 33 33 INTEGER , PUBLIC :: nyear_beg ! initial year (aa) 34 INTEGER , PUBLIC :: npyear ! Number of years read in CFC1112 file 34 35 35 36 REAL(wp), PUBLIC, DIMENSION(jpyear,jphem, jp_cfc) :: p_cfc ! partial hemispheric pressure for CFC … … 96 97 ! Temporal interpolation 97 98 ! ---------------------- 98 iyear_beg = nyear + ( nyear_res - 1900 - nyear_beg )99 iyear_beg = nyear - 1900 99 100 IF ( nmonth <= 6 ) THEN 100 iyear_beg = iyear_beg - 2 + nyear_beg101 iyear_beg = iyear_beg - 1 101 102 im1 = 6 - nmonth + 1 102 103 im2 = 6 + nmonth - 1 103 104 ELSE 104 iyear_beg = iyear_beg - 1 + nyear_beg105 105 im1 = 12 - nmonth + 7 106 106 im2 = nmonth - 7 -
trunk/NEMO/TOP_SRC/LOBSTER/par_lobster.F90
r1176 r2047 12 12 13 13 IMPLICIT NONE 14 PUBLIC15 14 16 15 #if defined key_lobster -
trunk/NEMO/TOP_SRC/MY_TRC/par_my_trc.F90
r1253 r2047 31 31 32 32 IMPLICIT NONE 33 PUBLIC34 33 35 INTEGER, P UBLIC, PARAMETER :: jp_lpc= jp_lobster + jp_pisces + jp_cfc + jp_c14b !:36 INTEGER, P UBLIC, PARAMETER :: jp_lpc_2d = jp_lobster_2d + jp_pisces_2d + jp_cfc_2d + jp_c14b_2d !:37 INTEGER, P UBLIC, PARAMETER :: jp_lpc_3d = jp_lobster_3d + jp_pisces_3d + jp_cfc_3d + jp_c14b_3d !:38 INTEGER, P UBLIC, PARAMETER :: jp_lpc_trd = jp_lobster_trd + jp_pisces_trd + jp_cfc_trd + jp_c14b_trd !:34 INTEGER, PARAMETER :: jp_lm = jp_lobster + jp_pisces + jp_cfc + jp_c14b !: 35 INTEGER, PARAMETER :: jp_lm_2d = jp_lobster_2d + jp_pisces_2d + jp_cfc_2d + jp_c14b_2d !: 36 INTEGER, PARAMETER :: jp_lm_3d = jp_lobster_3d + jp_pisces_3d + jp_cfc_3d + jp_c14b_3d !: 37 INTEGER, PARAMETER :: jp_lm_trd = jp_lobster_trd + jp_pisces_trd + jp_cfc_trd + jp_c14b_trd !: 39 38 40 39 #if defined key_my_trc … … 49 48 50 49 ! assign an index in trc arrays for each PTS prognostic variables 51 INTEGER, PUBLIC, PARAMETER :: jpmyt1 = jp_l pc+ 1 !: 1st MY_TRC tracer52 INTEGER, PUBLIC, PARAMETER :: jpmyt2 = jp_l pc+ 2 !: 2nd MY_TRC tracer50 INTEGER, PUBLIC, PARAMETER :: jpmyt1 = jp_lm + 1 !: 1st MY_TRC tracer 51 INTEGER, PUBLIC, PARAMETER :: jpmyt2 = jp_lm + 2 !: 2nd MY_TRC tracer 53 52 54 53 #else … … 64 63 65 64 ! Starting/ending PISCES do-loop indices (N.B. no PISCES : jpl_pcs < jpf_pcs the do-loop are never done) 66 INTEGER, PUBLIC, PARAMETER :: jp_myt0 = jp_l pc+ 1 !: First index of MY_TRC passive tracers67 INTEGER, PUBLIC, PARAMETER :: jp_myt1 = jp_l pc+ jp_my_trc !: Last index of MY_TRC passive tracers68 INTEGER, PUBLIC, PARAMETER :: jp_myt0_2d = jp_l pc_2d + 1 !: First index of MY_TRC passive tracers69 INTEGER, PUBLIC, PARAMETER :: jp_myt1_2d = jp_l pc_2d + jp_my_trc_2d !: Last index of MY_TRC passive tracers70 INTEGER, PUBLIC, PARAMETER :: jp_myt0_3d = jp_l pc_3d + 1 !: First index of MY_TRC passive tracers71 INTEGER, PUBLIC, PARAMETER :: jp_myt1_3d = jp_l pc_3d + jp_my_trc_3d !: Last index of MY_TRC passive tracers72 INTEGER, PUBLIC, PARAMETER :: jp_myt0_trd = jp_l pc_trd + 1 !: First index of MY_TRC passive tracers73 INTEGER, PUBLIC, PARAMETER :: jp_myt1_trd = jp_l pc_trd + jp_my_trc_trd !: Last index of MY_TRC passive tracers65 INTEGER, PUBLIC, PARAMETER :: jp_myt0 = jp_lm + 1 !: First index of MY_TRC passive tracers 66 INTEGER, PUBLIC, PARAMETER :: jp_myt1 = jp_lm + jp_my_trc !: Last index of MY_TRC passive tracers 67 INTEGER, PUBLIC, PARAMETER :: jp_myt0_2d = jp_lm_2d + 1 !: First index of MY_TRC passive tracers 68 INTEGER, PUBLIC, PARAMETER :: jp_myt1_2d = jp_lm_2d + jp_my_trc_2d !: Last index of MY_TRC passive tracers 69 INTEGER, PUBLIC, PARAMETER :: jp_myt0_3d = jp_lm_3d + 1 !: First index of MY_TRC passive tracers 70 INTEGER, PUBLIC, PARAMETER :: jp_myt1_3d = jp_lm_3d + jp_my_trc_3d !: Last index of MY_TRC passive tracers 71 INTEGER, PUBLIC, PARAMETER :: jp_myt0_trd = jp_lm_trd + 1 !: First index of MY_TRC passive tracers 72 INTEGER, PUBLIC, PARAMETER :: jp_myt1_trd = jp_lm_trd + jp_my_trc_trd !: Last index of MY_TRC passive tracers 74 73 75 74 !!====================================================================== -
trunk/NEMO/TOP_SRC/PISCES/par_pisces.F90
r1152 r2047 10 10 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 11 11 !!---------------------------------------------------------------------- 12 USE par_lobster, ONLY : jp_l obster !: number of tracers in LOBSTER13 USE par_lobster, ONLY : jp_l obster_2d !: number of 2D diag in LOBSTER14 USE par_lobster, ONLY : jp_l obster_3d !: number of 3D diag in LOBSTER15 USE par_lobster, ONLY : jp_l obster_trd !: number of biological diag in LOBSTER12 USE par_lobster, ONLY : jp_lpobster !: number of tracers in LOBSTER 13 USE par_lobster, ONLY : jp_lpobster_2d !: number of 2D diag in LOBSTER 14 USE par_lobster, ONLY : jp_lpobster_3d !: number of 3D diag in LOBSTER 15 USE par_lobster, ONLY : jp_lpobster_trd !: number of biological diag in LOBSTER 16 16 17 17 IMPLICIT NONE 18 PUBLIC19 18 20 INTEGER, P UBLIC, PARAMETER :: jp_l = jp_lobster !: cumulative number of already defined TRC21 INTEGER, P UBLIC, PARAMETER :: jp_l_2d = jp_lobster_2d !:22 INTEGER, P UBLIC, PARAMETER :: jp_l_3d = jp_lobster_3d !:23 INTEGER, P UBLIC, PARAMETER :: jp_l_trd = jp_lobster_trd !:19 INTEGER, PARAMETER :: jp_lp = jp_lpobster !: cumulative number of already defined TRC 20 INTEGER, PARAMETER :: jp_lp_2d = jp_lpobster_2d !: 21 INTEGER, PARAMETER :: jp_lp_3d = jp_lpobster_3d !: 22 INTEGER, PARAMETER :: jp_lp_trd = jp_lpobster_trd !: 24 23 25 24 #if defined key_pisces && defined key_kriest … … 37 36 ! WARNING: be carefull about the order when reading the restart 38 37 ! !!gm this warning should be obsolet with IOM 39 INTEGER, PUBLIC, PARAMETER :: jpdic = jp_l + 1 !: dissolved inoganic carbon concentration40 INTEGER, PUBLIC, PARAMETER :: jptal = jp_l + 2 !: total alkalinity41 INTEGER, PUBLIC, PARAMETER :: jpoxy = jp_l + 3 !: oxygen carbon concentration42 INTEGER, PUBLIC, PARAMETER :: jpcal = jp_l + 4 !: calcite concentration43 INTEGER, PUBLIC, PARAMETER :: jppo4 = jp_l + 5 !: phosphate concentration44 INTEGER, PUBLIC, PARAMETER :: jppoc = jp_l + 6 !: small particulate organic phosphate concentration45 INTEGER, PUBLIC, PARAMETER :: jpsil = jp_l + 7 !: silicate concentration46 INTEGER, PUBLIC, PARAMETER :: jpphy = jp_l + 8 !: phytoplancton concentration47 INTEGER, PUBLIC, PARAMETER :: jpzoo = jp_l + 9 !: zooplancton concentration48 INTEGER, PUBLIC, PARAMETER :: jpdoc = jp_l + 10 !: dissolved organic carbon concentration49 INTEGER, PUBLIC, PARAMETER :: jpdia = jp_l + 11 !: Diatoms Concentration50 INTEGER, PUBLIC, PARAMETER :: jpmes = jp_l + 12 !: Mesozooplankton Concentration51 INTEGER, PUBLIC, PARAMETER :: jpbsi = jp_l + 13 !: (big) Silicate Concentration52 INTEGER, PUBLIC, PARAMETER :: jpfer = jp_l + 14 !: Iron Concentration53 INTEGER, PUBLIC, PARAMETER :: jpnum = jp_l + 15 !: Big iron particles Concentration54 INTEGER, PUBLIC, PARAMETER :: jpsfe = jp_l + 16 !: number of particulate organic phosphate concentration55 INTEGER, PUBLIC, PARAMETER :: jpdfe = jp_l + 17 !: Diatoms iron Concentration56 INTEGER, PUBLIC, PARAMETER :: jpdsi = jp_l + 18 !: Diatoms Silicate Concentration57 INTEGER, PUBLIC, PARAMETER :: jpnfe = jp_l + 19 !: Nano iron Concentration58 INTEGER, PUBLIC, PARAMETER :: jpnch = jp_l + 20 !: Nano Chlorophyll Concentration59 INTEGER, PUBLIC, PARAMETER :: jpdch = jp_l + 21 !: Diatoms Chlorophyll Concentration60 INTEGER, PUBLIC, PARAMETER :: jpno3 = jp_l + 22 !: Nitrates Concentration61 INTEGER, PUBLIC, PARAMETER :: jpnh4 = jp_l + 23 !: Ammonium Concentration38 INTEGER, PUBLIC, PARAMETER :: jpdic = jp_lp + 1 !: dissolved inoganic carbon concentration 39 INTEGER, PUBLIC, PARAMETER :: jptal = jp_lp + 2 !: total alkalinity 40 INTEGER, PUBLIC, PARAMETER :: jpoxy = jp_lp + 3 !: oxygen carbon concentration 41 INTEGER, PUBLIC, PARAMETER :: jpcal = jp_lp + 4 !: calcite concentration 42 INTEGER, PUBLIC, PARAMETER :: jppo4 = jp_lp + 5 !: phosphate concentration 43 INTEGER, PUBLIC, PARAMETER :: jppoc = jp_lp + 6 !: small particulate organic phosphate concentration 44 INTEGER, PUBLIC, PARAMETER :: jpsil = jp_lp + 7 !: silicate concentration 45 INTEGER, PUBLIC, PARAMETER :: jpphy = jp_lp + 8 !: phytoplancton concentration 46 INTEGER, PUBLIC, PARAMETER :: jpzoo = jp_lp + 9 !: zooplancton concentration 47 INTEGER, PUBLIC, PARAMETER :: jpdoc = jp_lp + 10 !: dissolved organic carbon concentration 48 INTEGER, PUBLIC, PARAMETER :: jpdia = jp_lp + 11 !: Diatoms Concentration 49 INTEGER, PUBLIC, PARAMETER :: jpmes = jp_lp + 12 !: Mesozooplankton Concentration 50 INTEGER, PUBLIC, PARAMETER :: jpbsi = jp_lp + 13 !: (big) Silicate Concentration 51 INTEGER, PUBLIC, PARAMETER :: jpfer = jp_lp + 14 !: Iron Concentration 52 INTEGER, PUBLIC, PARAMETER :: jpnum = jp_lp + 15 !: Big iron particles Concentration 53 INTEGER, PUBLIC, PARAMETER :: jpsfe = jp_lp + 16 !: number of particulate organic phosphate concentration 54 INTEGER, PUBLIC, PARAMETER :: jpdfe = jp_lp + 17 !: Diatoms iron Concentration 55 INTEGER, PUBLIC, PARAMETER :: jpdsi = jp_lp + 18 !: Diatoms Silicate Concentration 56 INTEGER, PUBLIC, PARAMETER :: jpnfe = jp_lp + 19 !: Nano iron Concentration 57 INTEGER, PUBLIC, PARAMETER :: jpnch = jp_lp + 20 !: Nano Chlorophyll Concentration 58 INTEGER, PUBLIC, PARAMETER :: jpdch = jp_lp + 21 !: Diatoms Chlorophyll Concentration 59 INTEGER, PUBLIC, PARAMETER :: jpno3 = jp_lp + 22 !: Nitrates Concentration 60 INTEGER, PUBLIC, PARAMETER :: jpnh4 = jp_lp + 23 !: Ammonium Concentration 62 61 63 62 #elif defined key_pisces … … 75 74 ! WARNING: be carefull about the order when reading the restart 76 75 ! !!gm this warning should be obsolet with IOM 77 INTEGER, PUBLIC, PARAMETER :: jpdic = jp_l + 1 !: dissolved inoganic carbon concentration78 INTEGER, PUBLIC, PARAMETER :: jptal = jp_l + 2 !: total alkalinity79 INTEGER, PUBLIC, PARAMETER :: jpoxy = jp_l + 3 !: oxygen carbon concentration80 INTEGER, PUBLIC, PARAMETER :: jpcal = jp_l + 4 !: calcite concentration81 INTEGER, PUBLIC, PARAMETER :: jppo4 = jp_l + 5 !: phosphate concentration82 INTEGER, PUBLIC, PARAMETER :: jppoc = jp_l + 6 !: small particulate organic phosphate concentration83 INTEGER, PUBLIC, PARAMETER :: jpsil = jp_l + 7 !: silicate concentration84 INTEGER, PUBLIC, PARAMETER :: jpphy = jp_l + 8 !: phytoplancton concentration85 INTEGER, PUBLIC, PARAMETER :: jpzoo = jp_l + 9 !: zooplancton concentration86 INTEGER, PUBLIC, PARAMETER :: jpdoc = jp_l + 10 !: dissolved organic carbon concentration87 INTEGER, PUBLIC, PARAMETER :: jpdia = jp_l + 11 !: Diatoms Concentration88 INTEGER, PUBLIC, PARAMETER :: jpmes = jp_l + 12 !: Mesozooplankton Concentration89 INTEGER, PUBLIC, PARAMETER :: jpbsi = jp_l + 13 !: (big) Silicate Concentration90 INTEGER, PUBLIC, PARAMETER :: jpfer = jp_l + 14 !: Iron Concentration91 INTEGER, PUBLIC, PARAMETER :: jpbfe = jp_l + 15 !: Big iron particles Concentration92 INTEGER, PUBLIC, PARAMETER :: jpgoc = jp_l + 16 !: big particulate organic phosphate concentration93 INTEGER, PUBLIC, PARAMETER :: jpsfe = jp_l + 17 !: Small iron particles Concentration94 INTEGER, PUBLIC, PARAMETER :: jpdfe = jp_l + 18 !: Diatoms iron Concentration95 INTEGER, PUBLIC, PARAMETER :: jpdsi = jp_l + 19 !: Diatoms Silicate Concentration96 INTEGER, PUBLIC, PARAMETER :: jpnfe = jp_l + 20 !: Nano iron Concentration97 INTEGER, PUBLIC, PARAMETER :: jpnch = jp_l + 21 !: Nano Chlorophyll Concentration98 INTEGER, PUBLIC, PARAMETER :: jpdch = jp_l + 22 !: Diatoms Chlorophyll Concentration99 INTEGER, PUBLIC, PARAMETER :: jpno3 = jp_l + 23 !: Nitrates Concentration100 INTEGER, PUBLIC, PARAMETER :: jpnh4 = jp_l + 24 !: Ammonium Concentration76 INTEGER, PUBLIC, PARAMETER :: jpdic = jp_lp + 1 !: dissolved inoganic carbon concentration 77 INTEGER, PUBLIC, PARAMETER :: jptal = jp_lp + 2 !: total alkalinity 78 INTEGER, PUBLIC, PARAMETER :: jpoxy = jp_lp + 3 !: oxygen carbon concentration 79 INTEGER, PUBLIC, PARAMETER :: jpcal = jp_lp + 4 !: calcite concentration 80 INTEGER, PUBLIC, PARAMETER :: jppo4 = jp_lp + 5 !: phosphate concentration 81 INTEGER, PUBLIC, PARAMETER :: jppoc = jp_lp + 6 !: small particulate organic phosphate concentration 82 INTEGER, PUBLIC, PARAMETER :: jpsil = jp_lp + 7 !: silicate concentration 83 INTEGER, PUBLIC, PARAMETER :: jpphy = jp_lp + 8 !: phytoplancton concentration 84 INTEGER, PUBLIC, PARAMETER :: jpzoo = jp_lp + 9 !: zooplancton concentration 85 INTEGER, PUBLIC, PARAMETER :: jpdoc = jp_lp + 10 !: dissolved organic carbon concentration 86 INTEGER, PUBLIC, PARAMETER :: jpdia = jp_lp + 11 !: Diatoms Concentration 87 INTEGER, PUBLIC, PARAMETER :: jpmes = jp_lp + 12 !: Mesozooplankton Concentration 88 INTEGER, PUBLIC, PARAMETER :: jpbsi = jp_lp + 13 !: (big) Silicate Concentration 89 INTEGER, PUBLIC, PARAMETER :: jpfer = jp_lp + 14 !: Iron Concentration 90 INTEGER, PUBLIC, PARAMETER :: jpbfe = jp_lp + 15 !: Big iron particles Concentration 91 INTEGER, PUBLIC, PARAMETER :: jpgoc = jp_lp + 16 !: big particulate organic phosphate concentration 92 INTEGER, PUBLIC, PARAMETER :: jpsfe = jp_lp + 17 !: Small iron particles Concentration 93 INTEGER, PUBLIC, PARAMETER :: jpdfe = jp_lp + 18 !: Diatoms iron Concentration 94 INTEGER, PUBLIC, PARAMETER :: jpdsi = jp_lp + 19 !: Diatoms Silicate Concentration 95 INTEGER, PUBLIC, PARAMETER :: jpnfe = jp_lp + 20 !: Nano iron Concentration 96 INTEGER, PUBLIC, PARAMETER :: jpnch = jp_lp + 21 !: Nano Chlorophyll Concentration 97 INTEGER, PUBLIC, PARAMETER :: jpdch = jp_lp + 22 !: Diatoms Chlorophyll Concentration 98 INTEGER, PUBLIC, PARAMETER :: jpno3 = jp_lp + 23 !: Nitrates Concentration 99 INTEGER, PUBLIC, PARAMETER :: jpnh4 = jp_lp + 24 !: Ammonium Concentration 101 100 102 101 #else … … 113 112 114 113 ! Starting/ending PISCES do-loop indices (N.B. no PISCES : jpl_pcs < jpf_pcs the do-loop are never done) 115 INTEGER, PUBLIC, PARAMETER :: jp_pcs0 = jp_l + 1 !: First index of PISCES tracers116 INTEGER, PUBLIC, PARAMETER :: jp_pcs1 = jp_l + jp_pisces !: Last index of PISCES tracers117 INTEGER, PUBLIC, PARAMETER :: jp_pcs0_2d = jp_l _2d + 1 !: First index of 2D diag118 INTEGER, PUBLIC, PARAMETER :: jp_pcs1_2d = jp_l _2d + jp_pisces_2d !: Last index of 2D diag119 INTEGER, PUBLIC, PARAMETER :: jp_pcs0_3d = jp_l _3d + 1 !: First index of 3D diag120 INTEGER, PUBLIC, PARAMETER :: jp_pcs1_3d = jp_l _3d + jp_pisces_3d !: Last index of 3d diag121 INTEGER, PUBLIC, PARAMETER :: jp_pcs0_trd = jp_l _trd + 1 !: First index of bio diag122 INTEGER, PUBLIC, PARAMETER :: jp_pcs1_trd = jp_l _trd + jp_pisces_trd !: Last index of bio diag114 INTEGER, PUBLIC, PARAMETER :: jp_pcs0 = jp_lp + 1 !: First index of PISCES tracers 115 INTEGER, PUBLIC, PARAMETER :: jp_pcs1 = jp_lp + jp_pisces !: Last index of PISCES tracers 116 INTEGER, PUBLIC, PARAMETER :: jp_pcs0_2d = jp_lp_2d + 1 !: First index of 2D diag 117 INTEGER, PUBLIC, PARAMETER :: jp_pcs1_2d = jp_lp_2d + jp_pisces_2d !: Last index of 2D diag 118 INTEGER, PUBLIC, PARAMETER :: jp_pcs0_3d = jp_lp_3d + 1 !: First index of 3D diag 119 INTEGER, PUBLIC, PARAMETER :: jp_pcs1_3d = jp_lp_3d + jp_pisces_3d !: Last index of 3d diag 120 INTEGER, PUBLIC, PARAMETER :: jp_pcs0_trd = jp_lp_trd + 1 !: First index of bio diag 121 INTEGER, PUBLIC, PARAMETER :: jp_pcs1_trd = jp_lp_trd + jp_pisces_trd !: Last index of bio diag 123 122 124 123 -
trunk/NEMO/TOP_SRC/par_trc.F90
r1836 r2047 23 23 24 24 IMPLICIT NONE 25 PUBLIC26 25 27 26 ! Passive tracers : Total size 28 27 ! --------------- ! total number of passive tracers, of 2d and 3d output and trend arrays 29 INTEGER, PUBLIC, PARAMETER :: jptra = jp_lobster + jp_pisces + jp_cfc+ jp_my_trc30 INTEGER, PUBLIC, PARAMETER :: jpdia2d = jp_lobster_2d + jp_pisces_2d + jp_cfc_2d+ jp_my_trc_2d31 INTEGER, PUBLIC, PARAMETER :: jpdia3d = jp_lobster_3d + jp_pisces_3d + jp_cfc_3d+ jp_my_trc_3d28 INTEGER, PUBLIC, PARAMETER :: jptra = jp_lobster + jp_pisces + jp_cfc + jp_c14b + jp_my_trc 29 INTEGER, PUBLIC, PARAMETER :: jpdia2d = jp_lobster_2d + jp_pisces_2d + jp_cfc_2d + jp_c14b_2d + jp_my_trc_2d 30 INTEGER, PUBLIC, PARAMETER :: jpdia3d = jp_lobster_3d + jp_pisces_3d + jp_cfc_3d + jp_c14b_3d + jp_my_trc_3d 32 31 ! ! total number of sms diagnostic arrays 33 INTEGER, PUBLIC, PARAMETER :: jpdiabio = jp_lobster_trd + jp_pisces_trd + jp_cfc_trd + jp_my_trc_trd32 INTEGER, PUBLIC, PARAMETER :: jpdiabio = jp_lobster_trd + jp_pisces_trd + jp_cfc_trd + jp_c14b_trd + jp_my_trc_trd 34 33 35 34 ! 1D configuration ("key_c1d") … … 43 42 # if defined key_trc_diatrd 44 43 ! Passive tracers : size for TRP trends diagnotics (used if 'key_trc_diatrd' defined) 45 INTEGER, P ARAMETER :: jptrc_xad = 1 !: x- horizontal advection46 INTEGER, P ARAMETER :: jptrc_yad = 2 !: y- horizontal advection47 INTEGER, P ARAMETER :: jptrc_zad = 3 !: z- vertical advection48 INTEGER, P ARAMETER :: jptrc_xdf = 4 !: lateral diffusion49 INTEGER, P ARAMETER :: jptrc_ydf = 5 !: lateral diffusion50 INTEGER, P ARAMETER :: jptrc_zdf = 6 !: vertical diffusion (Kz)51 INTEGER, P ARAMETER :: jptrc_sbc = 7 !: surface boundary condition44 INTEGER, PUBLIC, PARAMETER :: jptrc_xad = 1 !: x- horizontal advection 45 INTEGER, PUBLIC, PARAMETER :: jptrc_yad = 2 !: y- horizontal advection 46 INTEGER, PUBLIC, PARAMETER :: jptrc_zad = 3 !: z- vertical advection 47 INTEGER, PUBLIC, PARAMETER :: jptrc_xdf = 4 !: lateral diffusion 48 INTEGER, PUBLIC, PARAMETER :: jptrc_ydf = 5 !: lateral diffusion 49 INTEGER, PUBLIC, PARAMETER :: jptrc_zdf = 6 !: vertical diffusion (Kz) 50 INTEGER, PUBLIC, PARAMETER :: jptrc_sbc = 7 !: surface boundary condition 52 51 #if ! defined key_trcldf_eiv && ! defined key_trcdmp 53 INTEGER, P ARAMETER :: jpdiatrc = 7 !: trends: 3*(advection + diffusion ) + sbc52 INTEGER, PUBLIC, PARAMETER :: jpdiatrc = 7 !: trends: 3*(advection + diffusion ) + sbc 54 53 #endif 55 54 #if defined key_trcldf_eiv && defined key_trcdmp 56 INTEGER, P ARAMETER :: jptrc_xei = 8 !: x- horiz. EIV advection57 INTEGER, P ARAMETER :: jptrc_yei = 9 !: y- horiz. EIV advection58 INTEGER, P ARAMETER :: jptrc_zei = 10 !: z- vert. EIV advection59 INTEGER, P ARAMETER :: jptrc_dmp = 11 !: damping60 INTEGER, P ARAMETER :: jpdiatrc = 11 !: trends: 3*(advection + diffusion + eiv ) + sbc + damping55 INTEGER, PUBLIC, PARAMETER :: jptrc_xei = 8 !: x- horiz. EIV advection 56 INTEGER, PUBLIC, PARAMETER :: jptrc_yei = 9 !: y- horiz. EIV advection 57 INTEGER, PUBLIC, PARAMETER :: jptrc_zei = 10 !: z- vert. EIV advection 58 INTEGER, PUBLIC, PARAMETER :: jptrc_dmp = 11 !: damping 59 INTEGER, PUBLIC, PARAMETER :: jpdiatrc = 11 !: trends: 3*(advection + diffusion + eiv ) + sbc + damping 61 60 #endif 62 61 #if defined key_trcldf_eiv && ! defined key_trcdmp 63 INTEGER, P ARAMETER :: jptrc_xei = 8 !: x- horiz. EIV advection64 INTEGER, P ARAMETER :: jptrc_yei = 9 !: y- horiz. EIV advection65 INTEGER, P ARAMETER :: jptrc_zei = 10 !: z- vert. EIV advection66 INTEGER, P ARAMETER :: jpdiatrc = 10 !: trends: 3*(advection + diffusion + eiv ) + sbc62 INTEGER, PUBLIC, PARAMETER :: jptrc_xei = 8 !: x- horiz. EIV advection 63 INTEGER, PUBLIC, PARAMETER :: jptrc_yei = 9 !: y- horiz. EIV advection 64 INTEGER, PUBLIC, PARAMETER :: jptrc_zei = 10 !: z- vert. EIV advection 65 INTEGER, PUBLIC, PARAMETER :: jpdiatrc = 10 !: trends: 3*(advection + diffusion + eiv ) + sbc 67 66 #endif 68 67 #if ! defined key_trcldf_eiv && defined key_trcdmp 69 INTEGER, P ARAMETER :: jptrc_dmp = 8 !: damping70 INTEGER, P ARAMETER :: jpdiatrc = 8 !: trends: 3*(advection + diffusion ) + sbc + damping68 INTEGER, PUBLIC, PARAMETER :: jptrc_dmp = 8 !: damping 69 INTEGER, PUBLIC, PARAMETER :: jpdiatrc = 8 !: trends: 3*(advection + diffusion ) + sbc + damping 71 70 #endif 72 71 #endif
Note: See TracChangeset
for help on using the changeset viewer.
