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 – NEMO

Changeset 11367


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

More print statements including fix

Location:
NEMO/branches/UKMO/NEMO_4.0_GO8_coupled_iodef/src
Files:
11 edited

Legend:

Unmodified
Added
Removed

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

    r11082 r11367  
    107107 
    108108      CASE ( np_dynALL )           !==  all dynamical processes  ==! 
     109 
     110      IF(narea == 68) THEN 
     111         WRITE(968,*) 'in icedyn before ice_dyn_rhg: u_ice = ',SUM(u_ice) 
     112         WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 
     113         WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 
     114         WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 
     115         WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 
     116         WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 
     117         WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 
     118         WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 
     119      ENDIF 
     120 
    109121         ! 
    110122         CALL ice_dyn_rhg   ( kt )                                          ! -- rheology   
     123 
     124      IF(narea == 68) THEN 
     125         WRITE(968,*) 'in icedyn before ice_dyn_adv: u_ice = ',SUM(u_ice) 
     126         WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 
     127         WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 
     128         WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 
     129         WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 
     130         WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 
     131         WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 
     132         WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 
     133      ENDIF 
     134 
    111135         CALL ice_dyn_adv   ( kt )                                          ! -- advection of ice 
     136 
     137      IF(narea == 68) THEN 
     138         WRITE(968,*) 'a_i after ice_dyn_adv = ',a_i(54,1) 
     139      ENDIF 
     140 
    112141         CALL ice_dyn_rdgrft( kt )                                          ! -- ridging/rafting  
     142 
     143      IF(narea == 68) THEN 
     144         WRITE(968,*) 'a_i after ice_dyn_rdgrft = ',a_i(54,1) 
     145      ENDIF 
     146 
    113147         CALL ice_cor       ( kt , 1 )                                      ! -- Corrections 
     148 
     149      IF(narea == 68) THEN 
     150         WRITE(968,*) 'a_i after ice_cor = ',a_i(54,1) 
     151      ENDIF 
     152 
    114153         ! 
    115154      CASE ( np_dynRHGADV  )       !==  no ridge/raft & no corrections ==! 

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

    r11082 r11367  
    7575         WRITE(numout,*) '~~~~~~~~~~~' 
    7676      ENDIF 
     77 
     78      IF(narea == 68) THEN 
     79         WRITE(968,*) 'in icedyn_adv: u_ice = ',SUM(u_ice) 
     80         WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 
     81         WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 
     82         WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 
     83         WRITE(968,*) 'h_ip(54,0:2) = ',SUM(h_ip) 
     84         WRITE(968,*) 'ato_i(54,0:2) = ',SUM(ato_i) 
     85         WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 
     86         WRITE(968,*) 'v_s(54,0:2) = ',SUM(v_s) 
     87         WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 
     88         WRITE(968,*) 'oa_i(54,0:2) = ',SUM(oa_i) 
     89         WRITE(968,*) 'a_i(54,0:2) = ',SUM(a_i) 
     90         WRITE(968,*) 'a_ip(54,0:2) = ',SUM(a_ip) 
     91         WRITE(968,*) 'v_ip(54,0:2) = ',SUM(v_ip) 
     92         WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 
     93         WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 
     94      ENDIF 
     95 
    7796      ! 
    7897      !---------------! 

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

    r10888 r11367  
    8585         CALL lbc_lnk_multi( 'icesbc', utau_ice, 'U', -1., vtau_ice, 'V', -1. ) 
    8686      ENDIF 
     87 
     88      IF(narea == 68) THEN 
     89         WRITE(968,*) 'utau_ice in icesbc = ',utau_ice(54,1) 
     90         WRITE(968,*) 'vtau_ice in icesbc = ',vtau_ice(54,1) 
     91      ENDIF 
     92 
    8793      ! 
    8894      IF( ln_timing )   CALL timing_stop('ice_sbc') 

  • 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 

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

    r11355 r11367  
    125125      ENDIF 
    126126      CALL lbc_lnk( 'icethd', zfric, 'T',  1. ) 
     127 
     128      IF(narea == 68) THEN 
     129         WRITE(968,*) 'sum qlead in icethd before loop = ',SUM(qlead) 
     130         WRITE(968,*) 'sum tmask = ',SUM(tmask) 
     131         WRITE(968,*) 'sum at_i = ',SUM(at_i) 
     132         WRITE(968,*) 'sum at_i_b = ',SUM(at_i_b) 
     133         WRITE(968,*) 'sum qsr_oce = ',SUM(qsr_oce) 
     134         WRITE(968,*) 'sum frq_m = ',SUM(frq_m) 
     135         WRITE(968,*) 'sum qns_oce = ',SUM(qns_oce) 
     136         WRITE(968,*) 'sum qemp_oce = ',SUM(qemp_oce) 
     137         WRITE(968,*) 'sum e3t_m = ',SUM(e3t_m) 
     138         WRITE(968,*) 'sum t_bo = ',SUM(t_bo) 
     139         WRITE(968,*) 'sum sst_m = ',SUM(sst_m) 
     140         WRITE(968,*) 'sum zfric = ',SUM(zfric) 
     141      ENDIF 
    127142      ! 
    128143      !--------------------------------------------------------------------! 
    … …  
    174189         END DO 
    175190      END DO 
     191 
     192      IF(narea == 68) THEN 
     193         WRITE(968,*) 'sum qlead in icethd after loop = ',SUM(qlead) 
     194      ENDIF 
    176195       
    177196      ! In case we bypass open-water ice formation 

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

    r11355 r11367  
    115115 
    116116      at_i(:,:) = SUM( a_i, dim=3 ) 
     117 
     118      IF(narea == 68) THEN 
     119         WRITE(968,*) 'utau_ice in icethd_do = ',utau_ice(54,1) 
     120         WRITE(968,*) 'vtau_ice in icethd_do = ',vtau_ice(54,1) 
     121      ENDIF 
     122 
    117123      !------------------------------------------------------------------------------! 
    118124      ! 1) Collection thickness of ice formed in leads and polynyas 
    … …  
    207213      END DO 
    208214 
    209       IF(narea == 419) THEN 
    210          WRITE(9419,*) 'npti in icetd_do = ',npti 
    211          WRITE(9419,*) 'max qlead = ',MAXVAL(  qlead(:,:) ) 
    212          WRITE(9419,*) 'min qlead = ',MINVAL(  qlead(:,:) ) 
    213          WRITE(9419,*) 'qlead(48,31) = ', qlead(48,31) 
    214          WRITE(9419,*) 'tau_icebfr(48,31) = ',tau_icebfr(48,31) 
     215      IF(narea == 68) THEN 
     216         WRITE(968,*) 'npti in icetd_do = ',npti 
     217         WRITE(968,*) 'max qlead = ',MAXVAL(  qlead(:,:) ) 
     218         WRITE(968,*) 'min qlead = ',MINVAL(  qlead(:,:) ) 
     219         WRITE(968,*) 'qlead(48,31) = ', qlead(48,31) 
     220         WRITE(968,*) 'tau_icebfr(48,31) = ',tau_icebfr(48,31) 
     221         WRITE(968,*) 'sum wfx_opw = ',SUM( wfx_opw ) 
    215222      ENDIF 
    216223 
    … …  
    239246         CALL tab_2d_1d( npti, nptidx(1:npti), sss_1d    (1:npti) , sss_m      ) 
    240247 
    241          IF(narea == 419) THEN 
    242             WRITE(9419,*) 'wfx_opw(48,31) = ',wfx_opw(48,31) 
    243             WRITE(9419,*) 'max wfx_opw_1d icethd_do after to 1d = ',MAXVAL(  ABS( wfx_opw_1d(:) )  ) 
     248         IF(narea == 68) THEN 
     249            WRITE(968,*) 'wfx_opw(48,31) = ',wfx_opw(48,31) 
     250            WRITE(968,*) 'max wfx_opw_1d icethd_do after to 1d = ',MAXVAL(  ABS( wfx_opw_1d(:) )  ) 
     251            WRITE(968,*) 'sum wfx_opw_1d = ',SUM( wfx_opw_1d ) 
    244252         ENDIF 
    245253 
    … …  
    309317         END DO 
    310318 
    311          IF(narea == 419) THEN 
    312             WRITE(9419,*) 'max wfx_opw_1d icethd_do after calculation = ',MAXVAL(  ABS( wfx_opw_1d(:) ) ) 
    313             WRITE(9419,*) 'max zv_newice icethd_do after calculation = ',MAXVAL(  ABS( zv_newice(:) )  ) 
    314             WRITE(9419,*) 'rhoi = ',rhoi 
    315             WRITE(9419,*) 'r1_rdtice = ',r1_rdtice 
    316             WRITE(9419,*) 'max qlead_1d icethd_do after calculation = ',MAXVAL(  ABS( qlead_1d(:) )  ) 
    317             WRITE(9419,*) 'max ze_newice icethd_do after calculation = ',MAXVAL(  ABS( ze_newice(:) )  ) 
    318             WRITE(9419,*) 'max t_bo_1d icethd_do after calculation = ',MAXVAL(  ABS( t_bo_1d(:) )  ) 
     319         IF(narea == 68) THEN 
     320            WRITE(968,*) 'sum wfx_opw_1d icethd_do after calculation = ',SUM( wfx_opw_1d ) 
     321            WRITE(968,*) 'sum zv_newice icethd_do after calculation = ',SUM(  zv_newice) 
     322            WRITE(968,*) 'rhoi = ',rhoi 
     323            WRITE(968,*) 'r1_rdtice = ',r1_rdtice 
     324            WRITE(968,*) 'sum qlead_1d icethd_do after calculation = ',SUM( qlead_1d) 
     325            WRITE(968,*) 'sum ze_newice icethd_do after calculation = ',SUM( ze_newice) 
     326            WRITE(968,*) 'sum t_bo_1d icethd_do after calculation = ',SUM( t_bo_1d) 
    319327         ENDIF 
    320328          

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

    r11355 r11367  
    145145         !      &           + wfx_opw(ji,jj) + wfx_dyn(ji,jj) + wfx_res(ji,jj) + wfx_lam(ji,jj) + wfx_pnd(ji,jj) 
    146146         ! wfx_snw(ji,jj) = wfx_snw_sni(ji,jj) + wfx_snw_dyn(ji,jj) + wfx_snw_sum(ji,jj) 
     147         WRITE(968,*) 'sum wfx_opw iceupdate before loop = ',SUM( wfx_opw ) 
    147148         WRITE(968,*) 'max emp_oce iceupdate before loop = ',MAXVAL(  ABS( emp_oce(:,:) )  ) 
    148149         WRITE(968,*) 'max wfx_err_sub iceupdate before loop = ',MAXVAL(  ABS( wfx_err_sub(:,:) )  ) 

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

    r11355 r11367  
    13791379      ENDIF 
    13801380 
    1381       IF(narea == 419) THEN 
    1382          WRITE(9419,*) 'k_ice, kt, k_fsbc = ',k_ice, kt, k_fsbc 
    1383          WRITE(9419,*) 'MOD( kt-1, k_fsbc ) = ',MOD( kt-1, k_fsbc ) 
     1381      IF(narea == 68) THEN 
     1382         WRITE(968,*) 'k_ice, kt, k_fsbc = ',k_ice, kt, k_fsbc 
     1383         WRITE(968,*) 'MOD( kt-1, k_fsbc ) = ',MOD( kt-1, k_fsbc ) 
    13841384      ENDIF 
    13851385       
    … …  
    14021402         ENDIF 
    14031403 
    1404          IF(narea == 419) THEN 
    1405             WRITE(9419,*) 'max zemp after coupling = ',MAXVAL(  ABS( zemp(:,:) )  ) 
     1404         IF(narea == 68) THEN 
     1405            WRITE(968,*) 'max zemp after coupling = ',MAXVAL(  ABS( zemp(:,:) )  ) 
    14061406         ENDIF 
    14071407 
    … …  
    14111411         IF( srcv(jpr_cal)%laction )     zemp(:,:) = zemp(:,:) - frcv(jpr_cal)%z3(:,:,1) 
    14121412 
    1413          IF(narea == 419) THEN 
    1414             WRITE(9419,*) 'max zemp after runoff = ',MAXVAL(  ABS( zemp(:,:) )  ) 
    1415             WRITE(9419,*) 'max emp after runoff = ',MAXVAL(  ABS( emp(:,:) )  ) 
     1413         IF(narea == 68) THEN 
     1414            WRITE(968,*) 'max zemp after runoff = ',MAXVAL(  ABS( zemp(:,:) )  ) 
     1415            WRITE(968,*) 'max emp after runoff = ',MAXVAL(  ABS( emp(:,:) )  ) 
    14161416         ENDIF 
    14171417  
    … …  
    14261426         ENDIF 
    14271427 
    1428          IF(narea == 419) THEN 
    1429             WRITE(9419,*) 'max emp after conversion from zemp = ',MAXVAL(  ABS( emp(:,:) )  ) 
     1428         IF(narea == 68) THEN 
     1429            WRITE(968,*) 'max emp after conversion from zemp = ',MAXVAL(  ABS( emp(:,:) )  ) 
    14301430         ENDIF 
    14311431 
    … …  
    15961596          
    15971597      ENDIF 
     1598 
     1599      IF(narea == 68) THEN 
     1600         WRITE(968,*) 'Has OASIS sent data = ',(  nrcvinfo(itx) == OASIS_Rcv ) 
     1601         WRITE(968,*) 'srcv(jpr_itx1)%laction = ',srcv(jpr_itx1)%laction 
     1602         WRITE(968,*) 'frcv(jpr_itx1)%z3(54,1,1) = ',frcv(jpr_itx1)%z3(54,1,1) 
     1603         WRITE(968,*) 'frcv(jpr_ity1)%z3(54,1,1) = ',frcv(jpr_ity1)%z3(54,1,1) 
     1604         WRITE(968,*) 'p_taui in sbc_cpl_ice_tau = ',p_taui(54,1) 
     1605         WRITE(968,*) 'p_tauj in sbc_cpl_ice_tau = ',p_tauj(54,1) 
     1606         WRITE(968,*) 'max p_taui in sbc_cpl_ice_tau = ',MAXVAL(ABS(p_taui(:,:))) 
     1607         WRITE(968,*) 'max p_tauj in sbc_cpl_ice_tau = ',MAXVAL(ABS(p_tauj(:,:))) 
     1608      ENDIF 
     1609 
    15981610      ! 
    15991611   END SUBROUTINE sbc_cpl_ice_tau 
    … …  
    17531765      ENDIF 
    17541766 
    1755       IF(narea == 419) THEN 
    1756          WRITE(9419,*) 'emp_tot(48,31) = ',emp_tot(48,31) 
    1757          WRITE(9419,*) 'emp_ice(48,31) = ',emp_ice(48,31) 
    1758          WRITE(9419,*) 'emp_oce(48,31) = ',emp_oce(48,31) 
    1759          WRITE(9419,*) 'sprecip(48,31) = ',sprecip(48,31) 
    1760          WRITE(9419,*) 'tprecip(48,31) = ',tprecip(48,31) 
    1761          WRITE(9419,*) 'evap_ice(48,31) = ',evap_ice(48,31,:) 
    1762          WRITE(9419,*) 'devap_ice(48,31) = ',devap_ice(48,31,:) 
     1767      IF(narea == 68) THEN 
     1768         WRITE(968,*) 'sum emp_tot = ',SUM(emp_tot) 
     1769         WRITE(968,*) 'emp_ice = ',SUM(emp_ice) 
     1770         WRITE(968,*) 'emp_oce = ',SUM(emp_oce) 
     1771         WRITE(968,*) 'sprecip = ',SUM(sprecip) 
     1772         WRITE(968,*) 'tprecip = ',SUM(tprecip) 
     1773         WRITE(968,*) 'evap_ice = ',SUM(evap_ice) 
     1774         WRITE(968,*) 'devap_ice = ',SUM(devap_ice) 
    17631775      ENDIF 
    17641776 
    … …  
    18141826      CASE( 'oce only' )         ! the required field is directly provided 
    18151827         zqns_tot(:,:) = frcv(jpr_qnsoce)%z3(:,:,1) 
     1828         DO jl=1,jpl 
     1829            zqns_ice(:,:,jl) = zqns_tot(:,:) * a_i(:,:,jl) 
     1830         ENDDO 
    18161831      CASE( 'conservative' )     ! the required fields are directly provided 
    18171832         zqns_tot(:,:) = frcv(jpr_qnsmix)%z3(:,:,1) 
    … …  
    18561871      zqns_oce = 0._wp 
    18571872      WHERE( ziceld /= 0._wp )   zqns_oce(:,:) = ( zqns_tot(:,:) - SUM( a_i * zqns_ice, dim=3 ) ) / ziceld(:,:) 
     1873 
     1874      IF(narea == 68) THEN 
     1875         WRITE(968,*) 'sum zqns_oce in sbccpl = ',SUM(zqns_oce) 
     1876         WRITE(968,*) 'sum zqns_tot = ',SUM(zqns_tot) 
     1877         WRITE(968,*) 'sum zqns_ice = ',SUM(zqns_ice) 
     1878         WRITE(968,*) 'sum ziceld = ',SUM(ziceld) 
     1879      ENDIF 
    18581880 
    18591881      ! Heat content per unit mass of snow (J/kg) 
    … …  
    19041926      ENDIF 
    19051927 
    1906       IF(narea == 419) THEN 
    1907          WRITE(9419,*) 'qns_tot(48,31) = ',qns_tot(48,31) 
    1908          WRITE(9419,*) 'qns_oce(48,31) = ',qns_oce(48,31) 
    1909          WRITE(9419,*) 'qns_ice(48,31) = ',qns_ice(48,31,:) 
    1910          WRITE(9419,*) 'qevap_ice(48,31) = ',qevap_ice(48,31,:) 
    1911          WRITE(9419,*) 'qprec_ice(48,31) = ',qprec_ice(48,31) 
    1912          WRITE(9419,*) 'qemp_oce(48,31) = ',qemp_oce(48,31) 
    1913          WRITE(9419,*) 'qemp_ice(48,31) = ',qemp_ice(48,31) 
     1928      IF(narea == 68) THEN 
     1929         WRITE(968,*) 'qns_tot) = ',SUM(qns_tot) 
     1930         WRITE(968,*) 'qns_oce) = ',SUM(qns_oce) 
     1931         WRITE(968,*) 'qns_ice) = ',SUM(qns_ice) 
     1932         WRITE(968,*) 'qevap_ice) = ',SUM(qevap_ice) 
     1933         WRITE(968,*) 'qprec_ice) = ',SUM(qprec_ice) 
     1934         WRITE(968,*) 'qemp_oce) = ',SUM(qemp_oce) 
     1935         WRITE(968,*) 'qemp_ice) = ',SUM(qemp_ice) 
    19141936      ENDIF 
    19151937 
    … …  
    19541976      CASE( 'oce only' ) 
    19551977         zqsr_tot(:,:  ) = MAX( 0._wp , frcv(jpr_qsroce)%z3(:,:,1) ) 
    1956          zqsr_ice(:,:,:) = 0.0 
     1978         DO jl=1,jpl 
     1979            zqsr_ice(:,:,jl) = zqsr_tot(:,:) * a_i(:,:,jl) 
     1980         ENDDO 
    19571981      CASE( 'conservative' ) 
    19581982         zqsr_tot(:,:  ) = frcv(jpr_qsrmix)%z3(:,:,1) 
    … …  
    19812005            END DO 
    19822006 
    1983             IF(narea == 419) THEN 
    1984                WRITE(9419,*) 'zqsr_ice(48,31) oce and ice = ',zqsr_ice(48,31,:) 
     2007            IF(narea == 68) THEN 
     2008               WRITE(968,*) 'zqsr_ice) oce and ice = ',SUM(zqsr_ice) 
    19852009            ENDIF 
    19862010 
    … …  
    20262050      ENDIF 
    20272051 
    2028       IF(narea == 419) THEN 
    2029          WRITE(9419,*) 'qsr_tot(48,31) = ',qsr_tot(48,31) 
    2030          WRITE(9419,*) 'qsr_ice(48,31) = ',qsr_ice(48,31,:) 
     2052      IF(narea == 68) THEN 
     2053         WRITE(968,*) 'qsr_tot) = ',SUM(qsr_tot) 
     2054         WRITE(968,*) 'qsr_ice) = ',SUM(qsr_ice) 
    20312055      ENDIF 
    20322056 
    … …  
    20542078      ENDIF 
    20552079 
    2056       IF(narea == 419) THEN 
    2057          WRITE(9419,*) 'dqns_ice(48,31) = ',dqns_ice(48,31,:) 
     2080      IF(narea == 68) THEN 
     2081         WRITE(968,*) 'dqns_ice) = ',SUM(dqns_ice) 
    20582082      ENDIF 
    20592083 
    … …  
    20682092      END SELECT 
    20692093 
    2070       IF(narea == 419) THEN 
    2071          WRITE(9419,*) 'qml_ice(48,31) = ',qml_ice(48,31,:) 
    2072          WRITE(9419,*) 'qcn_ice(48,31) = ',qcn_ice(48,31,:) 
     2094      IF(narea == 68) THEN 
     2095         WRITE(968,*) 'qml_ice) = ',SUM(qml_ice) 
     2096         WRITE(968,*) 'qcn_ice) = ',SUM(qcn_ice) 
    20732097      ENDIF 
    20742098 
    … …  
    20942118      ENDIF 
    20952119 
    2096       IF(narea == 419) THEN 
    2097          WRITE(9419,*) 'qtr_ice_top(48,31) = ',qtr_ice_top(48,31,:) 
     2120      IF(narea == 68) THEN 
     2121         WRITE(968,*) 'qtr_ice_top) = ',SUM(qtr_ice_top) 
    20982122      ENDIF 
    20992123 

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

    r11355 r11367  
    3535#if defined key_si3 
    3636   USE icestp         ! surface boundary condition: SI3 sea-ice model 
     37   USE ice 
    3738#endif 
    3839   USE sbcice_cice    ! surface boundary condition: CICE sea-ice model 
    … …  
    414415      IF( ln_wave     )   CALL sbc_wave( kt )            ! surface waves 
    415416 
    416       IF(narea == 419) THEN 
    417          WRITE(9419,*) 'nsbc, jp_purecpl = ',nsbc, jp_purecpl 
    418          WRITE(9419,*) 'max emp before sbc_cpl_rcv = ',MAXVAL(  ABS( emp(:,:) )  ) 
     417      IF(narea == 68) THEN 
     418         WRITE(968,*) 'in sbcmod before sbc_cpl_rcv: u_ice = ',SUM(u_ice) 
     419         WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 
     420         WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 
     421         WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 
     422         WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 
     423         WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 
     424         WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 
     425         WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 
    419426      ENDIF 
    420427 
    … …  
    440447      ! 
    441448      !                                            !==  Misc. Options  ==! 
    442       IF(narea == 419) THEN 
    443          WRITE(9419,*) 'nn_ice = ',nn_ice 
    444          WRITE(9419,*) 'max emp before ice_stp = ',MAXVAL(  ABS( emp(:,:) )  ) 
     449      IF(narea == 68) THEN 
     450         WRITE(968,*) 'in sbcmod before ice_stp: u_ice = ',SUM(u_ice) 
     451         WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 
     452         WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 
     453         WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 
     454         WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 
     455         WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 
     456         WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 
     457         WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 
    445458      ENDIF 
    446459      ! 

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

    r11355 r11367  
    7575   USE sbc_oce , ONLY : lk_oasis 
    7676   USE wet_dry        ! Wetting and drying setting   (wad_init routine) 
     77   USE ice 
    7778#if defined key_top 
    7879   USE trcini         ! passive tracer initialisation 
    … …  
    167168      istp = nit000 
    168169 
    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(:,:) )  ) 
     170      IF(narea == 68) THEN 
     171         WRITE(968,*) 'in step before any steps: u_ice = ',SUM(u_ice) 
     172         WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 
     173         WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 
     174         WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 
     175         WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 
     176         WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 
     177         WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 
     178         WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 
    174179      ENDIF 
    175180 

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

    r11355 r11367  
    3939   ! 
    4040   USE iom              ! xIOs server 
     41   USE ice 
    4142 
    4243   IMPLICIT NONE 
    … …  
    9495      IF( ln_timing )   CALL timing_start('stp') 
    9596 
    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(:,:) )  ) 
     97      IF(narea == 68) THEN 
     98         WRITE(968,*) 'in step before update IO: u_ice = ',SUM(u_ice) 
     99         WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 
     100         WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 
     101         WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 
     102         WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 
     103         WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 
     104         WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 
     105         WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 
    101106      ENDIF 
    102107 
    … …  
    115120      IF( ln_crs         )   CALL iom_setkt( kstp - nit000 + 1, TRIM(cxios_context)//"_crs" )   ! tell IOM we are at time step kstp 
    116121 
    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      IF(narea == 68) THEN 
     123         WRITE(968,*) 'in step before sbc: u_ice = ',SUM(u_ice) 
     124         WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 
     125         WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 
     126         WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 
     127         WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 
     128         WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 
     129         WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 
     130         WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 
    122131      ENDIF 
    123132 
    … …  
    139148      IF( ln_sto_eos ) CALL sto_pts( tsn  )          ! Random T/S fluctuations 
    140149 
    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(:,:) )  ) 
     150      IF(narea == 68) THEN 
     151         WRITE(968,*) 'in step before ocean physics: u_ice = ',SUM(u_ice) 
     152         WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 
     153         WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 
     154         WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 
     155         WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 
     156         WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 
     157         WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 
     158         WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 
    146159      ENDIF 
    147160 
    … …  
    180193      IF( l_ldfdyn_time                    )   CALL ldf_dyn( kstp )       ! eddy viscosity coeff.  
    181194 
    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(:,:) )  ) 
     195      IF(narea == 68) THEN 
     196         WRITE(968,*) 'in step before dynamics: u_ice = ',SUM(u_ice) 
     197         WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 
     198         WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 
     199         WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 
     200         WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 
     201         WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 
     202         WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 
     203         WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 
    187204      ENDIF 
    188205 
    … …  
    241258               &         CALL dyn_asm_inc   ( kstp )  ! apply dynamics assimilation increment 
    242259 
    243       IF(narea == 419) THEN 
    244          WRITE(9419,*) 'max ssha step before bdy_dyn3d_dmp = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     260      IF(narea == 68) THEN 
     261         WRITE(968,*) 'in step before bdy_dyn3d: u_ice = ',SUM(u_ice) 
     262         WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 
     263         WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 
     264         WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 
     265         WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 
     266         WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 
     267         WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 
     268         WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 
    245269      ENDIF 
    246270 
    … …  
    287311                         CALL dyn_spg       ( kstp )  ! surface pressure gradient 
    288312 
    289  
    290       IF(narea == 419) THEN 
    291          WRITE(9419,*) 'max ssha step before div_hor = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     313      IF(narea == 68) THEN 
     314         WRITE(968,*) 'in step before div_hor: u_ice = ',SUM(u_ice) 
     315         WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 
     316         WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 
     317         WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 
     318         WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 
     319         WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 
     320         WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 
     321         WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 
    292322      ENDIF 
    293323 
    … …  
    307337                         CALL dyn_zdf       ( kstp )  ! vertical diffusion 
    308338 
    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(:,:) )  ) 
     339      IF(narea == 68) THEN 
     340         WRITE(968,*) 'in step before cool skin: u_ice = ',SUM(u_ice) 
     341         WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 
     342         WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 
     343         WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 
     344         WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 
     345         WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 
     346         WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 
     347         WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 
    312348      ENDIF 
    313349 
    … …  
    317353      IF ( ln_diurnal )  CALL stp_diurnal( kstp ) 
    318354 
    319       IF(narea == 419) THEN 
    320          WRITE(9419,*) 'max sshn step before diagnostics = ',MAXVAL(  ABS( sshn(:,:) )  ) 
    321          WRITE(9419,*) 'max ssha = ',MAXVAL(  ABS( ssha(:,:) )  ) 
     355      IF(narea == 68) THEN 
     356         WRITE(968,*) 'in step before diagnostics: u_ice = ',SUM(u_ice) 
     357         WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 
     358         WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 
     359         WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 
     360         WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 
     361         WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 
     362         WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 
     363         WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 
    322364      ENDIF 
    323365       
    … …  
    343385#endif 
    344386 
    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(:,:) )  ) 
     387      IF(narea == 68) THEN 
     388         WRITE(968,*) 'in step before active tracers: u_ice = ',SUM(u_ice) 
     389         WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 
     390         WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 
     391         WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 
     392         WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 
     393         WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 
     394         WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 
     395         WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 
    348396      ENDIF 
    349397 
    … …  
    356404         & ln_trainc )   CALL tra_asm_inc   ( kstp )  ! apply tracer assimilation increment 
    357405 
    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(:,:) )  ) 
     406      IF(narea == 68) THEN 
     407         WRITE(968,*) 'in step before tra_sbc: u_ice = ',SUM(u_ice) 
     408         WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 
     409         WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 
     410         WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 
     411         WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 
     412         WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 
     413         WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 
     414         WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 
    361415      ENDIF 
    362416 
    363417                         CALL tra_sbc       ( kstp )  ! surface boundary condition 
    364418 
    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(:,:) )  ) 
     419      IF(narea == 68) THEN 
     420         WRITE(968,*) 'in step before tra_adv: u_ice = ',SUM(u_ice) 
     421         WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 
     422         WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 
     423         WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 
     424         WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 
     425         WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 
     426         WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 
     427         WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 
    368428      ENDIF 
    369429 
    … …  
    389449      IF( ln_zdfnpc  )   CALL tra_npc       ( kstp )  ! update after fields by non-penetrative convection 
    390450 
    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(:,:) )  ) 
     451      IF(narea == 68) THEN 
     452         WRITE(968,*) 'in step before tra_nxt: u_ice = ',SUM(u_ice) 
     453         WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 
     454         WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 
     455         WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 
     456         WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 
     457         WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 
     458         WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 
     459         WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 
    394460      ENDIF 
    395461 
    … …  
    413479                         CALL tra_nxt       ( kstp )  ! finalize (bcs) tracer fields at next time step and swap 
    414480 
    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(:,:) )  ) 
     481      IF(narea == 68) THEN 
     482         WRITE(968,*) 'in step before dyn_nxt: u_ice = ',SUM(u_ice) 
     483         WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 
     484         WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 
     485         WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 
     486         WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 
     487         WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 
     488         WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 
     489         WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 
    418490      ENDIF 
    419491 
    420492                         CALL dyn_nxt       ( kstp )  ! finalize (bcs) velocities at next time step and swap (always called after tra_nxt) 
    421493 
    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(:,:) )  ) 
     494      IF(narea == 68) THEN 
     495         WRITE(968,*) 'in step before ssh_swp: u_ice = ',SUM(u_ice) 
     496         WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 
     497         WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 
     498         WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 
     499         WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 
     500         WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 
     501         WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 
     502         WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 
    425503      ENDIF 
    426504 
    … …  
    435513      ! 
    436514 
    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(:,:) )  ) 
     515      IF(narea == 68) THEN 
     516         WRITE(968,*) 'in step before dia_hsb: u_ice = ',SUM(u_ice) 
     517         WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 
     518         WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 
     519         WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 
     520         WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 
     521         WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 
     522         WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 
     523         WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 
    440524      ENDIF 
    441525 
    … …  
    450534      IF( ln_sto_eos )   CALL sto_rst_write( kstp )   ! write restart file for stochastic parameters 
    451535 
    452       IF(narea == 419) THEN 
    453          WRITE(9419,*) 'max sshn step before agrif = ',MAXVAL(  ABS( sshn(:,:) )  ) 
     536      IF(narea == 68) THEN 
     537         WRITE(968,*) 'in step before agrif: u_ice = ',SUM(u_ice) 
     538         WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 
     539         WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 
     540         WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 
     541         WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 
     542         WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 
     543         WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 
     544         WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 
    454545      ENDIF 
    455546 
    … …  
    468559      !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 
    469560 
    470       IF(narea == 419) THEN 
    471          WRITE(9419,*) 'max sshn step before stp_ctl = ',MAXVAL(  ABS( sshn(:,:) )  ) 
     561      IF(narea == 68) THEN 
     562         WRITE(968,*) 'in step before stp_ctl: u_ice = ',SUM(u_ice) 
     563         WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 
     564         WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 
     565         WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 
     566         WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 
     567         WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 
     568         WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 
     569         WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 
    472570      ENDIF 
    473571 
Note: See TracChangeset for help on using the changeset viewer.