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

Timestamp:

2010-04-26T12:49:28+02:00 (14 years ago)

Author:

gm

Comment:

--

ticket/0663_MLF

@@ -50 +50 @@
 
    [[BR]]
+
    '''trasbc.F90''' : change in the heat and freshwater/salt fluxes 
       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.
@@ -77 +78 @@
 
 
-
-'''(2) The Asselin filter''' remove the forcing from the filter (dynnxt (?), tranxt, traqsr)
+   '''sshwzv.F90 : ''' applied the freshwater flux (emp) as the mean of 2 consecutive time-steps (kt+1/2 and kt-1/2)
+
+   [[BR]]
+
+
+
+'''(2) The Asselin filter''' remove the forcing from the filter (tranxt and sshwzv)
+      CAUTION : the sign of the correction applied to the asselin filter has to be checked.
+
+   [[BR]]
+
+
+   '''tranxt.F90 : ''' remove the forcing term from the Asselin filter
+
+      after the application of the classical Asselin filter, the part applied to the forcing term is removed :
+ 
+        at jk=1 :     - atfp * rdt * ( qns_b - qns )     on  temperature   (vvl case)    
+
+        at jk=1 :     - atfp * rdt * ( emps_b - emps )     on  salinity   (vvl case)    
+
+      if lk_vvl=F, the part applied to the concentration dilution term has also to be removed.
+
+      Solar penetration part :  
+
+      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.
+   [[BR]]
+
+
+   '''sshwzv.F90 : '''  a additional term appears in the filter (the last part of the RHS) :
+
+      sshn = ssha + atfp * ( sshb -2 sshn + ssha ) - atfp * rdt * ( emp_b - emp ) / rau0
+
+   [[BR]]
+
+   '''N.B. :'''  We choose not to remove the momentum forcing from the Asselin time filter. Therefore dynnxt.F90 is not changed.
 
 [[BR]]