Ignore:
Timestamp:
2020-02-28T16:55:11+01:00 (11 months ago)
Author:
davestorkey
Message:

Preparation for new timestepping scheme #2390.
Main changes:

  1. Initial euler timestep now handled in stp and not in TRA/DYN routines.
  2. Renaming of all timestep parameters. In summary, the namelist parameter is now rn_Dt and the current timestep is rDt (and rDt_ice, rDt_trc etc).
  3. Renaming of a few miscellaneous parameters, eg. atfp → rn_atfp (namelist parameter used everywhere) and rau0 → rho0.

This version gives bit-comparable results to the previous version of the trunk.

Location:
NEMO/trunk/src/OCE/FLO
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/OCE/FLO/flo4rk.F90

    r12377 r12489  
    130130         ! computation of Runge-Kutta factor 
    131131         DO jfl = 1, jpnfl 
    132             zrkxfl(jfl,jind) = rdt*zufl(jfl) 
    133             zrkyfl(jfl,jind) = rdt*zvfl(jfl) 
    134             zrkzfl(jfl,jind) = rdt*zwfl(jfl) 
     132            zrkxfl(jfl,jind) = rn_Dt*zufl(jfl) 
     133            zrkyfl(jfl,jind) = rn_Dt*zvfl(jfl) 
     134            zrkzfl(jfl,jind) = rn_Dt*zwfl(jfl) 
    135135         END DO 
    136136         IF( jind /= 4 ) THEN 
  • NEMO/trunk/src/OCE/FLO/floblk.F90

    r12377 r12489  
    233233            ! test to know if the "age" of the float is not bigger than the  
    234234            ! time step 
    235             IF( zagenewfl(jfl) > rdt ) THEN 
    236                zttfl(jfl) = (rdt-zagefl(jfl)) / zvol 
    237                zagenewfl(jfl) = rdt 
     235            IF( zagenewfl(jfl) > rn_Dt ) THEN 
     236               zttfl(jfl) = (rn_Dt-zagefl(jfl)) / zvol 
     237               zagenewfl(jfl) = rn_Dt 
    238238            ENDIF 
    239239             
     
    340340         ifin = 1 
    341341         DO jfl = 1, jpnfl 
    342             IF( zagefl(jfl) < rdt )   ifin = 0 
     342            IF( zagefl(jfl) < rn_Dt )   ifin = 0 
    343343            tpifl(jfl) = zgifl(jfl) + 0.5 
    344344            tpjfl(jfl) = zgjfl(jfl) + 0.5 
     
    347347         ifin = 1 
    348348         DO jfl = 1, jpnfl 
    349             IF( zagefl(jfl) < rdt )   ifin = 0 
     349            IF( zagefl(jfl) < rn_Dt )   ifin = 0 
    350350            tpifl(jfl) = zgifl(jfl) + 0.5 
    351351            tpjfl(jfl) = zgjfl(jfl) + 0.5 
  • NEMO/trunk/src/OCE/FLO/flowri.F90

    r12377 r12489  
    122122               ztem(jfl) = ts(iafloc,ibfloc,icfl,jp_tem,Kmm) 
    123123               zsal (jfl) = ts(iafloc,ibfloc,icfl,jp_sal,Kmm) 
    124                zrho (jfl) = (rhd(iafloc,ibfloc,icfl)+1)*rau0 
     124               zrho (jfl) = (rhd(iafloc,ibfloc,icfl)+1)*rho0 
    125125 
    126126            ENDIF 
     
    142142            ztem(jfl) = ts(iafloc,ibfloc,icfl,jp_tem,Kmm) 
    143143            zsal(jfl) = ts(iafloc,ibfloc,icfl,jp_sal,Kmm) 
    144             zrho(jfl) = (rhd(iafloc,ibfloc,icfl)+1)*rau0 
     144            zrho(jfl) = (rhd(iafloc,ibfloc,icfl)+1)*rho0 
    145145           
    146146         ENDIF 
     
    245245            !------------------------------- 
    246246            irec =  INT( (kt-nn_it000+1)/nn_writefl ) +1 
    247             ztime = ( kt-nn_it000 + 1 ) * rdt 
     247            ztime = ( kt-nn_it000 + 1 ) * rn_Dt 
    248248 
    249249            CALL flioputv( numflo , 'time_counter', ztime , start=(/irec/) ) 
Note: See TracChangeset for help on using the changeset viewer.