Changeset 6509 for branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM
- Timestamp:
- 2016-05-05T12:07:34+02:00 (8 years ago)
- 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 77 77 !! 78 78 INTEGER :: ji, jj, jn, jl, jk 79 REAL(wp) :: rlx !! relaxation time (1 day)79 REAL(wp) :: rlx !! relaxation time (1 day) 80 80 !!---------------------------------------------------------------------- 81 81 ! 82 82 IF( nn_timing == 1 ) CALL timing_start('trc_sms_idtra') 83 83 ! 84 rlx = 10./(60. * 60. * 24.) !! relaxation time (1/10 day) 84 85 IF (kt == nittrc000) THEN 85 86 IF(lwp) WRITE(numout,*) ' trcsms_idtra :' 86 87 IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~~~~~~~' 87 88 IF(lwp) WRITE(numout,*) ' - idtra decay factor : ', FDEC 89 IF(lwp) WRITE(numout,*) ' - relaxation time : ', rlx 88 90 # if defined key_debug_medusa 89 91 CALL flush(numout) … … 93 95 94 96 ! 95 rlx = 1/(60 * 60 * 24) !! relaxation time (1 day)96 97 inv_idtra(:,:,:) = 0.0 !! init the inventory 97 98 qtr_idtra(:,:,:) = 0.0 !! init the air-sea flux … … 106 107 !! First, a crude version. will be much inproved later. 107 108 qtr_idtra(ji,jj,jl) = rlx * (1. - trb(ji,jj,1,jn)) * tmask(ji,jj,1) * & 108 fse3t(ji,jj,1) / rdt!! Air-sea Flux109 fse3t(ji,jj,1) !! Air-sea Flux 109 110 110 111 !! 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 998 998 endif 999 999 # 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 1001 1002 # endif 1002 1003 !! IF(lwp) WRITE(numout,*) ' MEDUSA nyear =', nyear
Note: See TracChangeset
for help on using the changeset viewer.