Changeset 2097
- Timestamp:
- 2010-09-15T14:53:03+02:00 (13 years ago)
- Location:
- branches/DEV_r1879_FCM/NEMOGCM/NEMO
- Files:
-
- 14 edited
-
OPA_SRC/OBC/obcfla.F90 (modified) (10 diffs)
-
OPA_SRC/SBC/cpl_oasis3.F90 (modified) (1 diff)
-
OPA_SRC/SBC/sbccpl.F90 (modified) (5 diffs)
-
OPA_SRC/SBC/sbcice_lim_2.F90 (modified) (1 diff)
-
TOP_SRC/C14b/par_c14b.F90 (modified) (3 diffs)
-
TOP_SRC/CFC/par_cfc.F90 (modified) (3 diffs)
-
TOP_SRC/CFC/trcctl_cfc.F90 (modified) (2 diffs)
-
TOP_SRC/CFC/trcini_cfc.F90 (modified) (5 diffs)
-
TOP_SRC/CFC/trcsms_cfc.F90 (modified) (2 diffs)
-
TOP_SRC/LOBSTER/par_lobster.F90 (modified) (1 diff)
-
TOP_SRC/MY_TRC/par_my_trc.F90 (modified) (3 diffs)
-
TOP_SRC/MY_TRC/trcctl_my_trc.F90 (modified) (1 diff)
-
TOP_SRC/PISCES/par_pisces.F90 (modified) (4 diffs)
-
TOP_SRC/par_trc.F90 (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
branches/DEV_r1879_FCM/NEMOGCM/NEMO/OPA_SRC/OBC/obcfla.F90
r2007 r2097 78 78 !!------------------------------------------------------------------------------ 79 79 !! * Local declaration 80 INTEGER :: ji, jj , jk! dummy loop indices80 INTEGER :: ji, jj ! dummy loop indices 81 81 !!------------------------------------------------------------------------------ 82 82 … … 85 85 ua_e(ji,jj) = ( ubtfoe(jj) * hur(ji,jj) + sqrt( grav*hur(ji,jj) ) & 86 86 & * ( ( sshn_e(ji,jj) + sshn_e(ji+1,jj) ) * 0.5 & 87 & - sshfoe(jj) ) ) * uemsk(jj, jk)87 & - sshfoe(jj) ) ) * uemsk(jj,1) 88 88 END DO 89 89 END DO … … 95 95 ssha_e(ji,jj) = ssha_e(ji,jj) * ( 1. - temsk(jj,1) ) & 96 96 & + temsk(jj,1) * sshfoe(jj) 97 va_e(ji,jj) = vbtfoe(jj) * hvr(ji,jj) * uemsk(jj, jk)97 va_e(ji,jj) = vbtfoe(jj) * hvr(ji,jj) * uemsk(jj,1) 98 98 END DO 99 99 END DO … … 114 114 !!------------------------------------------------------------------------------ 115 115 !! * Local declaration 116 INTEGER :: ji, jj , jk! dummy loop indices116 INTEGER :: ji, jj ! dummy loop indices 117 117 !!------------------------------------------------------------------------------ 118 118 … … 121 121 ua_e(ji,jj) = ( ubtfow(jj) * hur(ji,jj) - sqrt( grav * hur(ji,jj) ) & 122 122 & * ( ( sshn_e(ji,jj) + sshn_e(ji+1,jj) ) * 0.5 & 123 & - sshfow(jj) ) ) * uwmsk(jj, jk)124 va_e(ji,jj) = vbtfow(jj) * hvr(ji,jj) * uwmsk(jj, jk)123 & - sshfow(jj) ) ) * uwmsk(jj,1) 124 va_e(ji,jj) = vbtfow(jj) * hvr(ji,jj) * uwmsk(jj,1) 125 125 END DO 126 126 DO jj = 1, jpj … … 147 147 !!------------------------------------------------------------------------------ 148 148 !! * Local declaration 149 INTEGER :: ji, jj , jk! dummy loop indices149 INTEGER :: ji, jj ! dummy loop indices 150 150 !!------------------------------------------------------------------------------ 151 151 … … 154 154 va_e(ji,jj) = ( vbtfon(ji) * hvr(ji,jj) + sqrt( grav * hvr(ji,jj) ) & 155 155 & * ( ( sshn_e(ji,jj) + sshn_e(ji,jj+1) ) * 0.5 & 156 & - sshfon(ji) ) ) * vnmsk(ji, jk)156 & - sshfon(ji) ) ) * vnmsk(ji,1) 157 157 END DO 158 158 END DO … … 164 164 ssha_e(ji,jj) = ssha_e(ji,jj) * ( 1. - tnmsk(ji,1) ) & 165 165 & + sshfon(ji) * tnmsk(ji,1) 166 ua_e(ji,jj) = ubtfon(ji) * hur(ji,jj) * vnmsk(ji, jk)166 ua_e(ji,jj) = ubtfon(ji) * hur(ji,jj) * vnmsk(ji,1) 167 167 END DO 168 168 END DO … … 182 182 !!------------------------------------------------------------------------------ 183 183 !! * Local declaration 184 INTEGER :: ji, jj , jk! dummy loop indices184 INTEGER :: ji, jj ! dummy loop indices 185 185 186 186 !!------------------------------------------------------------------------------ … … 190 190 va_e(ji,jj) = ( vbtfos(ji) * hvr(ji,jj) - sqrt( grav * hvr(ji,jj) ) & 191 191 & * ( ( sshn_e(ji,jj) + sshn_e(ji,jj+1) ) * 0.5 & 192 & - sshfos(ji) ) ) * vsmsk(ji, jk)193 ua_e(ji,jj) = ubtfos(ji) * hur(ji,jj) * vsmsk(ji, jk)192 & - sshfos(ji) ) ) * vsmsk(ji,1) 193 ua_e(ji,jj) = ubtfos(ji) * hur(ji,jj) * vsmsk(ji,1) 194 194 END DO 195 195 DO ji = 1, jpi -
branches/DEV_r1879_FCM/NEMOGCM/NEMO/OPA_SRC/SBC/cpl_oasis3.F90
r1698 r2097 27 27 USE mod_prism_proto ! OASIS3 prism module 28 28 USE mod_prism_def_partition_proto! OASIS3 prism module for partitioning 29 USE mod_prism_grids_writing ! OASIS3 prism module for writing grid files30 29 USE mod_prism_put_proto ! OASIS3 prism module for snding 31 30 USE mod_prism_get_proto ! OASIS3 prism module for receiving 32 USE mod_prism_grids_writing ! OASIS3 prism module for writing grids33 31 USE par_oce ! ocean parameters 34 32 USE dom_oce ! ocean space and time domain -
branches/DEV_r1879_FCM/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90
r2045 r2097 152 152 INTEGER , DIMENSION( jprcv) :: nrcvinfo ! OASIS info argument 153 153 154 #if ! defined key_lim2 && ! defined key_lim3 155 ! quick patch to be able to run the coupled model without sea-ice... 156 INTEGER, PARAMETER :: jpl = 1 157 REAL(wp), DIMENSION(jpi,jpj ) :: hicif, hsnif, u_ice, v_ice,fr1_i0,fr2_i0 158 REAL(wp), DIMENSION(jpi,jpj,jpl) :: tn_ice, alb_ice 159 REAL(wp) :: lfus 160 #endif 161 154 162 !! Substitution 155 163 # include "vectopt_loop_substitute.h90" … … 253 261 254 262 ! default definitions of srcv 255 srcv(:)%laction = .FALSE. ; srcv(:)%clgrid = 'T' ; srcv(:)%nsgn = 1 263 srcv(:)%laction = .FALSE. ; srcv(:)%clgrid = 'T' ; srcv(:)%nsgn = 1. 256 264 257 265 ! ! ------------------------- ! … … 438 446 439 447 ! default definitions of nsnd 440 ssnd(:)%laction = .FALSE. ; ssnd(:)%clgrid = 'T' ; ssnd(:)%nsgn = 1 448 ssnd(:)%laction = .FALSE. ; ssnd(:)%clgrid = 'T' ; ssnd(:)%nsgn = 1. 441 449 442 450 ! ! ------------------------- ! … … 491 499 ssnd(jps_ocz1)%clname = 'O_OCurz1' ; ssnd(jps_ivz1)%clname = 'O_IVelz1' 492 500 ! 493 ssnd(jps_ocx1:jps_ivz1)%nsgn = -1 501 ssnd(jps_ocx1:jps_ivz1)%nsgn = -1. ! vectors: change of the sign at the north fold 494 502 495 503 IF( cn_snd_crt(4) /= 'T' ) CALL ctl_stop( 'cn_snd_crt(4) must be equal to T' ) … … 713 721 ! ! non solar heat flux over the ocean (qns) 714 722 IF( srcv(jpr_qnsoce)%laction ) qns(:,:) = frcv(:,:,jpr_qnsoce) 715 IF( srcv(jpr_qnsmix)%laction ) qns(:,:) = frcv(:,:,jpr_qnsmix) 716 qns(:,:) = qns(:,:) - frcv(:,:,jpr_snow) * lfus ! add the latent heat of solid precip. melting 723 IF( srcv(jpr_qnsmix)%laction ) qns(:,:) = frcv(:,:,jpr_qnsmix) 724 ! add the latent heat of solid precip. melting 725 IF( srcv(jpr_snow )%laction ) qns(:,:) = qns(:,:) - frcv(:,:,jpr_snow) * lfus 717 726 718 727 ! ! solar flux over the ocean (qsr) -
branches/DEV_r1879_FCM/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim_2.F90
r2066 r2097 180 180 #endif 181 181 #if defined key_coupled 182 IF( ksbc == 5 ) CALL sbc_cpl_ice_flx( frld ,&182 IF( ksbc == 5 ) CALL sbc_cpl_ice_flx( reshape( frld, (/jpi,jpj,1/) ), & 183 183 & qns_tot, qns_ice, qsr_tot , qsr_ice, & 184 184 & emp_tot, emp_ice, dqns_ice, sprecip, & -
branches/DEV_r1879_FCM/NEMOGCM/NEMO/TOP_SRC/C14b/par_c14b.F90
r1252 r2097 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 !!====================================================================== -
branches/DEV_r1879_FCM/NEMOGCM/NEMO/TOP_SRC/CFC/par_cfc.F90
r1253 r2097 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 !!====================================================================== -
branches/DEV_r1879_FCM/NEMOGCM/NEMO/TOP_SRC/CFC/trcctl_cfc.F90
r2007 r2097 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 -
branches/DEV_r1879_FCM/NEMOGCM/NEMO/TOP_SRC/CFC/trcini_cfc.F90
r1581 r2097 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 -
branches/DEV_r1879_FCM/NEMOGCM/NEMO/TOP_SRC/CFC/trcsms_cfc.F90
r1459 r2097 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 -
branches/DEV_r1879_FCM/NEMOGCM/NEMO/TOP_SRC/LOBSTER/par_lobster.F90
r1176 r2097 12 12 13 13 IMPLICIT NONE 14 PUBLIC15 14 16 15 #if defined key_lobster -
branches/DEV_r1879_FCM/NEMOGCM/NEMO/TOP_SRC/MY_TRC/par_my_trc.F90
r1253 r2097 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 !!====================================================================== -
branches/DEV_r1879_FCM/NEMOGCM/NEMO/TOP_SRC/MY_TRC/trcctl_my_trc.F90
r1264 r2097 40 40 IF(lwp) WRITE(numout,*) ' use COLOR tracer ' 41 41 42 DO j n= 1, jp_my_trc42 DO jl = 1, jp_my_trc 43 43 jn = jp_myt0 + jl - 1 44 44 WRITE(ctrcnm(jn),'(a,i2.2)') 'CLR',jn -
branches/DEV_r1879_FCM/NEMOGCM/NEMO/TOP_SRC/PISCES/par_pisces.F90
r1152 r2097 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_lobster !: cumulative number of already defined TRC 20 INTEGER, PARAMETER :: jp_lp_2d = jp_lobster_2d !: 21 INTEGER, PARAMETER :: jp_lp_3d = jp_lobster_3d !: 22 INTEGER, PARAMETER :: jp_lp_trd = jp_lobster_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 -
branches/DEV_r1879_FCM/NEMOGCM/NEMO/TOP_SRC/par_trc.F90
r2007 r2097 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.
