Changeset 1004


Ignore:
Timestamp:
2008-05-30T11:47:14+02:00 (13 years ago)
Author:
cetlod
Message:

Update CFC modules to take into account the re-organization of TOP initialization phase, see ticket 168

Location:
trunk/NEMO/TOP_SRC/CFC
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/TOP_SRC/CFC/par_cfc.F90

    r933 r1004  
    2424   LOGICAL, PUBLIC, PARAMETER ::   lk_cfc     = .TRUE.      !: CFC flag  
    2525   INTEGER, PUBLIC, PARAMETER ::   jp_cfc     =  2          !: number of passive tracers 
    26    INTEGER, PUBLIC, PARAMETER ::   jp_cfc_2d  =  1          !: additional 2d output arrays ('key_trc_diaadd') 
    27    INTEGER, PUBLIC, PARAMETER ::   jp_cfc_3d  =  1          !: additional 3d output arrays ('key_trc_diaadd') 
     26   INTEGER, PUBLIC, PARAMETER ::   jp_cfc_2d  =  0          !: additional 2d output arrays ('key_trc_diaadd') 
     27   INTEGER, PUBLIC, PARAMETER ::   jp_cfc_3d  =  0          !: additional 3d output arrays ('key_trc_diaadd') 
    2828   INTEGER, PUBLIC, PARAMETER ::   jp_cfc_trd =  0          !: number of sms trends for CFC 
    2929    
  • trunk/NEMO/TOP_SRC/CFC/trcini_cfc.F90

    r933 r1004  
    5252      IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~' 
    5353 
     54 
    5455      ! Initialization of boundaries conditions 
    5556      ! ---------------------------------------  
     
    7374         ENDIF 
    7475         DO jn = jp_cfc0, jp_cfc1 
    75             qint(:,:,jn) = 0.e0 
     76            trn(:,:,:,jn) = 0.e0 
     77            qint(:,: ,jn) = 0.e0 
    7678         END DO 
    7779      ENDIF 
     
    134136      END DO 
    135137      ! 
     138 
     139      IF(lwp) WRITE(numout,*) 'Initialization of CFC tracers done' 
     140      IF(lwp) WRITE(numout,*) ' ' 
     141 
    136142   END SUBROUTINE trc_ini_cfc 
    137143    
  • trunk/NEMO/TOP_SRC/CFC/trclsm_cfc.F90

    r933 r1004  
    4242      !! ** input   :   Namelist namcfc 
    4343      !!---------------------------------------------------------------------- 
    44       CHARACTER (len=32) ::   clname = 'namelist.trc.sms' 
     44      CHARACTER (len=32) ::   clname = 'namelist.sms.cfc' 
    4545      INTEGER ::   numnat 
    4646      !! 
  • trunk/NEMO/TOP_SRC/CFC/trcsms_cfc.F90

    r933 r1004  
    3636   REAL(wp), PUBLIC, DIMENSION(jpi,jpj      ,jp_cfc0:jp_cfc1) ::   qint     ! flux function 
    3737 
    38    REAL(wp), DIMENSION(jp_cfc) ::   soa1, soa2, soa3, soa4   ! coefficient for solubility of CFC [mol/l/atm] 
    39    REAL(wp), DIMENSION(jp_cfc) ::   sob1, sob2, sob3         !    "               " 
    40    REAL(wp), DIMENSION(jp_cfc) ::   sca1, sca2, sca3, sca4   ! coefficients for schmidt number in degre Celcius 
     38   REAL(wp), DIMENSION(jp_cfc0:jp_cfc1) ::   soa1, soa2, soa3, soa4   ! coefficient for solubility of CFC [mol/l/atm] 
     39   REAL(wp), DIMENSION(jp_cfc0:jp_cfc1) ::   sob1, sob2, sob3         !    "               " 
     40   REAL(wp), DIMENSION(jp_cfc0:jp_cfc1) ::   sca1, sca2, sca3, sca4   ! coefficients for schmidt number in degre Celcius 
    4141       
    4242   !                          ! coefficients for conversion 
     
    4545   REAL(wp) ::   xconv3 = 1.0e+3       ! conversion from mol/l/atm to mol/m3/atm 
    4646   REAL(wp) ::   xconv4 = 1.0e-12      ! conversion from mol/m3/atm to mol/m3/pptv  
    47  
    48    REAL(wp) ::   xrhoa  = 1.22         ! Air density kg/m3 
    49    REAL(wp) ::   xcd    = 1.5e-3       ! drag coefficient 
    5047 
    5148   !! * Substitutions 
     
    8885      REAL(wp) ::   zca_cfc   ! concentration at equilibrium 
    8986      REAL(wp) ::   zak_cfc   ! transfert coefficients 
    90  
    91       REAL(wp) ::   ztx, zty, ztau 
    92       REAL(wp), DIMENSION(jpi,jpj)        :: zws     ! wind speed 
    9387 
    9488      REAL(wp), DIMENSION(jphem,jp_cfc)   ::   zpatm       ! atmospheric function 
     
    110104      ENDIF 
    111105      iyear_end = iyear_beg + 1 
    112  
    113  
    114       ! Estimation of wind speed as a function of wind stress 
    115 !CDIR NOVERRCHK 
    116       DO jj = 1, jpj 
    117 !CDIR NOVERRCHK 
    118          DO ji = 1, jpi 
    119             ztx  = utau(ji,jj) * umask(ji,jj,1) 
    120             zty  = vtau(ji,jj) * vmask(ji,jj,1) 
    121             ztau = SQRT( ztx * ztx + zty * zty ) 
    122             zws(ji,jj) = SQRT ( ztau / ( xrhoa * xcd ) ) 
    123          ENDDO 
    124       ENDDO 
    125106 
    126107      !                                                         !------------! 
     
    164145 
    165146               !    speed transfert : formulae of wanninkhof 1992 
    166                zv2     = zws(ji,jj) * zws(ji,jj) 
     147               zv2     = wndm(ji,jj) * wndm(ji,jj) 
    167148               zsch    = zsch / 660. 
    168149               zak_cfc = ( 0.39 * xconv2 * zv2 / SQRT(zsch) ) * tmask(ji,jj,1) 
     
    189170   END SUBROUTINE trc_sms_cfc 
    190171 
    191  
    192172   SUBROUTINE trc_cfc_cst 
    193173      !!--------------------------------------------------------------------- 
     
    196176      !! ** Purpose : sets constants for CFC model 
    197177      !!--------------------------------------------------------------------- 
    198       INTEGER ::   jn 
    199       !!--------------------------------------------------------------------- 
    200  
    201       DO jn = jp_cfc0, jp_cfc1 
    202          IF ( jn == jp11 ) THEN 
    203             ! coefficient for solubility of CFC11 in  mol/l/atm 
    204             soa1(jn) = -229.9261  
    205             soa2(jn) =  319.6552 
    206             soa3(jn) =  119.4471 
    207             soa4(jn) =   -1.39165 
    208             sob1(jn) =   -0.142382 
    209             sob2(jn) =    0.091459 
    210             sob3(jn) =   -0.0157274 
    211              
    212             ! coefficients for schmidt number in degre Celcius 
    213             sca1(jn) = 3501.8 
    214             sca2(jn) = -210.31 
    215             sca3(jn) =    6.1851 
    216             sca4(jn) =   -0.07513 
    217  
    218          ELSE IF( jn == jp12 ) THEN 
    219  
    220             ! coefficient for solubility of CFC12 in  mol/l/atm 
    221             soa1(jn) = -218.0971 
    222             soa2(jn) =  298.9702 
    223             soa3(jn) =  113.8049 
    224             soa4(jn) =   -1.39165 
    225             sob1(jn) =   -0.143566 
    226             sob2(jn) =    0.091015 
    227             sob3(jn) =   -0.0153924 
    228                          
    229             ! coefficients for schmidt number in degre Celcius 
    230             sca1(jn) =  3845.4  
    231             sca2(jn) =  -228.95 
    232             sca3(jn) =     6.1908  
    233             sca4(jn) =    -0.067430 
    234          ENDIF 
    235  
    236          WRITE(numout,*) 'coefficient for solubility of tracer',ctrcnm(jn) 
    237          WRITE(numout,*) soa1(jn), soa2(jn),soa3(jn), soa4(jn), & 
    238             &           sob1(jn), sob2(jn),sob3(jn) 
    239          WRITE(numout,*)  
    240          WRITE(numout,*) 'coefficient for schmidt of tracer',ctrcnm(jn) 
    241          WRITE(numout,*) sca1(jn), sca2(jn),sca3(jn), sca4(jn) 
    242       END DO 
    243       ! 
     178 
     179      ! coefficient for solubility of CFC11/CFC12 in mol/l/atm 
     180 
     181        soa1(jp11) = -229.9261  
     182        soa2(jp11) =  319.6552 
     183        soa3(jp11) =  119.4471 
     184        soa4(jp11) =  -1.39165 
     185        sob1(jp11) =  -0.142382 
     186        sob2(jp11) =   0.091459 
     187        sob3(jp11) =  -0.0157274 
     188 
     189        soa1(jp12) = -218.0971 
     190        soa2(jp12) =  298.9702 
     191        soa3(jp12) =  113.8049 
     192        soa4(jp12) =  -1.39165 
     193        sob1(jp12) =  -0.143566 
     194        sob2(jp12) =   0.091015 
     195        sob3(jp12) =  -0.0153924 
     196 
     197          
     198      ! coefficients for schmidt number in degre Celcius 
     199        sca1(jp11) = 3501.8 
     200        sca2(jp11) = -210.31 
     201        sca3(jp11) =  6.1851 
     202        sca4(jp11) = -0.07513 
     203 
     204        sca1(jp12) =  3845.4  
     205        sca2(jp12) =  -228.95 
     206        sca3(jp12) =  6.1908  
     207        sca4(jp12) =  -0.067430 
     208 
    244209   END SUBROUTINE trc_cfc_cst 
    245210    
Note: See TracChangeset for help on using the changeset viewer.