New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 6509 for branches/NERC – NEMO

Changeset 6509 for branches/NERC


Ignore:
Timestamp:
2016-05-05T12:07:34+02:00 (8 years ago)
Author:
jpalmier
Message:

JPALM -- 05-05-2016 -- update the pre-industrial pCO2 value

-- increase the ideal-tracer relax time at the surface

Location:
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/IDTRA/trcsms_idtra.F90

    r6466 r6509  
    7777      !! 
    7878      INTEGER                ::   ji, jj, jn, jl, jk 
    79       REAL(wp)               ::   rlx   !! relaxation time (1 day) 
     79      REAL(wp)               ::   rlx                 !! relaxation time (1 day) 
    8080      !!---------------------------------------------------------------------- 
    8181      ! 
    8282      IF( nn_timing == 1 )  CALL timing_start('trc_sms_idtra') 
    8383      ! 
     84      rlx = 10./(60. * 60. * 24.)                              !! relaxation time (1/10 day) 
    8485      IF (kt == nittrc000) THEN 
    8586         IF(lwp) WRITE(numout,*) '   trcsms_idtra :' 
    8687         IF(lwp) WRITE(numout,*) '   ~~~~~~~~~~~~~~~~~' 
    8788         IF(lwp) WRITE(numout,*) '   - idtra decay factor : ', FDEC 
     89         IF(lwp) WRITE(numout,*) '   - relaxation time    : ', rlx 
    8890# if defined key_debug_medusa 
    8991         CALL flush(numout) 
     
    9395 
    9496         ! 
    95       rlx = 1/(60 * 60 * 24)                                   !! relaxation time (1 day) 
    9697      inv_idtra(:,:,:) = 0.0                                   !! init the inventory 
    9798      qtr_idtra(:,:,:) = 0.0                                   !! init the air-sea flux 
     
    106107         !! First, a crude version. will be much inproved later. 
    107108             qtr_idtra(ji,jj,jl)  = rlx * (1. - trb(ji,jj,1,jn)) * tmask(ji,jj,1) *   &  
    108                                   fse3t(ji,jj,1) / rdt             !! Air-sea Flux 
     109                                  fse3t(ji,jj,1)                  !! Air-sea Flux 
    109110 
    110111           !! DEBUG-TEST : Set flux equal to 0, see if it induces the pb we see in the MED   
  • branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90

    r6466 r6509  
    998998      endif 
    999999#  if defined key_axy_pi_co2 
    1000       f_pco2a = hist_pco2(1) 
     1000      !! f_pco2a = hist_pco2(1) 
     1001      f_pco2a = 284.725          !! OCMIP pre-industrial pCO2 
    10011002#  endif 
    10021003      !! IF(lwp) WRITE(numout,*) ' MEDUSA nyear     =', nyear 
Note: See TracChangeset for help on using the changeset viewer.