Changeset 768 for branches/dev_001_GM/NEMO/TOP_SRC/CFC
- Timestamp:
- 2007-12-16T15:46:18+01:00 (17 years ago)
- Location:
- branches/dev_001_GM/NEMO/TOP_SRC/CFC
- Files:
-
- 3 edited
- 1 moved
Legend:
- Unmodified
- Added
- Removed
-
branches/dev_001_GM/NEMO/TOP_SRC/CFC/par_CFC.F90
r765 r768 16 16 PUBLIC 17 17 18 INTEGER, PUBLIC, PARAMETER :: jp_lp = jp_lobster + jp_pisces !: cumulative number of passive tracers 19 18 20 #if defined key_cfc 19 21 !!--------------------------------------------------------------------- … … 27 29 28 30 ! assign an index in trc arrays for each CFC prognostic variables 29 INTEGER, PUBLIC, PARAMETER :: jp_lp = jp_lobster + jp_pisces !: cumulative number of passive tracers30 31 INTEGER, PUBLIC, PARAMETER :: jp11 = jp_lp + 1 !: CFC-11 31 32 INTEGER, PUBLIC, PARAMETER :: jp12 = jp_lp + 2 !: CFC-12 … … 41 42 #endif 42 43 43 ! Starting/ending CFC do-loop indices (N.B. no CFC : jp l_cfc < jpf_cfcthe do-loop are never done)44 INTEGER, PUBLIC, PARAMETER :: jp f_cfc= jp_lp + 1 !: First index of CFC tracers45 INTEGER, PUBLIC, PARAMETER :: jp l_cfc= jp_lp + jp_cfc !: Last index of CFC tracers44 ! Starting/ending CFC do-loop indices (N.B. no CFC : jp_cfc0 > jp_cfc1 the do-loop are never done) 45 INTEGER, PUBLIC, PARAMETER :: jp_cfc0 = jp_lp + 1 !: First index of CFC tracers 46 INTEGER, PUBLIC, PARAMETER :: jp_cfc1 = jp_lp + jp_cfc !: Last index of CFC tracers 46 47 47 48 !!====================================================================== -
branches/dev_001_GM/NEMO/TOP_SRC/CFC/trccfc.F90
r765 r768 25 25 PUBLIC trc_cfc ! called in ??? 26 26 27 INTEGER , P ARAMETER :: jpyear = 100 ! temporal parameter28 INTEGER , P ARAMETER :: jphem = 2 ! parameter for the 2 hemispheres27 INTEGER , PUBLIC, PARAMETER :: jpyear = 100 ! temporal parameter 28 INTEGER , PUBLIC, PARAMETER :: jphem = 2 ! parameter for the 2 hemispheres 29 29 INTEGER , PUBLIC :: ndate_beg ! initial calendar date (aammjj) for CFC 30 30 INTEGER , PUBLIC :: nyear_res ! restoring time constant (year) 31 31 INTEGER , PUBLIC :: nyear_beg ! initial year (aa) 32 32 33 REAL(wp), PUBLIC, DIMENSION(jpyear,jphem, jp f_cfc:jpl_cfc) :: p_cfc ! partial hemispheric pressure for CFC33 REAL(wp), PUBLIC, DIMENSION(jpyear,jphem, jp_cfc0:jp_cfc1) :: p_cfc ! partial hemispheric pressure for CFC 34 34 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: xphem ! spatial interpolation factor for patm 35 REAL(wp), PUBLIC, DIMENSION(jpi,jpj ,jp f_cfc:jpl_cfc) :: qtr ! input function36 REAL(wp), PUBLIC, DIMENSION(jpi,jpj ,jp f_cfc:jpl_cfc) :: qint ! flux function35 REAL(wp), PUBLIC, DIMENSION(jpi,jpj ,jp_cfc0:jp_cfc1) :: qtr ! input function 36 REAL(wp), PUBLIC, DIMENSION(jpi,jpj ,jp_cfc0:jp_cfc1) :: qint ! flux function 37 37 38 38 REAL(wp), DIMENSION(jptra) :: soa1, soa2, soa3, soa4 ! coefficient for solubility of CFC [mol/l/atm] … … 50 50 !!---------------------------------------------------------------------- 51 51 !! NEMO/TOP 2.0 , LOCEAN-IPSL (2007) 52 !! $Id :$52 !! $Id$ 53 53 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 54 54 !!---------------------------------------------------------------------- … … 107 107 108 108 ! !------------! 109 DO jn = jp f_cfc, jpl_cfc! CFC loop !109 DO jn = jp_cfc0, jp_cfc1 ! CFC loop ! 110 110 ! !------------! 111 111 ! time interpolation at time kt … … 180 180 !!--------------------------------------------------------------------- 181 181 182 DO jn = 1, jptra182 DO jn = jp_cfc0, jp_cfc1 183 183 IF ( jn == jp11 ) THEN 184 184 ! coefficient for solubility of CFC11 in mol/l/atm … … 214 214 sca4(jn) = -0.067430 215 215 ENDIF 216 ENDDO 217 218 DO jn = 1, jptra 216 219 217 WRITE(numout,*) 'coefficient for solubility of tracer',ctrcnm(jn) 220 218 WRITE(numout,*) soa1(jn), soa2(jn),soa3(jn), soa4(jn), & -
branches/dev_001_GM/NEMO/TOP_SRC/CFC/trcini_cfc.F90
r766 r768 12 12 !! trc_ini_cfc : CFC model initialisation 13 13 !!---------------------------------------------------------------------- 14 USE oce_trc ! Ocean variables 14 15 USE par_trc ! TOP parameters 16 USE trc ! TOP variables 15 17 USE trccfc ! CFC sms trends 16 18 … … 28 30 !!---------------------------------------------------------------------- 29 31 !! NEMO/TOP 2.0 , LOCEAN-IPSL (2007) 30 !! $Id :$32 !! $Id$ 31 33 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 32 34 !!---------------------------------------------------------------------- … … 52 54 ! Initialization of boundaries conditions 53 55 ! --------------------------------------- 54 pp_cfc(:,:,:) = 0.e055 56 qtr (:,:,:) = 0.e0 56 57 xphem (:,:) = 0.e0 57 DO jn = 1, jptra58 DO jn = jp_cfc0, jp_cfc1 58 59 DO jm = 1, jphem 59 60 DO jl = 1, jpyear … … 71 72 WRITE(numout,*) 'Initialization de qint ; No restart : qint equal zero ' 72 73 ENDIF 73 DO jn = 1, jptra74 DO jn = jp_cfc0, jp_cfc1 74 75 qint(:,:,jn) = 0.e0 75 76 END DO … … 88 89 REWIND(inum) 89 90 90 91 ! Skip over 1st six descriptor lines 92 !----------------------------------- 93 DO jm = 1, 6 91 DO jm = 1, 6 ! Skip over 1st six descriptor lines 94 92 READ(inum,'(1x)') 95 93 END DO 96 94 97 98 ! Read file 99 ! --------- 100 DO jn = 31, 98 95 DO jn = 31, 98 ! Read file 101 96 READ(inum,*) zyy, p_cfc(jn,1,jp11), p_cfc(jn,1,jp12), & 102 97 & p_cfc(jn,2,jp11), p_cfc(jn,2,jp12) … … 106 101 END DO 107 102 108 p_cfc(32,1:2,jp11) = 5.e-4 103 p_cfc(32,1:2,jp11) = 5.e-4 ! modify the values of the first years 109 104 p_cfc(33,1:2,jp11) = 8.e-4 110 105 p_cfc(34,1:2,jp11) = 1.e-6 … … 115 110 p_cfc(39,1:2,jp11) = 1.e-2 116 111 117 118 IF(lwp) THEN 112 IF(lwp) THEN ! Control print 119 113 WRITE(numout,*) 120 114 WRITE(numout,*) ' Year p11HN p11HS p12HN p12HS ' -
branches/dev_001_GM/NEMO/TOP_SRC/CFC/trclsm_cfc.F90
r764 r768 1 MODULE trclsm_cfc 2 !!====================================================================== 3 !! *** MODULE trclsm_cfc *** 4 !! TOP : initialisation of some run parameters for CFC chemical model 5 !!====================================================================== 6 !! History : 2.0 ! 2007-12 (C. Ethe, G. Madec) from trclsm.cfc.h90 1 7 !!---------------------------------------------------------------------- 2 !! *** trclsm.cfc.h90 *** 3 !! TOP : Definition some run parameter for CFC chemical model 8 #if defined key_cfc 4 9 !!---------------------------------------------------------------------- 5 !! History : 1.0 ! 2003-08 (C. Ethe) Original code10 !! 'key_cfc' CFC tracers 6 11 !!---------------------------------------------------------------------- 12 !! trc_lsm_cfc : CFC model initialisation 13 !!---------------------------------------------------------------------- 14 USE oce_trc ! Ocean variables 15 USE par_trc ! TOP parameters 16 USE trc ! TOP variables 17 USE sms ! sms trends 18 USE trccfc ! CFC specific variable 19 20 IMPLICIT NONE 21 PRIVATE 22 23 PUBLIC trc_lsm_cfc ! called by trclsm.F90 module 7 24 8 25 !!---------------------------------------------------------------------- 9 !! NEMO/TOP 1.0 , LOCEAN-IPSL (2005)10 !! $Id 26 !! NEMO/TOP 2.0 , LOCEAN-IPSL (2007) 27 !! $Id$ 11 28 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 12 29 !!---------------------------------------------------------------------- … … 14 31 CONTAINS 15 32 16 SUBROUTINE trc_lsm 33 SUBROUTINE trc_lsm_cfc 17 34 !!------------------------------------------------------------------- 18 !! *** ROUTINE trc_lsm ***35 !! *** ROUTINE trc_lsm_cfc *** 19 36 !! 20 37 !! ** Purpose : Definition some run parameter for CFC model … … 50 67 IF(lwp) WRITE(numout,*) ' initial year (aa) nyear_beg = ', nyear_beg 51 68 ! 52 END SUBROUTINE trc_lsm 69 END SUBROUTINE trc_lsm_cfc 70 71 #else 72 !!---------------------------------------------------------------------- 73 !! Dummy module : No CFC 74 !!---------------------------------------------------------------------- 75 CONTAINS 76 SUBROUTINE trc_lsm_cfc ! Empty routine 77 END SUBROUTINE trc_lsm_cfc 78 #endif 79 80 !!====================================================================== 81 END MODULE trclsm_cfc
Note: See TracChangeset
for help on using the changeset viewer.