Changeset 1007


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

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

Location:
trunk/NEMO/TOP_SRC/PISCES
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/TOP_SRC/PISCES/p4zflx.F90

    r935 r1007  
    3333     atcco2 = 278.           !: 
    3434 
    35    REAL(wp) ::   & 
    36      xrhoa  = 1.22   ,   &   !: Air density kg/m3 
    37      xcd    = 1.5e-3         !: drag coefficient 
    38  
    3935   REAL(wp) :: & 
    4036     tco2flx = 0.            !: Total flux of carbon per year 
     
    6056      INTEGER, INTENT(in) :: kt 
    6157      INTEGER  ::   ji, jj, jrorr, nspyr 
    62       REAL(wp) ::   zttc, ztx, zty, ztau, zws, zkgwan 
     58      REAL(wp) ::   zttc, zws, zkgwan 
    6359      REAL(wp) ::   zfld, zflu, zfld16, zflu16, zfact 
    6460      REAL(wp) ::   zph, zah2, zbot, zdic, zalk, zschmitto2, zalka, zschmittco2 
     
    131127            zschmitto2  = 1953.4 - 128.0 * zttc + 3.9918 * zttc**2 - 0.050091 * zttc**3 
    132128 
    133             ! Estimation of wind speed as a function of wind stress 
    134             ztx  = utau(ji,jj) * umask(ji,jj,1) 
    135             zty  = vtau(ji,jj) * vmask(ji,jj,1) 
    136             ztau = SQRT( ztx * ztx + zty * zty ) 
    137             zws  = SQRT ( ztau / ( xrhoa * xcd ) ) 
     129            !  wind speed  
     130            zws  = wndm(ji,jj) 
    138131 
    139132      ! Compute the piston velocity for O2 and CO2 
  • trunk/NEMO/TOP_SRC/PISCES/trcini_pisces.F90

    r935 r1007  
    1818   USE par_trc         ! TOP parameters 
    1919   USE sms             ! Source Minus Sink variables 
     20   USE trc 
    2021   USE oce_trc         ! ocean variables 
    2122   USE trp_trc         ! 
    2223   USE p4zche  
    23  
     24   USE lib_mpp 
    2425 
    2526   IMPLICIT NONE 
     
    2829   PUBLIC   trc_ini_pisces   ! called by trcini.F90 module 
    2930 
     31   !! * Module variables 
     32   REAL(wp) :: & 
     33      sco2   =  2.312e-3         , & 
     34      alka0  =  2.423e-3         , & 
     35      oxyg0  =  177.6e-6         , & 
     36      po4    =  2.174e-6         , & 
     37      bioma0 =  1.000e-8         , & 
     38      silic1 =  91.65e-6         , & 
     39      no3    =  31.04e-6 * 7.6 
    3040 
    3141#  include "domzgr_substitute.h90" 
     
    4555      !! ** Purpose :   Initialisation of the PISCES biochemical model 
    4656      !!---------------------------------------------------------------------- 
    47       INTEGER ::   jk 
     57      INTEGER ::   ji, jj, jk, jn 
     58      REAL(wp) ::  caralk, bicarb, co3 
     59 
    4860 
    4961      !!---------------------------------------------------------------------- 
     
    6577 
    6678 
    67       !---------------------------------------------------------------------- 
    68       ! Initialize biological variables  
    69       !---------------------------------------------------------------------- 
     79 
    7080      ! Set biological ratios 
    7181      ! --------------------- 
     
    8090      ndayflxtr = 0      !  Initialize a counter for the computation of chemistry 
    8191 
    82       IF(lwp) WRITE(numout,*) ' Initialisation of PISCES done' 
     92      ! Initialization of tracer concentration in case of  no restart  
     93      !-------------------------------------------------------------- 
     94      IF( .NOT. lrsttr ) THEN   
     95          
     96         trn(:,:,:,jpdic) = sco2 
     97         trn(:,:,:,jpdoc) = bioma0 
     98         trn(:,:,:,jptal) = alka0 
     99         trn(:,:,:,jpoxy) = oxyg0 
     100         trn(:,:,:,jpcal) = bioma0 
     101         trn(:,:,:,jppo4) = po4 / po4r 
     102         trn(:,:,:,jppoc) = bioma0 
     103#  if ! defined key_kriest 
     104         trn(:,:,:,jpgoc) = bioma0 
     105         trn(:,:,:,jpbfe) = bioma0 * 5.e-6 
     106#  else 
     107         trn(:,:,:,jpnum) = bioma0 / ( 6. * xkr_massp ) 
     108#  endif 
     109         trn(:,:,:,jpsil) = silic1 
     110         trn(:,:,:,jpbsi) = bioma0 * 0.15 
     111         trn(:,:,:,jpdsi) = bioma0 * 5.e-6 
     112         trn(:,:,:,jpphy) = bioma0 
     113         trn(:,:,:,jpdia) = bioma0 
     114         trn(:,:,:,jpzoo) = bioma0 
     115         trn(:,:,:,jpmes) = bioma0 
     116         trn(:,:,:,jpfer) = 0.6E-9 
     117         trn(:,:,:,jpsfe) = bioma0 * 5.e-6 
     118         trn(:,:,:,jpdfe) = bioma0 * 5.e-6 
     119         trn(:,:,:,jpnfe) = bioma0 * 5.e-6 
     120         trn(:,:,:,jpnch) = bioma0 * 12. / 55. 
     121         trn(:,:,:,jpdch) = bioma0 * 12. / 55. 
     122         trn(:,:,:,jpno3) = no3 
     123         trn(:,:,:,jpnh4) = bioma0 
     124 
     125         ! Initialization of chemical variables of the carbon cycle 
     126         ! -------------------------------------------------------- 
     127         DO jk = 1, jpk 
     128            DO jj = 1, jpj 
     129               DO ji = 1, jpi 
     130                  caralk = trn(ji,jj,jk,jptal) - borat(ji,jj,jk) / (  1. + 1.E-8 / ( rtrn + akb3(ji,jj,jk) )  ) 
     131                  co3    = ( caralk - trn(ji,jj,jk,jpdic) ) *        tmask(ji,jj,jk)   & 
     132                     &   +                  0.5e-3          * ( 1. - tmask(ji,jj,jk) ) 
     133                  bicarb = ( 2. * trn(ji,jj,jk,jpdic) - caralk ) 
     134                  hi(ji,jj,jk) = ( ak23(ji,jj,jk) * bicarb / co3 ) *        tmask(ji,jj,jk)   & 
     135                     &         +            1.e-9                  * ( 1. - tmask(ji,jj,jk) ) 
     136               END DO 
     137            END DO 
     138         END DO 
     139          
     140      ENDIF 
     141 
     142 
     143 
     144      ! initialize the half saturation constant for silicate 
     145      ! ---------------------------------------------------- 
     146      xksi(:,:)    = 2.e-6 
     147      xksimax(:,:) = xksi(:,:) 
     148 
     149      IF(lwp) WRITE(numout,*) 'Initialization of PISCES tracers done' 
     150      IF(lwp) WRITE(numout,*) ' ' 
     151 
    83152      ! 
    84153   END SUBROUTINE trc_ini_pisces 
  • trunk/NEMO/TOP_SRC/PISCES/trclsm_pisces.F90

    r935 r1007  
    6161      !                               ! Open the namelist file 
    6262      !                               ! ---------------------- 
    63       clname ='namelist.trc.sms' 
     63      clname ='namelist.sms.pisces' 
    6464      CALL ctlopn( numnat, clname, 'OLD', 'FORMATTED', 'SEQUENTIAL',   & 
    6565         &           1, numout, .FALSE., 1 ) 
Note: See TracChangeset for help on using the changeset viewer.