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 12475 for NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/icedyn_rhg.F90 – NEMO

Ignore:
Timestamp:
2020-02-27T12:27:50+01:00 (4 years ago)
Author:
dancopsey
Message:
  • Add more print statements.
  • Move away from using snow to ice diagnostics and use a new snow to pond one instead.
File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/icedyn_rhg.F90

    r11081 r12475  
    6363      ! controls 
    6464      IF( ln_timing    )   CALL timing_start('icedyn_rhg')                                                             ! timing 
     65      write(numout,*)'ice_dyn_rhg 1: u_ice = ',u_ice(3,4) 
    6566      IF( ln_icediachk )   CALL ice_cons_hsm(0, 'icedyn_rhg', rdiag_v, rdiag_s, rdiag_t, rdiag_fv, rdiag_fs, rdiag_ft) ! conservation 
    6667      ! 
     
    7778         END DO 
    7879      ENDIF 
     80      write(numout,*)'ice_dyn_rhg 2: u_ice = ',u_ice(3,4) 
    7981      ! 
    8082      !--------------! 
     
    8688         !                             !------------------------! 
    8789         CALL ice_dyn_rhg_evp( kt, stress1_i, stress2_i, stress12_i, shear_i, divu_i, delta_i ) 
     90         write(numout,*)'ice_dyn_rhg 3: u_ice = ',u_ice(3,4) 
    8891         !          
    8992      END SELECT 
     
    9194      IF( lrst_ice ) THEN                       !* write EVP fields in the restart file 
    9295         IF( ln_rhg_EVP )   CALL rhg_evp_rst( 'WRITE', kt ) 
     96         write(numout,*)'ice_dyn_rhg 4: u_ice = ',u_ice(3,4) 
    9397      ENDIF 
    9498      ! 
    9599      ! controls 
    96100      IF( ln_icediachk )   CALL ice_cons_hsm(1, 'icedyn_rhg', rdiag_v, rdiag_s, rdiag_t, rdiag_fv, rdiag_fs, rdiag_ft) ! conservation 
     101      write(numout,*)'ice_dyn_rhg 5: u_ice = ',u_ice(3,4) 
    97102      IF( ln_ctl       )   CALL ice_prt3D   ('icedyn_rhg')                                                             ! prints 
    98       IF( ln_timing    )   CALL timing_stop ('icedyn_rhg')                                                             ! timing 
     103      write(numout,*)'ice_dyn_rhg 6: u_ice = ',u_ice(3,4) 
     104      IF( ln_timing    )   CALL timing_stop ('icedyn_rhg')   
     105      write(numout,*)'ice_dyn_rhg 7: u_ice = ',u_ice(3,4)                                                           ! timing 
    99106      ! 
    100107   END SUBROUTINE ice_dyn_rhg 
Note: See TracChangeset for help on using the changeset viewer.