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/nemogcm.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/nemogcm.F90

    r10888 r11355  
    150150      CALL mpp_max( 'nemogcm', nstop ) 
    151151 
     152      IF(narea == 419) THEN 
     153         OPEN(9419, file='crash_stats_419.out', status='replace') 
     154         WRITE(9419,*) 'Opening output file' 
     155      ENDIF 
     156 
     157      IF(narea == 68) THEN 
     158         OPEN(968, file='crash_stats_68.out', status='replace') 
     159         WRITE(968,*) 'Opening output file' 
     160      ENDIF 
     161 
    152162      IF(lwp) WRITE(numout,cform_aaa)   ! Flag AAAAAAA 
    153163 
     
    156166      !                            !-----------------------! 
    157167      istp = nit000 
     168 
     169      IF(narea == 419) THEN 
     170         WRITE(9419,*) 'max sshn nemogcm before any steps = ',MAXVAL(  ABS( sshn(:,:) )  ) 
     171         WRITE(9419,*) 'max ssha = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     172         WRITE(9419,*) 'max emp_b = ',MAXVAL(  ABS( emp_b(:,:) )  ) 
     173         WRITE(9419,*) 'max emp = ',MAXVAL(  ABS( emp(:,:) )  ) 
     174      ENDIF 
     175 
     176      ! Initialise emp to zero 
     177      emp(:,:) = 0.0 
     178 
    158179      ! 
    159180#if defined key_c1d 
     
    218239      !                            !------------------------! 
    219240      IF(lwp) WRITE(numout,cform_aaa)        ! Flag AAAAAAA 
     241      CLOSE(unit=9419) 
     242      CLOSE(unit=968) 
    220243      ! 
    221244      IF( nstop /= 0 .AND. lwp ) THEN        ! error print 
Note: See TracChangeset for help on using the changeset viewer.