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.
Changeset 11355 for NEMO/branches/UKMO/NEMO_4.0_GO8_coupled_iodef/src/OCE/step.F90 – NEMO

Ignore:
Timestamp:
2019-07-26T11:32:13+02:00 (5 years ago)
Author:
dancopsey
Message:

Add lots of print statements printing out values of files throughout the ocean_ice timestep.. New output will be in files crash_stats_419.out and crash_stats_68.out.

File:
1 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/UKMO/NEMO_4.0_GO8_coupled_iodef/src/OCE/step.F90

    r11110 r11355  
    9393      ! 
    9494      IF( ln_timing )   CALL timing_start('stp') 
     95 
     96      IF(narea == 419) THEN 
     97         WRITE(9419,*) 'max sshn step before update IO = ',MAXVAL(  ABS( sshn(:,:) )  ) 
     98         WRITE(9419,*) 'max ssha = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     99         WRITE(9419,*) 'max emp_b = ',MAXVAL(  ABS( emp_b(:,:) )  ) 
     100         WRITE(9419,*) 'max emp = ',MAXVAL(  ABS( emp(:,:) )  ) 
     101      ENDIF 
     102 
    95103      ! 
    96104      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    … …  
    107115      IF( ln_crs         )   CALL iom_setkt( kstp - nit000 + 1, TRIM(cxios_context)//"_crs" )   ! tell IOM we are at time step kstp 
    108116 
     117      IF(narea == 419) THEN 
     118         WRITE(9419,*) 'max sshn step before external forcing = ',MAXVAL(  ABS( sshn(:,:) )  ) 
     119         WRITE(9419,*) 'max ssha = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     120         WRITE(9419,*) 'max emp_b = ',MAXVAL(  ABS( emp_b(:,:) )  ) 
     121         WRITE(9419,*) 'max emp = ',MAXVAL(  ABS( emp(:,:) )  ) 
     122      ENDIF 
     123 
    109124      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    110125      ! Update external forcing (tides, open boundaries, and surface boundary condition (including sea-ice) 
    … …  
    115130                         CALL sbc     ( kstp )                   ! Sea Boundary Condition (including sea-ice) 
    116131 
     132 
     133 
     134 
    117135      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    118136      ! Update stochastic parameters and random T/S fluctuations 
    … …  
    120138      IF( ln_sto_eos ) CALL sto_par( kstp )          ! Stochastic parameters 
    121139      IF( ln_sto_eos ) CALL sto_pts( tsn  )          ! Random T/S fluctuations 
     140 
     141      IF(narea == 419) THEN 
     142         WRITE(9419,*) 'max sshn step before ocean physics = ',MAXVAL(  ABS( sshn(:,:) )  ) 
     143         WRITE(9419,*) 'max ssha = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     144         WRITE(9419,*) 'max emp_b = ',MAXVAL(  ABS( emp_b(:,:) )  ) 
     145         WRITE(9419,*) 'max emp = ',MAXVAL(  ABS( emp(:,:) )  ) 
     146      ENDIF 
    122147 
    123148      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    … …  
    155180      IF( l_ldfdyn_time                    )   CALL ldf_dyn( kstp )       ! eddy viscosity coeff.  
    156181 
     182      IF(narea == 419) THEN 
     183         WRITE(9419,*) 'max sshn step before ocean dynamics = ',MAXVAL(  ABS( sshn(:,:) )  ) 
     184         WRITE(9419,*) 'max ssha = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     185         WRITE(9419,*) 'max emp_b = ',MAXVAL(  ABS( emp_b(:,:) )  ) 
     186         WRITE(9419,*) 'max emp = ',MAXVAL(  ABS( emp(:,:) )  ) 
     187      ENDIF 
     188 
    157189      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    158190      !  Ocean dynamics : hdiv, ssh, e3, u, v, w 
    … …  
    160192 
    161193                            CALL ssh_nxt       ( kstp )  ! after ssh (includes call to div_hor) 
     194      IF(narea == 419) THEN 
     195         WRITE(9419,*) 'max ssha step before dom_vvl_sf_nxt = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     196         CALL FLUSH(9419) 
     197      ENDIF 
     198 
     199      IF(narea == 68) THEN 
     200         WRITE(968,*) 'max ssha step before dom_vvl_sf_nxt = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     201         CALL FLUSH(968) 
     202      ENDIF 
     203 
     204 
    162205      IF( .NOT.ln_linssh )  CALL dom_vvl_sf_nxt( kstp )  ! after vertical scale factors  
     206 
     207      IF(narea == 419) THEN 
     208         WRITE(9419,*) 'max ssha step before wzv = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     209      ENDIF 
     210 
    163211                            CALL wzv           ( kstp )  ! now cross-level velocity  
     212 
     213      IF(narea == 419) THEN 
     214         WRITE(9419,*) 'max ssha step before wAimp = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     215      ENDIF 
     216 
    164217      IF( ln_zad_Aimp )     CALL wAimp         ( kstp )  ! Adaptive-implicit vertical advection partitioning 
    165218                            CALL eos    ( tsn, rhd, rhop, gdept_n(:,:,:) )  ! now in situ density for hpg computation 
     219 
     220      IF(narea == 419) THEN 
     221         WRITE(9419,*) 'max ssha step before simplification = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     222      ENDIF 
     223 
    166224                             
    167225!!jc: fs simplification 
    … …  
    182240      IF(  lk_asminc .AND. ln_asmiau .AND. ln_dyninc )   & 
    183241               &         CALL dyn_asm_inc   ( kstp )  ! apply dynamics assimilation increment 
     242 
     243      IF(narea == 419) THEN 
     244         WRITE(9419,*) 'max ssha step before bdy_dyn3d_dmp = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     245      ENDIF 
     246 
    184247      IF( ln_bdy     )   CALL bdy_dyn3d_dmp ( kstp )  ! bdy damping trends 
    185248#if defined key_agrif 
    … …  
    187250               &         CALL Agrif_Sponge_dyn        ! momentum sponge 
    188251#endif 
     252 
     253      IF(narea == 419) THEN 
     254         WRITE(9419,*) 'max ssha step before dom_adv = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     255      ENDIF 
     256 
    189257                         CALL dyn_adv       ( kstp )  ! advection (vector or flux form) 
     258 
     259      IF(narea == 419) THEN 
     260         WRITE(9419,*) 'max ssha step before dyn_vor = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     261      ENDIF 
     262 
    190263                         CALL dyn_vor       ( kstp )  ! vorticity term including Coriolis 
     264 
     265      IF(narea == 419) THEN 
     266         WRITE(9419,*) 'max ssha step before dyn_ldft = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     267      ENDIF 
     268 
    191269                         CALL dyn_ldf       ( kstp )  ! lateral mixing 
     270 
     271      IF(narea == 419) THEN 
     272         WRITE(9419,*) 'max ssha step before dyn_osm = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     273      ENDIF 
     274 
    192275      IF( ln_zdfosm  )   CALL dyn_osm       ( kstp )  ! OSMOSIS non-local velocity fluxes 
     276 
     277      IF(narea == 419) THEN 
     278         WRITE(9419,*) 'max ssha step before dyn_hpg = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     279      ENDIF 
     280 
    193281                         CALL dyn_hpg       ( kstp )  ! horizontal gradient of Hydrostatic pressure 
     282 
     283      IF(narea == 419) THEN 
     284         WRITE(9419,*) 'max ssha step before dyn_spg = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     285      ENDIF 
     286 
    194287                         CALL dyn_spg       ( kstp )  ! surface pressure gradient 
     288 
     289 
     290      IF(narea == 419) THEN 
     291         WRITE(9419,*) 'max ssha step before div_hor = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     292      ENDIF 
    195293 
    196294                                                      ! With split-explicit free surface, since now transports have been updated and ssha as well 
    … …  
    201299         IF( ln_zad_Aimp )  CALL wAimp      ( kstp )  ! Adaptive-implicit vertical advection partitioning 
    202300      ENDIF 
     301 
     302      IF(narea == 419) THEN 
     303         WRITE(9419,*) 'max ssha step before dyn_zdf = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     304      ENDIF 
     305 
    203306       
    204307                         CALL dyn_zdf       ( kstp )  ! vertical diffusion 
     308 
     309      IF(narea == 419) THEN 
     310         WRITE(9419,*) 'max sshn step before cool skin = ',MAXVAL(  ABS( sshn(:,:) )  ) 
     311         WRITE(9419,*) 'max ssha = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     312      ENDIF 
    205313 
    206314      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    … …  
    208316      !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<       
    209317      IF ( ln_diurnal )  CALL stp_diurnal( kstp ) 
     318 
     319      IF(narea == 419) THEN 
     320         WRITE(9419,*) 'max sshn step before diagnostics = ',MAXVAL(  ABS( sshn(:,:) )  ) 
     321         WRITE(9419,*) 'max ssha = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     322      ENDIF 
    210323       
    211324      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    … …  
    230343#endif 
    231344 
     345      IF(narea == 419) THEN 
     346         WRITE(9419,*) 'max sshn step before active tracers = ',MAXVAL(  ABS( sshn(:,:) )  ) 
     347         WRITE(9419,*) 'max ssha = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     348      ENDIF 
     349 
    232350      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    233351      ! Active tracers                               
    … …  
    237355      IF(  lk_asminc .AND. ln_asmiau .AND. & 
    238356         & ln_trainc )   CALL tra_asm_inc   ( kstp )  ! apply tracer assimilation increment 
     357 
     358      IF(narea == 419) THEN 
     359         WRITE(9419,*) 'max sshn step before tra_sbc = ',MAXVAL(  ABS( sshn(:,:) )  ) 
     360         WRITE(9419,*) 'max ssha = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     361      ENDIF 
     362 
    239363                         CALL tra_sbc       ( kstp )  ! surface boundary condition 
     364 
     365      IF(narea == 419) THEN 
     366         WRITE(9419,*) 'max sshn step after tra_sbc = ',MAXVAL(  ABS( sshn(:,:) )  ) 
     367         WRITE(9419,*) 'max ssha = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     368      ENDIF 
     369 
    240370      IF( ln_traqsr  )   CALL tra_qsr       ( kstp )  ! penetrative solar radiation qsr 
    241371      IF( ln_trabbc  )   CALL tra_bbc       ( kstp )  ! bottom heat flux 
    … …  
    258388                         CALL tra_zdf       ( kstp )  ! vertical mixing and after tracer fields 
    259389      IF( ln_zdfnpc  )   CALL tra_npc       ( kstp )  ! update after fields by non-penetrative convection 
     390 
     391      IF(narea == 419) THEN 
     392         WRITE(9419,*) 'max sshn step before boundary conditions = ',MAXVAL(  ABS( sshn(:,:) )  ) 
     393         WRITE(9419,*) 'max ssha = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     394      ENDIF 
    260395 
    261396      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    … …  
    277412!!jc2: dynnxt must be the latest call. e3t_b are indeed updated in that routine 
    278413                         CALL tra_nxt       ( kstp )  ! finalize (bcs) tracer fields at next time step and swap 
     414 
     415      IF(narea == 419) THEN 
     416         WRITE(9419,*) 'max sshn step before dyn_nxt = ',MAXVAL(  ABS( sshn(:,:) )  ) 
     417         WRITE(9419,*) 'max ssha = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     418      ENDIF 
     419 
    279420                         CALL dyn_nxt       ( kstp )  ! finalize (bcs) velocities at next time step and swap (always called after tra_nxt) 
     421 
     422      IF(narea == 419) THEN 
     423         WRITE(9419,*) 'max sshn step before ssh_swp = ',MAXVAL(  ABS( sshn(:,:) )  ) 
     424         WRITE(9419,*) 'max ssha = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     425      ENDIF 
     426 
    280427                         CALL ssh_swp       ( kstp )  ! swap of sea surface height 
     428 
     429      IF(narea == 419) THEN 
     430         WRITE(9419,*) 'max sshn step before dom_vvl_sf_swp = ',MAXVAL(  ABS( sshn(:,:) )  ) 
     431         WRITE(9419,*) 'max ssha = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     432      ENDIF 
     433 
    281434      IF(.NOT.ln_linssh) CALL dom_vvl_sf_swp( kstp )  ! swap of vertical scale factors 
    282435      ! 
     436 
     437      IF(narea == 419) THEN 
     438         WRITE(9419,*) 'max sshn step before dia_hsb = ',MAXVAL(  ABS( sshn(:,:) )  ) 
     439         WRITE(9419,*) 'max ssha = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     440      ENDIF 
     441 
    283442      IF( ln_diahsb  )   CALL dia_hsb       ( kstp )  ! - ML - global conservation diagnostics 
    284443 
    … …  
    291450      IF( ln_sto_eos )   CALL sto_rst_write( kstp )   ! write restart file for stochastic parameters 
    292451 
     452      IF(narea == 419) THEN 
     453         WRITE(9419,*) 'max sshn step before agrif = ',MAXVAL(  ABS( sshn(:,:) )  ) 
     454      ENDIF 
     455 
    293456#if defined key_agrif 
    294457      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    … …  
    304467      ! Control 
    305468      !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 
     469 
     470      IF(narea == 419) THEN 
     471         WRITE(9419,*) 'max sshn step before stp_ctl = ',MAXVAL(  ABS( sshn(:,:) )  ) 
     472      ENDIF 
     473 
    306474                         CALL stp_ctl      ( kstp, indic ) 
    307475                          
Note: See TracChangeset for help on using the changeset viewer.