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 11367 for NEMO/branches/UKMO/NEMO_4.0_GO8_coupled_iodef/src/ICE/icestp.F90 – NEMO

Ignore:
Timestamp:
2019-07-29T18:10:24+02:00 (5 years ago)
Author:
dancopsey
Message:

More print statements including fix

File:
1 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/UKMO/NEMO_4.0_GO8_coupled_iodef/src/ICE/icestp.F90

    r11355 r11367  
    146146                                        CALL store_fields             ! Store now ice values 
    147147 
    148       IF(narea == 419) THEN 
    149          WRITE(9419,*) 'max emp ice_stp before ice_sbc_tau = ',emp(48,31) 
    150          WRITE(9419,*) 'wfx_opw(48,31) = ',wfx_opw(48,31) 
    151          WRITE(9419,*) 'qlead(48,31) = ',qlead(48,31) 
     148      IF(narea == 68) THEN 
     149         WRITE(968,*) 'in icestp before ice_sbc_tau: u_ice = ',SUM(u_ice) 
     150         WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 
     151         WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 
     152         WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 
     153         WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 
     154         WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 
     155         WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 
     156         WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 
     157         WRITE(968,*) 'sum wfx_opw = ',SUM( wfx_opw ) 
    152158      ENDIF 
    153159         ! 
    … …  
    166172         ! 
    167173 
    168       IF(narea == 419) THEN 
    169          WRITE(9419,*) 'max emp ice_stp before ice_dyn = ',emp(48,31) 
    170          WRITE(9419,*) 'wfx_opw(48,31) = ',wfx_opw(48,31) 
    171          WRITE(9419,*) 'qlead(48,31) = ',qlead(48,31) 
     174      IF(narea == 68) THEN 
     175         WRITE(968,*) 'in icestp before ice_dyn: u_ice = ',SUM(u_ice) 
     176         WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 
     177         WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 
     178         WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 
     179         WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 
     180         WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 
     181         WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 
     182         WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 
     183         WRITE(968,*) 'sum wfx_opw = ',SUM( wfx_opw ) 
    172184      ENDIF 
    173185 
    … …  
    179191         ! 
    180192 
    181       IF(narea == 419) THEN 
    182          WRITE(9419,*) 'max emp ice_stp before ice_var_glo2eqv = ',emp(48,31) 
    183          WRITE(9419,*) 'wfx_opw(48,31) = ',wfx_opw(48,31) 
    184          WRITE(9419,*) 'qlead(48,31) = ',qlead(48,31) 
     193      IF(narea == 68) THEN 
     194         WRITE(968,*) 'max emp ice_stp before ice_var_glo2eqv = ',emp(54,1) 
     195         WRITE(968,*) 'sum wfx_opw = ',SUM( wfx_opw ) 
     196         WRITE(968,*) 'sum qlead = ',SUM(qlead) 
     197         WRITE(968,*) 'a_i(54,1) = ',a_i(54,1) 
    185198      ENDIF 
    186199 
    … …  
    204217         !------------------------------------------------------! 
    205218 
    206       IF(narea == 419) THEN 
    207          WRITE(9419,*) 'max emp ice_stp before ice_sbc_flx = ',emp(48,31) 
    208          WRITE(9419,*) 'wfx_opw(48,31) = ',wfx_opw(48,31) 
    209          WRITE(9419,*) 'qlead(48,31) = ',qlead(48,31) 
     219      IF(narea == 68) THEN 
     220         WRITE(968,*) 'max emp ice_stp before ice_sbc_flx = ',emp(54,1) 
     221         WRITE(968,*) 'sum wfx_opw = ',SUM( wfx_opw ) 
     222         WRITE(968,*) 'sum qlead = ',SUM(qlead) 
     223         WRITE(968,*) 'a_i(54,1) = ',a_i(54,1) 
    210224      ENDIF 
    211225 
    … …  
    215229         !----------------------------! 
    216230 
    217       IF(narea == 419) THEN 
    218          WRITE(9419,*) 'max emp ice_stp before ice_thd = ',emp(48,31) 
    219          WRITE(9419,*) 'wfx_opw(48,31) = ',wfx_opw(48,31) 
    220          WRITE(9419,*) 'qlead(48,31) = ',qlead(48,31) 
     231      IF(narea == 68) THEN 
     232         WRITE(968,*) 'max emp ice_stp before ice_thd = ',emp(54,1) 
     233         WRITE(968,*) 'sum wfx_opw = ',SUM( wfx_opw ) 
     234         WRITE(968,*) 'sum qlead = ',SUM(qlead) 
     235         WRITE(968,*) 'a_i(54,1) = ',a_i(54,1) 
    221236      ENDIF 
    222237 
    223238         IF( ln_icethd )                CALL ice_thd( kt )            ! -- Ice thermodynamics   
    224239 
    225       IF(narea == 419) THEN 
    226          WRITE(9419,*) 'max emp ice_stp before ice_cor = ',emp(48,31) 
    227          WRITE(9419,*) 'wfx_opw(48,31) = ',wfx_opw(48,31) 
     240      IF(narea == 68) THEN 
     241         WRITE(968,*) 'max emp ice_stp before ice_cor = ',emp(54,1) 
     242         WRITE(968,*) 'sum wfx_opw = ',SUM( wfx_opw ) 
     243         WRITE(968,*) 'a_i(54,1) = ',a_i(54,1) 
    228244      ENDIF 
    229245     
    … …  
    232248         ! 
    233249 
    234       IF(narea == 419) THEN 
    235          WRITE(9419,*) 'max emp ice_stp before ice_var_glo2eqv = ',emp(48,31) 
    236          WRITE(9419,*) 'wfx_opw(48,31) = ',wfx_opw(48,31) 
     250      IF(narea == 68) THEN 
     251         WRITE(968,*) 'max emp ice_stp before ice_var_glo2eqv = ',emp(54,1) 
     252         WRITE(968,*) 'sum wfx_opw = ',SUM( wfx_opw ) 
    237253      ENDIF 
    238254 
    … …  
    241257         ! 
    242258 
    243       IF(narea == 419) THEN 
    244          WRITE(9419,*) 'max emp ice_stp before ice_update_flx = ',emp(48,31) 
    245          WRITE(9419,*) 'wfx_opw(48,31) = ',wfx_opw(48,31) 
     259      IF(narea == 68) THEN 
     260         WRITE(968,*) 'max emp ice_stp before ice_update_flx = ',emp(54,1) 
     261         WRITE(968,*) 'sum wfx_opw = ',SUM( wfx_opw ) 
    246262      ENDIF 
    247263 
    … …  
    249265         ! 
    250266 
    251       IF(narea == 419) THEN 
    252          WRITE(9419,*) 'max emp ice_stp before ice_dia = ',emp(48,31) 
     267      IF(narea == 68) THEN 
     268         WRITE(968,*) 'max emp ice_stp before ice_dia = ',emp(54,1) 
    253269      ENDIF 
    254270 
    … …  
    261277         IF( ln_icectl )                CALL ice_ctl( kt )            ! -- alerts in case of model crash 
    262278 
    263       IF(narea == 419) THEN 
    264          WRITE(9419,*) 'max emp ice_stp at end = ',emp(48,31) 
     279      IF(narea == 68) THEN 
     280         WRITE(968,*) 'max emp ice_stp at end = ',emp(54,1) 
    265281      ENDIF 
    266282 
Note: See TracChangeset for help on using the changeset viewer.