Version 16 (modified by ctlod, 12 years ago) (diff)

Last edited Timestamp?

Sensitivity experiments performed with LIM 2.0


Abstract for the simulations performed over 10/500 or 1000 years (last ones are the first presented):

  • for PHC_40? based on tag nemo_v3_3 :
    - nn_sssr=2
    - alpha_precip=alpha_runoff=1. & no more +10% applied on evaporation
    - bug on sbcrnf.90 corrected
    

  • for PHC_42? based on PHC_16 but :
    - sea-ice albedo is set to 0.53 instead 0.50
    
  • for PHC_39? based on PHC_16 but :
    - Reduce the SSS damping parameter from 864mm/day to 166 mm/day
    
  • for PHC_38? based on PHC_16 but :
    - Remove the SSS damping under Sea-ice area using frld field in sbcssr.F90
    

  • for PHC_41? based on PHC_37 but :
    - damping term coefficient deds = -166mm/day instead -864mm/day 
    
  • for PHC_41? based on PHC_37 but :
    - tapering the damping term erp o -6 or 6 mm/day 
    
  • for PHC_37? based on PHC_36 but :
    - set nn_sssr = 2 instead 1 to take into account concentration/dilution effect through emps array.
    
  • for PHC_36? based on PHC_35 but :
    - remove the -9% on precipitation & runoff and +10% on evaporation field since the bug on sbcrnf.F90 is corrected
    
  • for PHC_35? based on PHC_16 but :
    - bug correction in sbcrnf.F90 module concerning missing test MOD(kt-1,nn_fsbc).
    
  • for PHC_34? based on PHC_16 but :
    - set now nn_sssr = 2 instead 1 and correct bu into sbcssr.F90 using (nn_sssr >= 1) instead (nn_sssr == 1)
    
  • for PHC_32? based on PHC_16 but :
    - apply air temperature correction on DFS4 data between 65°S & 80°S following differences with CLIO data.
    - The PHC_33 experience has a lower correction on air temperature than PHC_32.
    
  • for PHC_31? based on PHC_16 but :
    - set minimum of leads (acrit variable) to 2% in both hemispheres in the sea-ice namelist instead default value of 1.1O^-6
    


  • for PHC_16? based on PHC_15 but :
    - apply -9% on precipitation & runoff instead of -8% in PHC_15 experiment
    
  • for PHC_15?: based on PHC_11 but :
    just use "good" input fields for K1 and M2 waves
    
  • for PHC_12 based on PHC_11 but :
    few tests to improve transport of few straits:
    - fmask = 4.0 at the output of the Baltic sea
    - fmask = 0.5 at Gibraltar
    - fmask = 1.0 at Bab El Mandeb & add an identical value at one corner of the strait
    - one specificity of the Sound strait is an annual mean output flow, so to catch it in outputs,
    few modifications have been done in the diafwb.F90 module. further more an error on the location
    of the Sound strait has been also corrected
    - 
    
  • for PHC_11: based on PHC_10 but includes both:
    it seems that input fields were not the good ones
    - global Kz tides 
    - increased Kz tides in the ITF area
    
  • for PHC_10?: based on PHC_9 but includes both:
    - 3 waves lenghts penetration (RGB) through the traqsr.F90 module 
    - ocean color: climatology of a 2D chlorophyll field
    
  • for PHC_9?: based on the PHC_8 experiment
    change the way to compute wind stress & wind module:
    - wind module is now read and corresponds to ||U||clim, i.e. the one computed using 6hourly components 
    not the one using climatological components Uclim/Vclim 
    - wind stress: (||U||.u)clim & (||U||.v)clim fields are directly read instead computed. the drag coefficient is the one 
    computed using on the wind module read just above.
    - in doing this we don't substract anymore relative velocity of the ocean/sea-ice in the computation of the wind module 
    and the wind stress
    
  • for PHC_4?: based on the PHC_3 experiment
    just remove the TKE penetration in setting nn_etau= 0 in namelist:
    
  • from PHC_1 to PHC_8 :
    all these runs have been used for:
    - freshwater budget adjustement in playing with new namelist parameters alpha_precip, alpha_runoff
    - Gibraltar, Bab El Mandeb & Sound sraits transports adjustement ( based on scale factors or/and bondary condition fmask )
    - they all used the following computation for the wind module:  we add the wind module variance difference 
    (||U||clim)^2 - (Uclim^2 + Vclim^2) to the computation of the wind module where ||U||clim is the climatological wind module 
    and Uclim/Vclim are the climatological wind component
    


Results from the 10th year:

Runs PHC_1 PHC_2 PHC_3 PHC_5 PHC_6 PHC_7 PHC_8 PHC_9 PHC_10 PHC_11 PHC_15 PHC_16
Precip / -10% -12% / -8% -10% -8% -8% -8% -8% -8% -9%
Evap / / / +10% +5% +10% +10% +10% +10% +10% +10% +10%
Empold (Sv) +1.88 +0.26 -0.06 +1.37 +0.31 -0.25 +0.07 +0.13 +0.1 +0.12 +0.13 -0.02
SST/rt PHC3 / / / ++ + ++ ++ ?? ?? ?? ?? ??
SSH drift (m) / / -0.006 +0.12 +0.028 -0.022 +0.006 +0.01 +0.008 +0.01 +0.01 -0.002


Attachments (1)

Download all attachments as: .zip