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.
ticket/0663_MLF (diff) – NEMO

Changes between Version 36 and Version 37 of ticket/0663_MLF


Ignore:
Timestamp:
2010-09-06T17:44:50+02:00 (14 years ago)
Author:
mlelod
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • ticket/0663_MLF

    v36 v37  
    111111      Heat content trends due to penetrative solar radiation (3D) are stored at now and before time steps in qsr_trd_hc_n and qsr_trd_hc_b. 
    112112 
    113    '''sbcmod : ''' Swap and restart of these before new forcing fields 
    114  
    115       Add qsr_trd_hc_b, sbc_trd_hc_b and sbc_trd_sc_b in swap and restart read/write. 
     113   '''sbcmod : '''  
     114 
     115      Correction of a little bug in restart "put":  when writing restart files, the utau_b, vtau_b, qns_b variables etc... contain utau, vtau, qns. 
    116116 
    117117   '''trasbc : ''' Changes for modified forcing integration 
    118118 
    119       sbc_trd_(h/s)c_n is aded to ta in forward Euler time stepping case 
     119      Swap and restart of sbc_trd_(h/s)c 
     120 
     121      sbc_trd_(h/s)c_n is added to ta in forward Euler time stepping case 
    120122 
    121123      0.5 * ( sbc_trd_(h/s)c_b + sbc_trd_(h/s)c_n ) is aded to "ta" in the leap-frog time stepping case 
     
    123125   '''traqsr : ''' Changes for modified forcing integration 
    124126 
     127      Swap and restart of qsr_trd_hc 
     128 
    125129      The 3D heat content tendency term of penetrative solar flux is stored in "qsr_trd_hc_n" in all cases. To that purpose, the division by fse3t in the computation of etot3 has been removed. It is done in the next step, when heat content tendency is converted to temperature tendency and added to ta. 
    126130 
     
    129133      Then only qsr_trd_hc_n variable is aded to ta in the forward Euler time stepping case and 0.5 * ( qsr_trd_hc_b + qsr_trd_hc_n ) is added to ta in the leap-frog time stepping case. 
    130134 
    131    '''sshwzv : ''' Sea surface elevation time stepping 
    132       comment for the forward Euler time stepping case: "In the forward Euler time stepping case, the same formulation as in the leap-frog case can be used because emp_b field is initialized with the vlaues of emp field. Hence, 0.5 * ( emp + emp_b ) = emp" 
     135   '''sshwzv : '''  
    133136 
    134137      surface weighted averaging for the computation of sshf_a and sshf_b in vvl case (as for sshu and sshv)