Changes between Version 8 and Version 9 of ticket/0663_MLF


Ignore:
Timestamp:
2010-04-26T12:49:28+02:00 (11 years ago)
Author:
gm
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • ticket/0663_MLF

    v8 v9  
    5050 
    5151   [[BR]] 
     52 
    5253   '''trasbc.F90''' : change in the heat and freshwater/salt fluxes  
    5354      The thermohalin fluxes (qns, emp, emps) are now applied as the mean of 2 consecutive time-steps (kt+1/2 and kt-1/2). This involes the following changes : (1) incorporate the time mean of the forcing ; (2) put the IF-ELSE-ENDIF outside the DO-loop for vector optimisation ; (3) remove the key_zco optimisation (increase readability) as zco is less and less used ; (4) correct a bug in vvl case, where salt flux associated with freezing/melting were not taken into account. 
     
    7778 
    7879 
    79  
    80 '''(2) The Asselin filter''' remove the forcing from the filter (dynnxt (?), tranxt, traqsr) 
     80   '''sshwzv.F90 : ''' applied the freshwater flux (emp) as the mean of 2 consecutive time-steps (kt+1/2 and kt-1/2) 
     81 
     82   [[BR]] 
     83 
     84 
     85 
     86'''(2) The Asselin filter''' remove the forcing from the filter (tranxt and sshwzv) 
     87      CAUTION : the sign of the correction applied to the asselin filter has to be checked. 
     88 
     89   [[BR]] 
     90 
     91 
     92   '''tranxt.F90 : ''' remove the forcing term from the Asselin filter 
     93 
     94      after the application of the classical Asselin filter, the part applied to the forcing term is removed : 
     95  
     96        at jk=1 :     - atfp * rdt * ( qns_b - qns )     on  temperature   (vvl case)     
     97 
     98        at jk=1 :     - atfp * rdt * ( emps_b - emps )     on  salinity   (vvl case)     
     99 
     100      if lk_vvl=F, the part applied to the concentration dilution term has also to be removed. 
     101 
     102      Solar penetration part :   
     103 
     104      CAUTION 1 : the case where the light penetration is provided by the bio-model ( i.e. lk_qsr_bio = ln_qsr_bio =.true. is not yet treated. It probably requires changes in p4zopt.F90 and trcopt.F90. Not completely obvious since this depends on the TOP time-stepping strategy. 
     105   [[BR]] 
     106 
     107 
     108   '''sshwzv.F90 : '''  a additional term appears in the filter (the last part of the RHS) : 
     109 
     110      sshn = ssha + atfp * ( sshb -2 sshn + ssha ) - atfp * rdt * ( emp_b - emp ) / rau0 
     111 
     112   [[BR]] 
     113 
     114   '''N.B. :'''  We choose not to remove the momentum forcing from the Asselin time filter. Therefore dynnxt.F90 is not changed. 
    81115 
    82116[[BR]]