- Timestamp:
- 2016-04-11T17:54:56+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/IDTRA/trcsms_idtra.F90
r6295 r6466 74 74 !! 75 75 !!---------------------------------------------------------------------- 76 INTEGER, INTENT( in ) :: kt ! ocean time-step index 77 !! 78 INTEGER :: ji, jj, jn, jl, jk 76 INTEGER, INTENT( in ) :: kt ! ocean time-step index 77 !! 78 INTEGER :: ji, jj, jn, jl, jk 79 REAL(wp) :: rlx !! relaxation time (1 day) 79 80 !!---------------------------------------------------------------------- 80 81 ! … … 92 93 93 94 ! 94 inv_idtra(:,:,:) = 0.0 !! init the inventory 95 qtr_idtra(:,:,:) = 0.0 !! init the air-sea flux 95 rlx = 1/(60 * 60 * 24) !! relaxation time (1 day) 96 inv_idtra(:,:,:) = 0.0 !! init the inventory 97 qtr_idtra(:,:,:) = 0.0 !! init the air-sea flux 96 98 DO jl = 1, jp_idtra 97 99 jn = jp_idtra0 + jl - 1 … … 103 105 104 106 !! First, a crude version. will be much inproved later. 105 qtr_idtra(ji,jj,jl) = (1. - trb(ji,jj,1,jn)) * tmask(ji,jj,1) * &107 qtr_idtra(ji,jj,jl) = rlx * (1. - trb(ji,jj,1,jn)) * tmask(ji,jj,1) * & 106 108 fse3t(ji,jj,1) / rdt !! Air-sea Flux 109 110 !! DEBUG-TEST : Set flux equal to 0, see if it induces the pb we see in the MED 111 !! qtr_idtra(ji,jj,jl) = 0.0 107 112 ENDDO 108 113 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.