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 6466 for branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/IDTRA – NEMO

Ignore:
Timestamp:
2016-04-11T17:54:56+02:00 (8 years ago)
Author:
jpalmier
Message:

JPALM -- 11-04-2016 -- add dust deposition input through namelist

-- relax time to IDTRA surface flux

File:
1 edited

Legend:

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

    r6295 r6466  
    7474      !!               
    7575      !!---------------------------------------------------------------------- 
    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) 
    7980      !!---------------------------------------------------------------------- 
    8081      ! 
     
    9293 
    9394         ! 
    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 
    9698      DO jl = 1, jp_idtra 
    9799         jn = jp_idtra0 + jl - 1 
     
    103105 
    104106         !! 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) *   &  
    106108                                  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 
    107112           ENDDO 
    108113         ENDDO 
Note: See TracChangeset for help on using the changeset viewer.