Ignore:
Timestamp:
2017-03-15T09:04:30+01:00 (4 years ago)
Author:
cbricaud
Message:

code cleaning and correct bug for wn computing in vvl case

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/TOP_SRC/oce_trc.F90

    r7521 r7795  
    147147 
    148148#if defined key_vvl  
    149    USE crs , ONLY :   e3t        =>  e3t_n_crs         !: vertical scale factors at t- 
    150    USE crs , ONLY :   e3u        =>  e3u_n_crs         !: vertical scale factors at u- 
    151    USE crs , ONLY :   e3v        =>  e3v_n_crs         !: vertical scale factors v- 
    152    USE crs , ONLY :   e3w        =>  e3w_n_crs         !: w-points (m) 
    153149   USE crs , ONLY :   e3t_b      =>  e3t_b_crs         !: vertical scale factors at t- 
    154    USE crs , ONLY :   e3u_b      =>  e3u_b_crs         !: vertical scale factors at u- 
    155    USE crs , ONLY :   e3v_b      =>  e3v_b_crs         !: vertical scale factors v- 
    156    USE crs , ONLY :   e3w_b      =>  e3w_b_crs         !: w-points (m) 
    157150   USE crs , ONLY :   e3t_n      =>  e3t_n_crs         !: vertical scale factors at t- 
    158151   USE crs , ONLY :   e3u_n      =>  e3u_n_crs         !: vertical scale factors at u- 
     
    160153   USE crs , ONLY :   e3w_n      =>  e3w_n_crs         !: w-points (m) 
    161154   USE crs , ONLY :   e3t_a      =>  e3t_a_crs         !: vertical scale factors at t- 
    162    USE crs , ONLY :   e3u_a      =>  e3u_a_crs         !: vertical scale factors at u- 
    163    USE crs , ONLY :   e3v_a      =>  e3v_a_crs         !: vertical scale factors v- 
    164    USE crs , ONLY :   e3w_a      =>  e3w_a_crs         !: w-points (m) 
    165    USE crs , ONLY :   gdept      =>  gdept_n_crs       !: depth of t-points (m) 
    166    USE crs , ONLY :   gdept_crs  =>  gdept_n_crs       !: depth of t-points (m) 
    167155   USE crs , ONLY :   gdept_n    =>  gdept_n_crs       !: depth of t-points (m) 
    168    USE crs , ONLY :   gdepw_n    =>  gdepw_n_crs       !: depth of t-points (m) 
    169    USE crs , ONLY :   e3t_max_crs => e3t_max_n_crs 
    170    USE crs , ONLY :   e3u_max_crs => e3u_max_n_crs 
    171    USE crs , ONLY :   e3v_max_crs => e3v_max_n_crs 
    172    USE crs , ONLY :   e3w_max_crs => e3w_max_n_crs 
    173156#else 
    174157   USE crs , ONLY :   e3t        =>  e3t_0_crs         !: vertical scale factors at t- 
     
    199182   USE crs , ONLY :   e3w_max_crs => e3w_max_0_crs 
    200183#endif 
    201    USE crs , ONLY :   e3t_0        =>  e3t_0_crs         !: vertical scale factors at t- 
    202    USE crs , ONLY :   e3u_0        =>  e3u_0_crs         !: vertical scale factors at t- 
    203    USE crs , ONLY :   e3v_0        =>  e3v_0_crs         !: vertical scale factors at t- 
    204    USE crs , ONLY :   e3w_0        =>  e3w_0_crs         !: vertical scale factors at t- 
    205  
    206    USE crs , ONLY     :   ff         =>  ff_crs         !: f-points (m) 
    207    USE crs , ONLY     :   gdept_0    =>  gdept_0_crs       !: depth of t-points (m) 
    208184   USE dom_oce , ONLY :   gdept_1d   =>  gdept_1d      !: depth of t-points (m) 
    209 #if defined key_zco 
    210    USE crs , ONLY :   gdept      =>  gdept_0_crs       !: depth of t-points (m) 
    211    USE crs , ONLY :   gdepw      =>  gdepw_crs       !: depth of t-points (m) 
    212 #endif 
    213185  !* masks, bathymetry * 
    214186   USE crs , ONLY :   mbkt       =>   mbkt_crs       !: vertical index of the bottom last T- ocean level 
    215    USE crs , ONLY :   mbku       =>   mbku_crs       !: vertical index of the bottom last U- ocean level 
    216    USE crs , ONLY :   mbkv       =>   mbkv_crs       !: vertical index of the bottom last V- ocean level 
    217187   USE crs , ONLY :   tmask_i    =>   tmask_i_crs    !: Interior mask at t-points 
    218188   USE crs , ONLY :   tmask      =>   tmask_crs      !: land/ocean mask at t-points 
    219    USE crs , ONLY :   umask      =>   umask_crs      !: land/ocean mask at u-points    
    220    USE crs , ONLY :   vmask      =>   vmask_crs      !: land/ocean mask at v-points  
    221    USE crs , ONLY :   fmask      =>   fmask_crs      !: land/ocean mask at f-points  
    222189 
    223190 !* ocean fields: here now and after fields * 
     
    226193   USE crs , ONLY :   wn      =>    wn_crs      !: vertical velocity (m s-1)   
    227194   USE crs , ONLY :   tsn     =>    tsn_crs     !: 4D array contaning ( tn, sn ) 
    228    USE crs , ONLY :   tsb     =>    tsb_crs     !: 4D array contaning ( tb, sb ) 
    229    USE crs , ONLY :   tsa     =>    tsa_crs     !: 4D array contaning ( ta, sa ) 
    230195   USE crs , ONLY :   rhop    =>    rhop_crs    !: potential volumic mass (kg m-3)  
    231    USE crs , ONLY :   rhd     =>    rhd_crs    !: in situ density anomalie rhd=(rho-rau0)/rau0 (no units) 
    232    USE crs , ONLY :   rn2b    =>    rb2_crs     !: in situ density anomalie rhd=(rho-rau0)/rau0 (no units) 
    233    USE crs , ONLY :   rab_n   =>    rab_crs_n     !: in situ density anomalie rhd=(rho-rau0)/rau0 (no units) 
    234196   USE crs , ONLY :   hdivn   =>    hdivn_crs   !: horizontal divergence (1/s) 
    235    USE crs , ONLY :   hdivb   =>    hdivb_crs   !: horizontal divergence (1/s) 
    236197   USE crs , ONLY :   sshb    =>    sshb_crs    !: sea surface height at t-point [m]    
    237198   USE crs , ONLY :   sshn    =>    sshn_crs    !: sea surface height at t-point [m]    
     
    239200 
    240201   !* surface fluxes * 
    241    USE crs , ONLY :   utau       =>    utau_crs       !: i-surface stress component 
    242    USE crs , ONLY :   vtau       =>    vtau_crs       !: j-surface stress component 
    243    USE crs , ONLY :   wndm       =>    wndm_crs       !: 10m wind speed  
    244202   USE crs , ONLY :   qsr        =>    qsr_crs        !: penetrative solar radiation (w m-2)   
    245203   USE crs , ONLY :   emp        =>    emp_crs        !: freshwater budget: volume flux               [Kg/m2/s] 
    246204   USE crs , ONLY :   emp_b      =>    emp_b_crs      !: freshwater budget: volume flux               [Kg/m2/s] 
    247    USE crs , ONLY :   sfx        =>    sfx_crs        !: freshwater budget: concentration/dillution   [Kg/m2/s] 
    248205   USE crs , ONLY :   fmmflx     =>    fmmflx_crs     !: freshwater budget: volume flux               [Kg/m2/s] 
    249206   USE crs , ONLY :   rnf        =>    rnf_crs        !: river runoff   [Kg/m2/s] 
     
    251208   USE crs , ONLY :   nk_rnf     =>    nk_rnf_crs     !: depth of runoff in model level 
    252209   USE crs , ONLY :   fr_i       =>    fr_i_crs       !: ice fraction (between 0 to 1) 
     210 
    253211   USE trcnam_trp , ONLY :  aht0     =>   rn_ahtrc_0        !: horizontal eddy diffusivity for tracers (m2/s) 
    254212   USE crs , ONLY :  ahtu     =>   ahtu_crs        !: lateral diffusivity coef. at u-points  
     
    258216   USE crs , ONLY :  r_fact_lap     =>  r_fact_lap_crs        !: enhanced zonal diffusivity coefficient 
    259217   USE ldftra_oce , ONLY :  rldf     =>   rldf 
     218 
    260219   USE crs , ONLY :  trc_i => trc_i_crs 
    261220   USE crs , ONLY :  trc_o => trc_o_crs 
     221 
    262222   USE crs , ONLY :   avt        =>   avt_crs         !: vert. diffusivity coef. at w-point for temp   
    263 #if defined key_zdfddm 
    264    USE crs , ONLY :   avs        =>   avs_crs         !: salinity vertical diffusivity coeff. at w-point 
    265 #endif 
    266  
    267    USE trc_oce, ONLY : nn_dttrc 
    268  
    269    USE crs , ONLY :   nmln        =>   nmln_crs        !: number of level in the mixed layer 
    270223   USE crs , ONLY :   hmld        =>   hmld_crs        !: mixing layer depth (turbocline) 
    271224   USE crs , ONLY :   hmlp        =>   hmlp_crs        !: mixed layer depth  (rho=rho0+zdcrit) (m) 
    272    USE crs , ONLY :   hmlpt       =>   hmlpt_crs       !: mixed layer depth at t-points (m) 
    273225 
    274226  !* direction of lateral diffusion * 
Note: See TracChangeset for help on using the changeset viewer.