Changeset 2052 for branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/CFC
- Timestamp:
- 2010-08-13T11:15:17+02:00 (14 years ago)
- Location:
- branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/CFC
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/CFC/par_cfc.F90
r2038 r2052 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 … … 34 33 LOGICAL, PUBLIC, PARAMETER :: lk_cfc = .TRUE. !: CFC flag 35 34 INTEGER, PUBLIC, PARAMETER :: jp_cfc = 2 !: number of passive tracers 36 INTEGER, PUBLIC, PARAMETER :: jp_cfc_2d = 2 !: additional 2d output arrays ('key_ diatrc')37 INTEGER, PUBLIC, PARAMETER :: jp_cfc_3d = 0 !: additional 3d output arrays ('key_ diatrc')35 INTEGER, PUBLIC, PARAMETER :: jp_cfc_2d = 2 !: additional 2d output arrays ('key_trc_diaadd') 36 INTEGER, PUBLIC, PARAMETER :: jp_cfc_3d = 0 !: additional 3d output arrays ('key_trc_diaadd') 38 37 INTEGER, PUBLIC, PARAMETER :: jp_cfc_trd = 0 !: number of sms trends for 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 !!====================================================================== -
branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/CFC/trcini_cfc.F90
r2038 r2052 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 … … 47 49 !!---------------------------------------------------------------------- 48 50 49 50 51 IF(lwp) WRITE(numout,*) 51 52 IF(lwp) WRITE(numout,*) ' trc_ini_cfc: initialisation of CFC chemical model' … … 56 57 ! --------------------------------------- 57 58 xphem (:,:) = 0.e0 58 DO jl = 1, jp_cfc 59 jn = jp_cfc0 + jl - 1 60 DO jm = 1, jphem 61 DO js = 1, jpyear 62 p_cfc(js,jm,jn) = 0.0 63 END DO 64 END DO 65 END DO 66 59 p_cfc(:,:,:) = 0.e0 67 60 68 61 ! Initialization of qint in case of no restart … … 74 67 WRITE(numout,*) 'Initialization de qint ; No restart : qint equal zero ' 75 68 ENDIF 69 qint_cfc(:,:,:) = 0.e0 76 70 DO jl = 1, jp_cfc 77 71 jn = jp_cfc0 + jl - 1 78 72 trn (:,:,:,jn) = 0.e0 79 qint_cfc(:,: ,jn) = 0.e080 73 END DO 81 74 ENDIF … … 89 82 IF(lwp) WRITE(numout,*) 'read of formatted file cfc1112atm' 90 83 91 CALL ctl_opn( inum, 'cfc1112.atm', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. )84 CALL ctl_opn( inum, clname, 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. ) 92 85 REWIND(inum) 93 86 … … 147 140 148 141 END SUBROUTINE trc_ini_cfc 149 142 150 143 #else 151 144 !!---------------------------------------------------------------------- -
branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/CFC/trcsms_cfc.F90
r2038 r2052 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
Note: See TracChangeset
for help on using the changeset viewer.