Opened 3 months ago

Last modified 3 months ago

#2491 assigned Defect

output time step in istate

Reported by: rblod Owned by: rblod
Priority: low Milestone: Unscheduled
Component: TOP Version: trunk
Severity: minor Keywords:
Cc:

Description

Context

When outputting the initial state (nn_istate=1), the velocities doesn't correspond to the initial state (easy to see when starting from rest, velocities in output.ini.nc are not zero)

Analysis

diawri_state is called inside dia_wri so writes the now time step (Nnn). Velocities at time step now are corrected in dynspg_ts so doesn't correspond to the initial state anymore

Recommendation

Inside dia_wri, call diawri_state with before time step as argument : mod(kmm+1,3)+1 should do it
If assigned to me, I can do it …

Commit History (0)

(No commits)

Change History (3)

comment:1 Changed 3 months ago by smasson

  • Owner changed from systeam to rblod
  • Status changed from new to assigned

enjoy!

comment:2 Changed 3 months ago by jchanut

you can also simply revert to correct initial barotropic velocities before writing:

      DO jk = 1, jpkm1
         puu(:,:,jk,Kmm) = ( puu(:,:,jk,Kmm) - un_adv(:,:)*r1_hu(:,:,Kmm) + puu_b(:,:,Kmm) ) * umask(:,:,jk)
         pvv(:,:,jk,Kmm) = ( pvv(:,:,jk,Kmm) - vn_adv(:,:)*r1_hv(:,:,Kmm) + pvv_b(:,:,Kmm) ) * vmask(:,:,jk)
      END DO

comment:3 Changed 3 months ago by clem

This is also the case for sea-ice. It writes the state of sea ice at the end of the 1st time step…

Note: See TracTickets for help on using tickets.