Changeset 12489 for NEMO/trunk/src/OFF


Ignore:
Timestamp:
2020-02-28T16:55:11+01:00 (9 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.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/OFF/dtadyn.F90

    r12377 r12489  
    546546      ENDIF 
    547547 
    548       ssh(:,:,Kmm) = ssh(:,:,Kmm) + atfp * ( ssh(:,:,Kbb) - 2 * ssh(:,:,Kmm) + ssh(:,:,Kaa))   
     548      ssh(:,:,Kmm) = ssh(:,:,Kmm) + rn_atfp * ( ssh(:,:,Kbb) - 2 * ssh(:,:,Kmm) + ssh(:,:,Kaa))   
    549549 
    550550      !! Do we also need to time filter e3t?? 
     
    622622      !!---------------------------------------------------------------------- 
    623623      ! 
    624       z2dt = 2._wp * rdt 
     624      z2dt = 2._wp * rn_Dt 
    625625      ! 
    626626      zhdiv(:,:) = 0._wp 
     
    629629      END DO 
    630630      !                                                ! Sea surface  elevation time-stepping 
    631       pssha(:,:) = ( psshb(:,:) - z2dt * ( r1_rau0 * pemp(:,:)  + zhdiv(:,:) ) ) * ssmask(:,:) 
     631      pssha(:,:) = ( psshb(:,:) - z2dt * ( r1_rho0 * pemp(:,:)  + zhdiv(:,:) ) ) * ssmask(:,:) 
    632632      !                                                 !  
    633633      !                                                 ! After acale factors at t-points ( z_star coordinate ) 
Note: See TracChangeset for help on using the changeset viewer.