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

Changeset 8892


Ignore:
Timestamp:
2017-12-04T15:09:10+01:00 (6 years ago)
Author:
frrh
Message:

Commit updates with debugging write statements.

Location:
branches/UKMO/dev_r8126_LIM3_couple/NEMOGCM/NEMO
Files:
11 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r8126_LIM3_couple/NEMOGCM/NEMO/LIM_SRC_3/ice.F90

    r8879 r8892  
    480480      ii = ii + 1 
    481481      ALLOCATE( u_ice_b(jpi,jpj) , v_ice_b(jpi,jpj) , at_i_b(jpi,jpj) , STAT=ierr(ii) ) 
    482        
     482at_i_b(:,:) = 0.0       
     483u_ice_b(:,:) = 0.0       
     484v_ice_b(:,:) = 0.0       
     485 
    483486      ! * Ice thickness distribution variables 
    484487      ii = ii + 1 
  • branches/UKMO/dev_r8126_LIM3_couple/NEMOGCM/NEMO/LIM_SRC_3/icecor.F90

    r8879 r8892  
    103103      !                             !----------------------------------------------------- 
    104104      IF ( jpl > 1 )   CALL ice_itd_reb( kt ) 
    105  
     105write(numout,*) "RSRH out of ice_itd_reb in icecor", kt ; flush(numout) 
    106106      !                             !----------------------------------------------------- 
    107107      CALL ice_var_zapsmall         !  Zap small values                                  ! 
     
    133133      CASE( 1 )                        !--- dyn trend diagnostics 
    134134         ! 
     135write(numout,*) "RSRH GG icecor CASE 1", kt ; flush(numout) 
    135136!!gm   here I think the number of ice cat is too small to use a SUM instruction... 
    136137         DO jj = 1, jpj 
     
    152153      CASE( 2 )                        !--- thermo trend diagnostics & ice aging 
    153154         ! 
     155write(numout,*) "RSRH GG icecor CASE 2", kt ; flush(numout) 
    154156         oa_i(:,:,:) = oa_i(:,:,:) + a_i(:,:,:) * rdt_ice   ! ice natural aging incrementation 
    155157         ! 
     
    173175         ! 
    174176      END SELECT 
     177write(numout,*) "RSRH HH icecor", kt ; flush(numout) 
    175178      ! 
    176179      ! controls 
     
    179182      IF( ln_icectl .AND. kn == 2 )   CALL ice_prt( kt, iiceprt, jiceprt, 2, ' - Final state - ' )                   ! prints 
    180183      IF( nn_timing == 1 )   CALL timing_stop ('icecor')                                                             ! timing 
     184write(numout,*) "RSRH leaving icecor", kt ; flush(numout) 
    181185      ! 
    182186   END SUBROUTINE ice_cor 
  • branches/UKMO/dev_r8126_LIM3_couple/NEMOGCM/NEMO/LIM_SRC_3/iceforcing.F90

    r8879 r8892  
    146146         ! 
    147147      CASE ( jp_purecpl )           !--- coupled formulation 
     148write(numout,*) "RSRH call cpl_ice_flx at_i_b",kt,glob_sum(at_i_b(:,:));flush(numout) 
     149write(numout,*) "RSRH call cpl_ice_flx alb_ice",kt,glob_sum(alb_ice(:,:));flush(numout) 
     150write(numout,*) "RSRH call cpl_ice_flx sst_m",kt,glob_sum(sst_m(:,:));flush(numout) 
     151 
    148152                                CALL sbc_cpl_ice_flx( picefr=at_i_b, palbi=alb_ice, psst=sst_m, pist=t_su, phs=h_s, phi=h_i ) 
     153write(numout,*) "RSRH call cpl_ice_flx",kt,glob_sum(at_i_b(:,:));flush(numout) 
    149154         IF( nn_iceflx == 2 )   CALL ice_flx_dist( t_su, alb_ice, qns_ice, qsr_ice, dqns_ice, evap_ice, devap_ice, nn_iceflx ) 
    150155         ! 
  • branches/UKMO/dev_r8126_LIM3_couple/NEMOGCM/NEMO/LIM_SRC_3/iceistate.F90

    r8879 r8892  
    365365         DO jl = 1, jpl 
    366366            a_ip_frac(:,:,jl) = rn_apnd * zswitch(:,:) * zfac 
     367WRITE(numout,*) "RSRH iceistate AA: a_ip_frac", glob_sum(a_ip_frac(:,:,:)) ; flush(numout) 
     368WRITE(numout,*) "RSRH iceistate AA: zswitch", glob_sum(zswitch(:,:)) ; flush(numout) 
     369WRITE(numout,*) "RSRH iceistate BB: h_ip pre calc", glob_sum(h_ip(:,:,:)) ; flush(numout) 
     370 
    367371            h_ip     (:,:,jl) = rn_hpnd * zswitch(:,:) * zfac 
     372WRITE(numout,*) "RSRH iceistate BB: h_ip post calc", glob_sum(h_ip(:,:,:)) ; flush(numout) 
    368373         END DO 
     374WRITE(numout,*) "RSRH iceistate CC: a_ip_frac pre a_ip calc", glob_sum(a_ip_frac(:,:,:)) ; flush(numout) 
     375WRITE(numout,*) "RSRH iceistate CC: a_i pre a_ip calc", glob_sum(a_i(:,:,:)) ; flush(numout) 
     376WRITE(numout,*) "RSRH iceistate CC: a_ip pre a_ip calc", glob_sum(a_ip(:,:,:)) ; flush(numout) 
    369377         a_ip(:,:,:) = a_ip_frac(:,:,:) * a_i (:,:,:)  
     378WRITE(numout,*) "RSRH iceistate DD: h_ip pre v_ip calc", glob_sum(h_ip(:,:,:)) ; flush(numout) 
     379WRITE(numout,*) "RSRH iceistate DD: a_ip pre v_ip calc", glob_sum(a_ip(:,:,:)) ; flush(numout) 
    370380         v_ip(:,:,:) = h_ip     (:,:,:) * a_ip(:,:,:) 
    371381 
  • branches/UKMO/dev_r8126_LIM3_couple/NEMOGCM/NEMO/LIM_SRC_3/iceitd.F90

    r8879 r8892  
    2626   USE in_out_manager ! I/O manager 
    2727   USE lib_mpp        ! MPP library 
    28    USE lib_fortran    ! fortran utilities (glob_sum + no signed zero) 
     28   USE lib_fortran    ! fortran utilities (sum + no signed zero) 
    2929   USE prtctl         ! Print control 
    3030 
     
    108108         zdhice(:,:) = 0._wp 
    109109         zhbnew(:,:) = 0._wp 
    110  
     110write(numout,*) "RSRH PE iceitd AA h_i", sum(h_i(:,:,:)); flush(numout) 
     111write(numout,*) "RSRH PE iceitd AA a_i", sum(a_i(:,:,:)); flush(numout) 
    111112         CALL tab_3d_2d( npti, nptidx(1:npti), h_i_2d (1:npti,1:jpl), h_i   ) 
    112113         CALL tab_3d_2d( npti, nptidx(1:npti), h_ib_2d(1:npti,1:jpl), h_i_b ) 
     
    197198         DO jl = 1, jpl 
    198199            ! 
     200write(numout,*) "RSRH PE iceitd BB h_i", sum(h_i(:,:,:)); flush(numout) 
     201write(numout,*) "RSRH PE iceitd BB a_i", sum(a_i(:,:,:)); flush(numout) 
    199202            CALL tab_2d_1d( npti, nptidx(1:npti), h_ib_1d(1:npti), h_i_b(:,:,jl) ) 
    200203            CALL tab_2d_1d( npti, nptidx(1:npti), h_i_1d (1:npti), h_i(:,:,jl)   ) 
     
    241244               END DO 
    242245               ! 
     246write(numout,*) "RSRH PE iceitd CC h_i", sum(h_i(:,:,:)); flush(numout) 
     247write(numout,*) "RSRH PE iceitd CC a_i", sum(a_i(:,:,:)); flush(numout) 
    243248               CALL tab_1d_2d( npti, nptidx(1:npti), h_i_1d (1:npti), h_i(:,:,jl)   ) 
    244249               CALL tab_1d_2d( npti, nptidx(1:npti), a_i_1d  (1:npti), a_i(:,:,jl)    ) 
     
    294299         ! 7) Make sure h_i >= minimum ice thickness hi_min 
    295300         !---------------------------------------------------------------------------------------------- 
     301write(numout,*) "RSRH PE iceitd DD h_i", sum(h_i(:,:,:)); flush(numout) 
     302write(numout,*) "RSRH PE iceitd DD a_i", sum(a_i(:,:,:)); flush(numout) 
     303write(numout,*) "RSRH PE iceitd DD a_ip", sum(a_ip(:,:,:)); flush(numout) 
    296304         CALL tab_2d_1d( npti, nptidx(1:npti), h_i_1d (1:npti), h_i(:,:,1)   ) 
    297305         CALL tab_2d_1d( npti, nptidx(1:npti), a_i_1d  (1:npti), a_i(:,:,1)    ) 
    298306         CALL tab_2d_1d( npti, nptidx(1:npti), a_ip_1d (1:npti), a_ip(:,:,1)   ) 
     307write(numout,*) "RSRH PE iceitd EE h_i", sum(h_i(:,:,:)); flush(numout) 
     308write(numout,*) "RSRH PE iceitd EE a_i", sum(a_i(:,:,:)); flush(numout) 
     309write(numout,*) "RSRH PE iceitd EE a_ip", sum(a_ip(:,:,:)); flush(numout) 
    299310          
    300311         DO ji = 1, npti 
     
    306317         END DO 
    307318         ! 
     319write(numout,*) "RSRH PE iceitd FF h_i", sum(h_i(:,:,:)); flush(numout) 
     320write(numout,*) "RSRH PE iceitd FF a_i", sum(a_i(:,:,:)); flush(numout) 
     321write(numout,*) "RSRH PE iceitd FF a_ip", sum(a_ip(:,:,:)); flush(numout) 
    308322         CALL tab_1d_2d( npti, nptidx(1:npti), h_i_1d (1:npti), h_i(:,:,1)   ) 
    309323         CALL tab_1d_2d( npti, nptidx(1:npti), a_i_1d  (1:npti), a_i(:,:,1)    ) 
    310324         CALL tab_1d_2d( npti, nptidx(1:npti), a_ip_1d (1:npti), a_ip(:,:,1)   ) 
     325write(numout,*) "RSRH PE iceitd GG h_i", sum(h_i(:,:,:)); flush(numout) 
     326write(numout,*) "RSRH PE iceitd GG a_i", sum(a_i(:,:,:)); flush(numout) 
     327write(numout,*) "RSRH PE iceitd GG a_ip", sum(a_ip(:,:,:)); flush(numout) 
    311328         ! 
    312329      ENDIF 
     
    527544      END WHERE 
    528545      ! 
     546write(numout,*) "RSRH PE iceitd HH h_i", sum(h_i(:,:,:)); flush(numout) 
     547write(numout,*) "RSRH PE iceitd HH a_i", sum(a_i(:,:,:)); flush(numout) 
     548write(numout,*) "RSRH PE iceitd HH a_ip", sum(a_ip(:,:,:)); flush(numout) 
    529549      CALL tab_2d_3d( npti, nptidx(1:npti), h_i_2d (1:npti,1:jpl), h_i  ) 
    530550      CALL tab_2d_3d( npti, nptidx(1:npti), a_i_2d (1:npti,1:jpl), a_i  ) 
     
    536556      CALL tab_2d_3d( npti, nptidx(1:npti), v_ip_2d(1:npti,1:jpl), v_ip ) 
    537557      CALL tab_2d_3d( npti, nptidx(1:npti), t_su_2d(1:npti,1:jpl), t_su ) 
     558write(numout,*) "RSRH PE iceitd II h_i", sum(h_i(:,:,:)); flush(numout) 
     559write(numout,*) "RSRH PE iceitd II a_i", sum(a_i(:,:,:)); flush(numout) 
     560write(numout,*) "RSRH PE iceitd II a_ip", sum(a_ip(:,:,:)); flush(numout) 
    538561      ! 
    539562   END SUBROUTINE ice_itd_shiftice 
     
    576599         END DO 
    577600         ! 
     601write(numout,*) "RSRH PE iceitd KK h_i",jl, sum(h_i(:,:,jl)); flush(numout) 
     602write(numout,*) "RSRH PE iceitd KK a_i",jl, sum(a_i(:,:,jl)); flush(numout) 
     603write(numout,*) "RSRH PE iceitd KK a_ip",jl, sum(a_ip(:,:,jl)); flush(numout) 
     604        ! 
    578605!!clem   CALL tab_2d_1d( npti, nptidx(1:npti), h_i_1d (1:npti), h_i(:,:,jl)   ) 
    579606         CALL tab_2d_1d( npti, nptidx(1:npti), a_i_1d  (1:npti), a_i(:,:,jl)    ) 
    580607         CALL tab_2d_1d( npti, nptidx(1:npti), v_i_1d  (1:npti), v_i(:,:,jl)    ) 
     608write(numout,*) "RSRH PE iceitd LL h_i",jl, sum(h_i(:,:,jl)); flush(numout) 
     609write(numout,*) "RSRH PE iceitd LL a_i",jl, sum(a_i(:,:,jl)); flush(numout) 
     610write(numout,*) "RSRH PE iceitd LL a_ip",jl, sum(a_ip(:,:jl)); flush(numout) 
    581611         ! 
    582612         DO ji = 1, npti 
     
    606636      !                       !----------------------------------------- 
    607637      DO jl = jpl-1, 1, -1    ! Identify thicknesses that are too small 
     638write(numout,*) "RSRH PE iceitd LL for jl=",jl ;flush(numout) 
    608639         !                    !----------------------------------------- 
    609640         npti = 0 ; nptidx(:) = 0 
     
    616647            END DO 
    617648         END DO 
    618          ! 
     649write(numout,*) "RSRH PE iceitd LL after 1st loop for jl=",jl ;flush(numout) 
     650write(numout,*) "RSRH PE iceitd get summing LL after 1st loop for jl=",jl ;flush(numout) 
     651         ! 
     652write(numout,*) "RSRH PE iceitd MM pre h_i",jl, sum(h_i(:,:,jl)); flush(numout) 
     653write(numout,*) "RSRH PE iceitd MM pre a_i",jl, sum(a_i(:,:,jl)); flush(numout) 
     654write(numout,*) "RSRH PE iceitd MM pre a_ip",jl, sum(a_ip(:,:jl)); flush(numout) 
    619655         CALL tab_2d_1d( npti, nptidx(1:npti), a_i_1d  (1:npti), a_i(:,:,jl+1)    ) ! jl+1 is ok 
    620656         CALL tab_2d_1d( npti, nptidx(1:npti), v_i_1d  (1:npti), v_i(:,:,jl+1)    ) ! jl+1 is ok 
     657write(numout,*) "RSRH PE iceitd MM h_i",jl, sum(h_i(:,:,jl)); flush(numout) 
     658write(numout,*) "RSRH PE iceitd MM a_i",jl, sum(a_i(:,:,jl)); flush(numout) 
     659write(numout,*) "RSRH PE iceitd MM a_ip",jl, sum(a_ip(:,:jl)); flush(numout) 
    621660         DO ji = 1, npti 
    622661            jdonor(ji,jl) = jl + 1 
     
    624663            zdvice(ji,jl) = v_i_1d(ji) 
    625664         END DO 
     665write(numout,*) "RSRH PE iceitd MM after 2nd loop for jl=",jl ;flush(numout) 
    626666         ! 
    627667         IF( npti > 0 ) THEN 
     
    632672            zdvice(1:npti,jl) = 0._wp 
    633673         ENDIF 
     674write(numout,*) "RSRH PE iceitd NN after 3rd loop for jl=",jl ;flush(numout) 
    634675         ! 
    635676      END DO 
     677write(numout,*) "RSRH PE iceitd NN out of NN loop for jl=",jl ;flush(numout) 
    636678      ! 
    637679   END SUBROUTINE ice_itd_reb 
  • branches/UKMO/dev_r8126_LIM3_couple/NEMOGCM/NEMO/LIM_SRC_3/icestp.F90

    r8879 r8892  
    107107      IF( MOD( kt-1, nn_fsbc ) == 0 ) THEN   ! --- Ice time step --- ! 
    108108         !                                   !-----------------------! 
    109          ! 
     109         !meaning 
    110110         kt_ice = kt                              ! -- Ice model time step 
    111111         ! 
     
    116116         v_oce(:,:) = ssv_m(:,:) 
    117117         ! 
     118write(numout,*) "RSRH ice_stp 00 before eos_stp ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout) 
     119write(numout,*) "RSRH ice_stp 00 qns,sfx,emp before eos_stp ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout) 
    118120         CALL eos_fzp( sss_m(:,:) , t_bo(:,:) )   ! -- freezing temperature [Kelvin] (set to rt0 over land) 
    119121         t_bo(:,:) = ( t_bo(:,:) + rt0 ) * tmask(:,:,1) + rt0 * ( 1._wp - tmask(:,:,1) ) 
     
    134136         CALL ice_rst_opn( kt )     ! Open Ice restart file (if necessary)  
    135137         ! 
     138write(numout,*) "RSRH ice_stp 00 before ice_dyn ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout) 
     139write(numout,*) "RSRH ice_stp 00 qns,sfx,emp before ice_dyn ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout) 
    136140         IF( ln_icedyn .AND. .NOT.lk_c1d )   & 
    137141            &                           CALL ice_dyn( kt )            ! -- Ice dynamics 
     142write(numout,*) "RSRH ice_stp 00 after ice_dyn ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout) 
     143write(numout,*) "RSRH ice_stp 00 qns,sfx,emp aft ice_dyn ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout) 
    138144         ! 
    139145         !                          !==  lateral boundary conditions  ==! 
     
    145151         ! 
    146152         !                          !==  previous lead fraction and ice volume for flux calculations 
     153write(numout,*) "RSRH ice_stp 00 be4 glo2eqv ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout) 
     154write(numout,*) "RSRH ice_stp 00 qns,sfx,emp b4 glo2eqv ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout) 
    147155         CALL ice_var_glo2eqv            ! h_i and h_s for ice albedo calculation 
     156write(numout,*) "RSRH ice_stp 00 be4 ice_var_agg ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout) 
    148157         CALL ice_var_agg(1)             ! at_i for coupling  
     158write(numout,*) "RSRH ice_stp 00 qns,sfx,emp b4 store_filds ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout) 
     159write(numout,*) "RSRH ice_stp 00 be4 store_filds ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout) 
    149160         CALL store_fields               ! Store now ice values 
    150161 
     
    162173         !    qprec_ice, qevap_ice 
    163174         !------------------------------------------------------! 
     175write(numout,*) "RSRH ice_stp 00 qns,sfx,emp before forc ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout) 
     176write(numout,*) "RSRH ice_stp 00 be4 forc ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout) 
    164177                                        CALL ice_forcing_flx( kt, ksbc ) 
    165178         !----------------------------! 
    166179         ! --- ice thermodynamics --- ! 
    167180         !----------------------------! 
     181write(numout,*) "RSRH ice_stp 00 qns,sfx,emp before thd ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout) 
     182write(numout,*) "RSRH ice_stp 00 be4 thd", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout) 
    168183         IF( ln_icethd )                CALL ice_thd( kt )            ! -- Ice thermodynamics       
    169184         ! 
     185write(numout,*) "RSRH ice_stp 00 qns,sfx,emp before icecor ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout) 
     186write(numout,*) "RSRH ice_stp 00 before icecor ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout) 
    170187         ! 
    171188         IF( ln_icethd )                CALL ice_cor( kt , 2 )        ! -- Corrections 
     189write(numout,*) "RSRH ice_stp AA qns,sfx,emp  ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout) 
     190write(numout,*) "RSRH ice_stp AA ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout) 
    172191# if defined key_agrif 
    173192         IF( .NOT. Agrif_Root() )       CALL agrif_update_lim3( kt ) 
    174193# endif 
    175194         CALL ice_var_glo2eqv        ! necessary calls (at least for coupling) 
     195write(numout,*) "RSRH ice_stp BB qns,sfx,emp  ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout) 
     196write(numout,*) "RSRH ice_stp BB ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout) 
    176197         CALL ice_var_agg( 2 )       ! necessary calls (at least for coupling) 
     198write(numout,*) "RSRH ice_stp CC qns,sfx,emp  ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout) 
     199write(numout,*) "RSRH ice_stp CC ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout) 
    177200                                     ! 
    178201!! clem: one should switch the calculation of the fluxes onto the parent grid but the following calls do not work 
     
    182205!!# endif 
    183206                                        CALL ice_update_flx( kt )     ! -- Update ocean surface mass, heat and salt fluxes 
     207write(numout,*) "RSRH ice_stp DD qns,sfx,emp  ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout) 
     208write(numout,*) "RSRH ice_stp DD ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout) 
    184209!!# if defined key_agrif 
    185210!!         IF( .NOT. Agrif_Root() )     CALL Agrif_ParentGrid_To_ChildGrid() 
    186211!!# endif 
    187212         IF( ln_icediahsb )             CALL ice_dia( kt )            ! -- Diagnostics outputs  
     213write(numout,*) "RSRH ice_stp EE qns,sfx,emp  ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout) 
     214write(numout,*) "RSRH ice_stp EE ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout) 
    188215         ! 
    189216                                        CALL ice_wri( kt )            ! -- Ice outputs  
     217write(numout,*) "RSRH ice_stp FF qns,sfx,emp  ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout) 
     218write(numout,*) "RSRH ice_stp FF ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout) 
    190219         ! 
    191220         IF( lrst_ice )                 CALL ice_rst_write( kt )      ! -- Ice restart file  
     221write(numout,*) "RSRH ice_stp HH qns,sfx,emp  ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout) 
     222write(numout,*) "RSRH ice_stp HH ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout) 
    192223         ! 
    193224         IF( ln_icectl )                CALL ice_ctl( kt )            ! -- alerts in case of model crash 
     225write(numout,*) "RSRH ice_stp II qns,sfx,emp  ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout) 
     226write(numout,*) "RSRH ice_stp II ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout) 
    194227         ! 
    195228      ENDIF   ! End sea-ice time step only 
     
    199232      !-------------------------! 
    200233      IF( ln_icedyn )                   CALL ice_update_tau( kt, ub(:,:,1), vb(:,:,1) )   ! -- update surface ocean stresses 
     234write(numout,*) "RSRH ice_stp JJ qns,sfx,emp  ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout) 
     235write(numout,*) "RSRH ice_stp JJ ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout) 
    201236!!gm   remark, the ocean-ice stress is not saved in ice diag call above .....  find a solution!!! 
    202237      ! 
     
    360395       
    361396      ! ice velocities & total concentration 
     397 
     398 
     399write(numout,*) "RSRH store_fields a_i_b sum", glob_sum(a_i_b(:,:,:));flush(numout) 
     400 
    362401      at_i_b(:,:)  = SUM( a_i_b(:,:,:), dim=3 ) 
    363402      u_ice_b(:,:) = u_ice(:,:) 
  • branches/UKMO/dev_r8126_LIM3_couple/NEMOGCM/NEMO/LIM_SRC_3/icethd.F90

    r8879 r8892  
    9999      ! controls 
    100100      IF( nn_timing == 1 )   CALL timing_start('icethd')                                                             ! timing 
     101 
     102write(numout,*) "RSRH icethd AA", kt; flush(numout) 
    101103      IF( ln_icediachk   )   CALL ice_cons_hsm(0, 'icethd', rdiag_v, rdiag_s, rdiag_t, rdiag_fv, rdiag_fs, rdiag_ft) ! conservation 
    102  
     104write(numout,*) "RSRH icethd BB", kt; flush(numout) 
     105 
     106write(numout,*) "RSRH icethd BB u_ice", kt, glob_sum(u_ice(:,:)); flush(numout) 
     107write(numout,*) "RSRH icethd BB v_ice", kt, glob_sum(v_ice(:,:)); flush(numout) 
     108write(numout,*) "RSRH icethd BB ssu_m", kt, glob_sum(ssu_m(:,:)); flush(numout) 
     109write(numout,*) "RSRH icethd BB ssv_m", kt, glob_sum(ssv_m(:,:)); flush(numout) 
    103110      IF( kt == nit000 .AND. lwp ) THEN 
    104111         WRITE(numout,*) 
     
    106113         WRITE(numout,*) '~~~~~~~' 
    107114      ENDIF 
     115write(numout,*) "RSRH icethd CC", kt; flush(numout) 
    108116       
    109117      CALL ice_var_glo2eqv 
     118write(numout,*) "RSRH icethd DD", kt; flush(numout) 
     119       
     120      ! RSRH. zfric needs to be properly initialised....Duuuuhh! 
     121      zfric(:,:) = 0.0  
    110122 
    111123      !---------------------------------------------! 
     
    113125      !---------------------------------------------! 
    114126      IF( ln_icedyn ) THEN 
     127write(numout,*) "RSRH icethd EE in", kt; flush(numout) 
    115128         zu_io(:,:) = u_ice(:,:) - ssu_m(:,:) 
    116129         zv_io(:,:) = v_ice(:,:) - ssv_m(:,:) 
     130 
     131write(numout,*) "RSRH icethd EE u_ice", kt, glob_sum(u_ice(:,:)); flush(numout) 
     132write(numout,*) "RSRH icethd EE v_ice", kt, glob_sum(v_ice(:,:)); flush(numout) 
     133write(numout,*) "RSRH icethd EE ssu_m", kt, glob_sum(ssu_m(:,:)); flush(numout) 
     134write(numout,*) "RSRH icethd EE ssv_m", kt, glob_sum(ssv_m(:,:)); flush(numout) 
     135write(numout,*) "RSRH icethd EE zu_io", kt, glob_sum(zu_io(:,:)); flush(numout) 
     136write(numout,*) "RSRH icethd EE zv_io", kt, glob_sum(zv_io(:,:)); flush(numout) 
     137 
     138write(numout,*) "RSRH PE icethd EE u_ice", kt, sum(u_ice(:,:)); flush(numout) 
     139write(numout,*) "RSRH PE icethd EE v_ice", kt, sum(v_ice(:,:)); flush(numout) 
     140write(numout,*) "RSRH PE icethd EE ssu_m", kt, sum(ssu_m(:,:)); flush(numout) 
     141write(numout,*) "RSRH PE icethd EE ssv_m", kt, sum(ssv_m(:,:)); flush(numout) 
     142write(numout,*) "RSRH PE icethd EE zu_io", kt, sum(zu_io(:,:)); flush(numout) 
     143write(numout,*) "RSRH PE icethd EE zv_io", kt, sum(zv_io(:,:)); flush(numout) 
     144     
     145 
    117146         DO jj = 2, jpjm1  
    118147            DO ji = fs_2, fs_jpim1 
     
    122151            END DO 
    123152         END DO 
     153write(numout,*) "RSRH icethd EE", kt; flush(numout) 
    124154      ELSE      !  if no ice dynamics => transmit directly the atmospheric stress to the ocean 
     155write(numout,*) "RSRH icethd FF in", kt; flush(numout) 
    125156         DO jj = 2, jpjm1 
    126157            DO ji = fs_2, fs_jpim1 
     
    131162         END DO 
    132163      ENDIF 
     164write(numout,*) "RSRH icethd GG", kt; flush(numout) 
     165write(numout,*) "RSRH icethd GG pre zfric link", kt; flush(numout) 
     166write(numout,*) "RSRH icethd GG sum pre zfric link", kt, glob_sum(zfric(:,:)); flush(numout) 
     167write(numout,*) "RSRH PE icethd GG sum pre zfric link SUM", kt, sum(zfric(:,:)); flush(numout) 
     168 
    133169      CALL lbc_lnk( zfric, 'T',  1. ) 
    134       ! 
     170 write(numout,*) "RSRH icethd GG after zfric link", kt; flush(numout) 
     171     ! 
    135172      ftr_ice(:,:,:) = 0._wp  ! initialization (part of solar radiation transmitted through the ice) 
     173 
     174write(numout,*) "RSRH GG at_i", glob_sum(at_i(:,:)); flush(numout) 
     175write(numout,*) "RSRH GG at_i_b", glob_sum(at_i_b(:,:)); flush(numout) 
     176write(numout,*) "RSRH GG qsr_oce", glob_sum(qsr_oce(:,:)); flush(numout) 
     177write(numout,*) "RSRH GG qns_oce", glob_sum(qns_oce(:,:)); flush(numout) 
     178write(numout,*) "RSRH GG frq_m", glob_sum(frq_m(:,:)); flush(numout) 
     179write(numout,*) "RSRH GG qemp_oce", glob_sum(qemp_oce(:,:)); flush(numout) 
     180write(numout,*) "RSRH GG e3t_m", glob_sum(e3t_m(:,:)); flush(numout) 
     181write(numout,*) "RSRH GG t_bo", glob_sum(t_bo(:,:)); flush(numout) 
     182write(numout,*) "RSRH GG sst_m", glob_sum(sst_m(:,:)); flush(numout) 
     183 
     184write(numout,*) "RSRH GG zfric", glob_sum(zfric(:,:)); flush(numout) 
    136185 
    137186      !--------------------------------------------------------------------! 
     
    185234         END DO 
    186235      END DO 
     236write(numout,*) "RSRH icethd HH", kt; flush(numout) 
    187237       
    188238      ! In case we bypass open-water ice formation 
     
    206256         &           - at_i (:,:) *  fhld(:,:)                                  ! heat flux taken during bottom growth/melt  
    207257                                                                                !    (fhld should be 0 while bott growth) 
     258write(numout,*) "RSRH icethd II", kt; flush(numout) 
    208259      !-------------------------------------------------------------------------------------------! 
    209260      ! Thermodynamic computation (only on grid points covered by ice) => loop over ice categories 
     
    211262      DO jl = 1, jpl 
    212263 
    213          ! select ice covered grid points 
     264 write(numout,*) "RSRH icethd JJ for Jl",jl, kt; flush(numout) 
     265        ! select ice covered grid points 
    214266         npti = 0 ; nptidx(:) = 0 
    215267         DO jj = 1, jpj 
     
    222274         END DO 
    223275 
     276write(numout,*) "RSRH icethd KK for Jl",jl, kt; flush(numout) 
    224277         IF( lk_mpp )         CALL mpp_ini_ice( npti , numout ) 
    225278 
     279write(numout,*) "RSRH icethd LL for Jl",jl, kt; flush(numout) 
    226280         IF( npti > 0 ) THEN  ! If there is no ice, do nothing. 
    227281            !                                                                 
    228282                              CALL ice_thd_1d2d( jl, 1 )            ! --- Move to 1D arrays --- ! 
    229283            !                                                       ! --- & Change units of e_i, e_s from J/m2 to J/m3 --- ! 
     284write(numout,*) "RSRH icethd MM for Jl",jl, kt; flush(numout) 
    230285            ! 
    231286            s_i_new   (1:npti) = 0._wp ; dh_s_tot (1:npti) = 0._wp  ! --- some init --- !  (important to have them here)  
     
    239294                              CALL ice_thd_ent( e_i_1d(1:npti,:) )         ! Ice enthalpy remapping 
    240295            ENDIF 
     296write(numout,*) "RSRH icethd NN for Jl",jl, kt; flush(numout) 
    241297            ! 
    242298                              CALL ice_thd_sal( ln_icedS )          ! --- Ice salinity --- !     
     299write(numout,*) "RSRH icethd OO for Jl",jl, kt; flush(numout) 
    243300            ! 
    244301                              CALL ice_thd_temp                     ! --- temperature update --- ! 
     302write(numout,*) "RSRH icethd PP for Jl",jl, kt; flush(numout) 
    245303            ! 
    246304!!gm please create a new logical (l_thd_lam or a better explicit name) set one for all in icestp.F90 module 
     
    254312            END IF 
    255313            ! 
     314write(numout,*) "RSRH icethd QQ for Jl",jl, kt; flush(numout) 
    256315            IF( ln_icedA )    CALL ice_thd_da                       ! --- lateral melting --- ! 
    257316            ! 
     317write(numout,*) "RSRH icethd RR for Jl",jl, kt; flush(numout) 
    258318                              CALL ice_thd_1d2d( jl, 2 )            ! --- Change units of e_i, e_s from J/m3 to J/m2 --- ! 
    259319            !                                                       ! --- & Move to 2D arrays --- ! 
     320write(numout,*) "RSRH icethd SS for Jl",jl, kt; flush(numout) 
    260321            ! 
    261322            IF( lk_mpp )      CALL mpp_comm_free( ncomm_ice ) !RB necessary ?? 
     
    267328 
    268329      IF( ln_icediachk )   CALL ice_cons_hsm(1, 'icethd', rdiag_v, rdiag_s, rdiag_t, rdiag_fv, rdiag_fs, rdiag_ft) 
     330write(numout,*) "RSRH icethd UU for Jl",jl, kt; flush(numout) 
    269331      ! 
    270332                           CALL ice_var_zapsmall           ! --- remove very small ice concentration (<1e-10) --- ! 
     333write(numout,*) "RSRH icethd VV for Jl",jl, kt; flush(numout) 
    271334      !                                                    !     & make sure at_i=SUM(a_i) & ato_i=1 where at_i=0 
    272335      !                    
  • branches/UKMO/dev_r8126_LIM3_couple/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90

    r8879 r8892  
    4444   ! 
    4545   USE in_out_manager ! I/O manager 
     46   use lib_fortran 
    4647   USE iom            ! NetCDF library 
    4748   USE lib_mpp        ! distribued memory computing library 
     
    11331134      ! 
    11341135      IF( ln_mixcpl )   zmsk(:,:) = 1. - xcplmask(:,:,0) 
     1136 
     1137write(numout,*) "RSRH in sbc_cpl_rcv",kt 
     1138write(numout,*) "RSRH in sbc_cpl_rcv qns, emp, sfx",kt, glob_sum(qns(:,:)), glob_sum(emp(:,:)),glob_sum( sfx(:,:));flush(numout) 
     1139 
    11351140      ! 
    11361141      !                                                      ! ======================================================= ! 
     
    12221227      ENDIF 
    12231228      !                                                      ! ========================= ! 
     1229write(numout,*) "RSRH in sbc_cpl_rcv BB qns, emp, sfx",kt, glob_sum(qns(:,:)), glob_sum(emp(:,:)),glob_sum( sfx(:,:));flush(numout) 
    12241230      !                                                      !    wind stress module     !   (taum) 
    12251231      !                                                      ! ========================= ! 
     
    12841290 
    12851291      !                                                      ! ================== ! 
     1292write(numout,*) "RSRH in sbc_cpl_rcv CC qns, emp, sfx",kt, glob_sum(qns(:,:)), glob_sum(emp(:,:)),glob_sum( sfx(:,:));flush(numout) 
    12861293      !                                                      ! atmosph. CO2 (ppm) ! 
    12871294      !                                                      ! ================== ! 
     
    14211428      ENDIF 
    14221429       
     1430write(numout,*) "RSRH in sbc_cpl_rcv QQ qns, emp, sfx",kt, glob_sum(qns(:,:)), glob_sum(emp(:,:)),glob_sum( sfx(:,:));flush(numout) 
    14231431      !                                                      ! ========================= ! 
    14241432      IF( k_ice <= 1 .AND. MOD( kt-1, k_fsbc ) == 0 ) THEN   !  heat & freshwater fluxes ! (Ocean only case) 
     
    14521460         ELSE                   ;   emp(:,:) =                              zemp(:,:) 
    14531461         ENDIF 
     1462write(numout,*) "RSRH in sbc_cpl_rcv VV qns, emp, sfx",kt, glob_sum(qns(:,:)), glob_sum(emp(:,:)),glob_sum( sfx(:,:));flush(numout) 
    14541463         ! 
    14551464         !                                                       ! non solar heat flux over the ocean (qns) 
     
    14881497         ! 
    14891498      ENDIF 
     1499write(numout,*) "RSRH in sbc_cpl_rcv ZZ qns, emp, sfx",kt, glob_sum(qns(:,:)), glob_sum(emp(:,:)),glob_sum( sfx(:,:));flush(numout) 
    14901500      ! 
    14911501      CALL wrk_dealloc( jpi,jpj, ztx, zty, zmsk, zemp, zqns, zqsr, ztx2, zty2 ) 
     
    17781788      CALL wrk_alloc( jpi,jpj,jpl, zqns_ice, zqsr_ice, zdqns_ice, zqevap_ice, zfrqsr_tr_i ) 
    17791789 
     1790write(numout,*) "RSRH in sbc_cpl_ice_flx AA qns, emp, sfx", glob_sum(qns(:,:)), glob_sum(emp(:,:)),glob_sum( sfx(:,:));flush(numout) 
     1791write(numout,*) "RSRH in sbc_cpl_ice_flx AA emp_tot, emp_ice, emp_oce", glob_sum(emp_tot(:,:)), glob_sum(emp_ice(:,:)), glob_sum(emp_oce(:,:));flush(numout) 
     1792 
     1793 
    17801794      IF (sn_rcv_emp%clcat == 'yes') THEN 
    17811795         CALL wrk_alloc( jpi, jpj, jpl, zevap_ice) 
     
    17981812      SELECT CASE( TRIM( sn_rcv_emp%cldes ) ) 
    17991813      CASE( 'conservative' )   ! received fields: jpr_rain, jpr_snow, jpr_ievp, jpr_tevp 
     1814write(numout,*) "RSRH in sbc_cpl_ice_flx BB con case" 
    18001815         zsprecip(:,:) =   frcv(jpr_snow)%z3(:,:,1)                  ! May need to ensure positive here 
    18011816         ztprecip(:,:) =   frcv(jpr_rain)%z3(:,:,1) + zsprecip(:,:)  ! May need to ensure positive here 
    18021817         zemp_tot(:,:) =   frcv(jpr_tevp)%z3(:,:,1) - ztprecip(:,:) 
    18031818         zemp_ice(:,:) = ( frcv(jpr_ievp)%z3(:,:,1) - frcv(jpr_snow)%z3(:,:,1) ) * picefr(:,:) 
     1819write(numout,*) "RSRH in sbc_cpl_ice_flx BB con casejpr_snow ",glob_sum(frcv(jpr_snow)%z3(:,:,1)); flush(numout) 
     1820write(numout,*) "RSRH in sbc_cpl_ice_flx BB con casejpr_rain ",glob_sum(frcv(jpr_rain)%z3(:,:,1)); flush(numout) 
     1821write(numout,*) "RSRH in sbc_cpl_ice_flx BB con casejpr_tevp ",glob_sum(frcv(jpr_tevp)%z3(:,:,1)); flush(numout) 
     1822write(numout,*) "RSRH in sbc_cpl_ice_flx BB con casejpr_ievp ",glob_sum(frcv(jpr_ievp)%z3(:,:,1)); flush(numout) 
     1823write(numout,*) "RSRH in sbc_cpl_ice_flx BB con case picefr ",glob_sum(picefr(:,:)); flush(numout) 
    18041824      CASE( 'oce and ice'   )   ! received fields: jpr_sbpr, jpr_semp, jpr_oemp, jpr_ievp 
     1825write(numout,*) "RSRH in sbc_cpl_ice_flx BB oce and ice case" 
    18051826         zemp_tot(:,:) = ziceld(:,:) * frcv(jpr_oemp)%z3(:,:,1) + picefr(:,:) * frcv(jpr_sbpr)%z3(:,:,1) 
    18061827         zemp_ice(:,:) = frcv(jpr_semp)%z3(:,:,1) * picefr(:,:) 
     
    18081829         ztprecip(:,:) = frcv(jpr_semp)%z3(:,:,1) - frcv(jpr_sbpr)%z3(:,:,1) + zsprecip(:,:) 
    18091830      END SELECT 
     1831write(numout,*) "RSRH in sbc_cpl_ice_flx BB what case?",TRIM( sn_rcv_emp%cldes ) 
     1832 
     1833write(numout,*) "RSRH in sbc_cpl_ice_flx BB emp_tot, emp_ice, emp_oce", glob_sum(emp_tot(:,:)), glob_sum(emp_ice(:,:)), glob_sum(emp_oce(:,:));flush(numout) 
     1834write(numout,*) "RSRH in sbc_cpl_ice_flx BB zemp_tot, zemp_ice", glob_sum(zemp_tot(:,:)), glob_sum(zemp_ice(:,:));flush(numout) 
    18101835 
    18111836#if defined key_lim3 
     
    19061931 
    19071932#endif 
     1933write(numout,*) "RSRH in sbc_cpl_ice_flx LL emp_tot, emp_ice, emp_oce", glob_sum(emp_tot(:,:)), glob_sum(emp_ice(:,:)), glob_sum(emp_oce(:,:));flush(numout) 
    19081934      ! outputs 
    19091935!!      IF( srcv(jpr_rnf)%laction )   CALL iom_put( 'runoffs' , rnf(:,:) * tmask(:,:,1)                                 )  ! runoff 
     
    19752001      ! Heat content per unit mass of rain (J/kg) 
    19762002      zcptrain(:,:) = rcp * ( SUM( (tn_ice(:,:,:) - rt0) * a_i(:,:,:), dim=3 ) + sst_m(:,:) * ziceld(:,:) )  
     2003write(numout,*) "RSRH in sbc_cpl_ice_flx NN zcptsnw,zcptn", glob_sum(zcptsnw(:,:)), glob_sum(zcptn(:,:));flush(numout) 
    19772004 
    19782005      ! --- enthalpy of snow precip over ice in J/m3 (to be used in 1D-thermo) --- ! 
     
    19912018!!    zqemp_ice(:,:) = -   frcv(jpr_ievp)%z3(:,:,1)        * picefr(:,:)      *   zcptsnw (:,:)   &        ! ice evap 
    19922019!!       &             +   zsprecip(:,:)                   * zsnw             * zqprec_ice(:,:) * r1_rhosn ! solid precip over ice 
     2020write(numout,*) "RSRH in sbc_cpl_ice_flx NN zsprecip,zcptrain", glob_sum(zsprecip(:,:)), glob_sum(zcptrain(:,:));flush(numout) 
    19932021       
    19942022      ! --- total non solar flux (including evap/precip) --- ! 
     
    19972025      ! --- in case both coupled/forced are active, we must mix values --- !  
    19982026      IF( ln_mixcpl ) THEN 
     2027write(numout,*) "RSRH in sbc_cpl_ice_flx OO.1 zqns_tot,zqns_ice, zqns_oce", glob_sum(zqns_tot(:,:)), glob_sum(zqns_ice(:,:,:)), glob_sum(zqns_oce(:,:));flush(numout) 
     2028write(numout,*) "RSRH in sbc_cpl_ice_flx OO.1 xcplmask", glob_sum(xcplmask(:,:,0));flush(numout) 
     2029write(numout,*) "RSRH in sbc_cpl_ice_flx OO.1 xcplmask", glob_sum(Zmsk(:,:));flush(numout) 
    19992030         qns_tot(:,:) = qns_tot(:,:) * xcplmask(:,:,0) + zqns_tot(:,:)* zmsk(:,:) 
    20002031         qns_oce(:,:) = qns_oce(:,:) * xcplmask(:,:,0) + zqns_oce(:,:)* zmsk(:,:) 
     2032write(numout,*) "RSRH in sbc_cpl_ice_flx OO.2 zqns_tot,zqns_ice, zqns_oce", glob_sum(zqns_tot(:,:)), glob_sum(zqns_ice(:,:,:)), glob_sum(zqns_oce(:,:));flush(numout) 
    20012033         DO jl=1,jpl 
    20022034            qns_ice  (:,:,jl) = qns_ice  (:,:,jl) * xcplmask(:,:,0) +  zqns_ice  (:,:,jl)* zmsk(:,:) 
    20032035            qevap_ice(:,:,jl) = qevap_ice(:,:,jl) * xcplmask(:,:,0) +  zqevap_ice(:,:,jl)* zmsk(:,:) 
    20042036         ENDDO 
     2037write(numout,*) "RSRH in sbc_cpl_ice_flx OO.3 zqns_tot,zqns_ice, zqns_oce", glob_sum(zqns_tot(:,:)), glob_sum(zqns_ice(:,:,:)), glob_sum(zqns_oce(:,:));flush(numout) 
    20052038         qprec_ice(:,:) = qprec_ice(:,:) * xcplmask(:,:,0) + zqprec_ice(:,:)* zmsk(:,:) 
    20062039         qemp_oce (:,:) =  qemp_oce(:,:) * xcplmask(:,:,0) +  zqemp_oce(:,:)* zmsk(:,:) 
    20072040         qemp_ice (:,:) =  qemp_ice(:,:) * xcplmask(:,:,0) +  zqemp_ice(:,:)* zmsk(:,:) 
     2041write(numout,*) "RSRH in sbc_cpl_ice_flx OO.4 zqns_tot,zqns_ice, zqns_oce", glob_sum(zqns_tot(:,:)), glob_sum(zqns_ice(:,:,:)), glob_sum(zqns_oce(:,:));flush(numout) 
    20082042      ELSE 
     2043write(numout,*) "RSRH in sbc_cpl_ice_flx PP zqns_tot,zqns_ice, zqns_oce", glob_sum(zqns_tot(:,:)), glob_sum(zqns_ice(:,:,:)), glob_sum(zqns_oce(:,:));flush(numout) 
    20092044         qns_tot  (:,:  ) = zqns_tot  (:,:  ) 
    20102045         qns_oce  (:,:  ) = zqns_oce  (:,:  ) 
     
    20152050         qemp_ice (:,:  ) = zqemp_ice (:,:  ) 
    20162051      ENDIF 
     2052write(numout,*) "RSRH in sbc_cpl_ice_flx PP qns_tot, qns_ice, qns_oce", glob_sum(qns_tot(:,:)), glob_sum(qns_ice(:,:)), glob_sum(qns_oce(:,:));flush(numout) 
    20172053 
    20182054#else 
     
    20492085      ! note: hflx for runoff and iceshelf are done in sbcrnf and sbcisf resp. 
    20502086      ! 
     2087write(numout,*) "RSRH in sbc_cpl_ice_flx MM qns, emp, sfx", glob_sum(qns(:,:)), glob_sum(emp(:,:)),glob_sum( sfx(:,:));flush(numout) 
    20512088      !                                                      ! ========================= ! 
    20522089      SELECT CASE( TRIM( sn_rcv_qsr%cldes ) )                !      solar heat fluxes    !   (qsr) 
     
    22022239         CALL wrk_dealloc( jpi,jpj,1,zevap_ice) 
    22032240      ENDIF 
     2241write(numout,*) "RSRH in sbc_cpl_ice_flx ZZ qns, emp, sfx", glob_sum(qns(:,:)), glob_sum(emp(:,:)),glob_sum( sfx(:,:));flush(numout) 
    22042242      ! 
    22052243      IF( nn_timing == 1 )   CALL timing_stop('sbc_cpl_ice_flx') 
  • branches/UKMO/dev_r8126_LIM3_couple/NEMOGCM/NEMO/OPA_SRC/SBC/sbcssm.F90

    r8879 r8892  
    2020   ! 
    2121   USE in_out_manager ! I/O manager 
     22   use lib_fortran 
    2223   USE prtctl         ! Print control 
    2324   USE iom            ! IOM library 
     
    5960      ! 
    6061      !                                        !* surface T-, U-, V- ocean level variables (T, S, depth, velocity) 
     62write(numout,*) "RSRH sbcssm tsb", kt, glob_sum(tsb(:,:,:,jp_tem)); flush(numout) 
     63write(numout,*) "RSRH sbcssm tsn", kt, glob_sum(tsn(:,:,:,jp_tem)); flush(numout) 
     64write(numout,*) "RSRH sbcssm tsa", kt, glob_sum(tsa(:,:,:,jp_tem)); flush(numout) 
     65write(numout,*) "RSRH PE sbcssm tsb", kt, sum(tsb(:,:,:,jp_tem)); flush(numout) 
     66write(numout,*) "RSRH PE sbcssm tsn", kt, sum(tsn(:,:,:,jp_tem)); flush(numout) 
     67write(numout,*) "RSRH PE sbcssm tsa", kt, sum(tsa(:,:,:,jp_tem)); flush(numout) 
    6168      DO jj = 1, jpj 
    6269         DO ji = 1, jpi 
    63             zts(ji,jj,jp_tem) = tsn(ji,jj,mikt(ji,jj),jp_tem) 
    64             zts(ji,jj,jp_sal) = tsn(ji,jj,mikt(ji,jj),jp_sal) 
     70            zts(ji,jj,jp_tem) = tsb(ji,jj,mikt(ji,jj),jp_tem) 
     71            zts(ji,jj,jp_sal) = tsb(ji,jj,mikt(ji,jj),jp_sal) 
    6572         END DO 
    6673      END DO 
     
    6875      IF( nn_fsbc == 1 ) THEN                             !   Instantaneous surface fields        ! 
    6976         !                                                ! ---------------------------------------- ! 
     77write(numout,*) "RSRH sbcssm nnfsbc=1", kt; flush(numout) 
     78write(numout,*) "RSRH sbcssm ub 00", kt, glob_sum(ub(:,:,1)); flush(numout) 
     79write(numout,*) "RSRH sbcssm vb 00", kt, glob_sum(vb(:,:,1)); flush(numout) 
     80write(numout,*) "RSRH sbcssm ssu_m 00", kt, glob_sum(ssu_m(:,:)); flush(numout) 
     81write(numout,*) "RSRH sbcssm ssv_m 00", kt, glob_sum(ssv_m(:,:)); flush(numout) 
     82write(numout,*) "RSRH sbcssm sst_m 00", kt, glob_sum(sst_m(:,:)); flush(numout) 
     83 
    7084         ssu_m(:,:) = ub(:,:,1) 
    7185         ssv_m(:,:) = vb(:,:,1) 
    72          IF( l_useCT )  THEN    ;   sst_m(:,:) = eos_pt_from_ct( zts(:,:,jp_tem), zts(:,:,jp_sal) ) 
    73          ELSE                    ;   sst_m(:,:) = zts(:,:,jp_tem) 
     86write(numout,*) "RSRH sbcssm ub 01", kt, glob_sum(ub(:,:,1)); flush(numout) 
     87write(numout,*) "RSRH sbcssm vb 01", kt, glob_sum(vb(:,:,1)); flush(numout) 
     88write(numout,*) "RSRH sbcssm ssu_m 01", kt, glob_sum(ssu_m(:,:)); flush(numout) 
     89write(numout,*) "RSRH sbcssm ssv_m 01", kt, glob_sum(ssv_m(:,:)); flush(numout) 
     90         IF( l_useCT )  THEN     
     91   sst_m(:,:) = eos_pt_from_ct( zts(:,:,jp_tem), zts(:,:,jp_sal) ) 
     92write(numout,*) "RSRH sbcssm sst_m AA", kt, glob_sum(sst_m(:,:)); flush(numout) 
     93         ELSE                     
     94   sst_m(:,:) = zts(:,:,jp_tem) 
     95write(numout,*) "RSRH sbcssm sst_m BB", kt, glob_sum(sst_m(:,:)); flush(numout) 
    7496         ENDIF 
    7597         sss_m(:,:) = zts(:,:,jp_sal) 
    7698         !                          ! removed inverse barometer ssh when Patm forcing is used (for sea-ice dynamics) 
    77          IF( ln_apr_dyn ) THEN   ;   ssh_m(:,:) = sshn(:,:) - 0.5 * ( ssh_ib(:,:) + ssh_ibb(:,:) ) 
    78          ELSE                    ;   ssh_m(:,:) = sshn(:,:) 
     99         IF( ln_apr_dyn ) THEN    
     100   ssh_m(:,:) = sshn(:,:) - 0.5 * ( ssh_ib(:,:) + ssh_ibb(:,:) ) 
     101         ELSE                     
     102   ssh_m(:,:) = sshn(:,:) 
    79103         ENDIF 
    80104         ! 
     
    84108         ! 
    85109      ELSE 
     110write(numout,*) "RSRH sbcssm nnfsbc NOT 1", kt; flush(numout) 
    86111         !                                                ! ----------------------------------------------- ! 
    87112         IF( kt == nit000 .AND. .NOT. l_ssm_mean ) THEN   !   Initialisation: 1st time-step, no input means ! 
    88             !                                             ! ----------------------------------------------- ! 
     113 write(numout,*) "RSRH sbcssm nnfsbc NOT 1 BB", kt; flush(numout) 
     114           !                                             ! ----------------------------------------------- ! 
    89115            IF(lwp) WRITE(numout,*) 
    90116            IF(lwp) WRITE(numout,*) 'sbc_ssm : mean fields initialised to instantaneous values' 
     
    107133            !                                             ! ---------------------------------------- ! 
    108134         ELSEIF( MOD( kt - 2 , nn_fsbc ) == 0 ) THEN      !   Initialisation: New mean computation   ! 
     135 write(numout,*) "RSRH sbcssm nnfsbc NOT 1 CC", kt; flush(numout) 
    109136            !                                             ! ---------------------------------------- ! 
    110137            ssu_m(:,:) = 0._wp     ! reset to zero ocean mean sbc fields 
     
    119146         !                                                !        Cumulate at each time step        ! 
    120147         !                                                ! ---------------------------------------- ! 
     148write(numout,*) "RSRH sbcssm DD add ub vb", kt; flush(numout) 
     149write(numout,*) "RSRH sbcssm  ub DD", kt, glob_sum(ub(:,:,1)); flush(numout) 
     150write(numout,*) "RSRH sbcssm  vb DD", kt, glob_sum(vb(:,:,1)); flush(numout) 
     151write(numout,*) "RSRH sbcssm  ssu_m DD", kt, glob_sum(ssu_m(:,:)); flush(numout) 
     152write(numout,*) "RSRH sbcssm  ssv_m DD", kt, glob_sum(ssv_m(:,:)); flush(numout) 
     153write(numout,*) "RSRH sbcssm  sst_m DD", kt, glob_sum(sst_m(:,:)); flush(numout) 
    121154         ssu_m(:,:) = ssu_m(:,:) + ub(:,:,1) 
    122155         ssv_m(:,:) = ssv_m(:,:) + vb(:,:,1) 
    123          IF( l_useCT )  THEN    ;   sst_m(:,:) = sst_m(:,:) + eos_pt_from_ct( zts(:,:,jp_tem), zts(:,:,jp_sal) ) 
    124          ELSE                    ;   sst_m(:,:) = sst_m(:,:) + zts(:,:,jp_tem) 
     156write(numout,*) "RSRH ub DD aft", kt, glob_sum(ub(:,:,1)); flush(numout) 
     157write(numout,*) "RSRH vb DD aft", kt, glob_sum(vb(:,:,1)); flush(numout) 
     158write(numout,*) "RSRH ssu_m DD aft", kt, glob_sum(ssu_m(:,:)); flush(numout) 
     159write(numout,*) "RSRH ssv_m DD aft", kt, glob_sum(ssv_m(:,:)); flush(numout) 
     160         IF( l_useCT )  THEN     
     161 write(numout,*) "RSRH sbcssm EE update sst_m", kt; flush(numout) 
     162              sst_m(:,:) = sst_m(:,:) + eos_pt_from_ct( zts(:,:,jp_tem), zts(:,:,jp_sal) ) 
     163         ELSE                     
     164 write(numout,*) "RSRH sbcssm FF update sst_m", kt; flush(numout) 
     165              sst_m(:,:) = sst_m(:,:) + zts(:,:,jp_tem) 
    125166         ENDIF 
    126167         sss_m(:,:) = sss_m(:,:) + zts(:,:,jp_sal) 
     
    145186            e3t_m(:,:) = e3t_m(:,:) * zcoef     ! mean vertical scale factor [m] 
    146187            frq_m(:,:) = frq_m(:,:) * zcoef     ! mean fraction of solar net radiation absorbed in the 1st T level [-] 
    147             ! 
     188write(numout,*) "RSRH ub EE aft", kt, glob_sum(ub(:,:,1)); flush(numout) 
     189write(numout,*) "RSRH vb EE aft", kt, glob_sum(vb(:,:,1)); flush(numout) 
     190write(numout,*) "RSRH ssu_m EE aft", kt, glob_sum(ssu_m(:,:)); flush(numout) 
     191write(numout,*) "RSRH ssv_m EE aft", kt, glob_sum(ssv_m(:,:)); flush(numout) 
     192write(numout,*) "RSRH sst_m EE aft", kt, glob_sum(sst_m(:,:)); flush(numout) 
     193           ! 
    148194         ENDIF 
    149195         !                                                ! ---------------------------------------- ! 
     
    235281         ENDIF 
    236282      ENDIF 
    237       ! 
     283 
     284write(numout,*) "RSRH ub FF aft", glob_sum(ub(:,:,1)); flush(numout) 
     285write(numout,*) "RSRH vb FF aft",  glob_sum(vb(:,:,1)); flush(numout) 
     286write(numout,*) "RSRH ssu_m FF aft",  glob_sum(ssu_m(:,:)); flush(numout) 
     287write(numout,*) "RSRH ssv_m FF aft",  glob_sum(ssv_m(:,:)); flush(numout) 
     288 
     289 
     290     ! 
    238291      IF( .NOT. l_ssm_mean ) THEN   ! default initialisation. needed by iceistate 
    239292         ! 
  • branches/UKMO/dev_r8126_LIM3_couple/NEMOGCM/NEMO/OPA_SRC/TRA/trasbc.F90

    r8877 r8892  
    2727   USE trd_oce        ! trends: ocean variables 
    2828   USE trdtra         ! trends manager: tracers  
     29use lib_fortran 
    2930   ! 
    3031   USE in_out_manager ! I/O manager 
     
    9192         ztrds(:,:,:) = tsa(:,:,:,jp_sal) 
    9293      ENDIF 
     94 
     95write(numout,*) "RSRH tra_sbc tsb, tsa, tsn AA", kt, glob_sum(tsb(:,:,:,1)), glob_sum(tsa(:,:,:,1)),glob_sum(tsn(:,:,:,1)); flush(numout) 
     96 
    9397      ! 
    9498!!gm  This should be moved into sbcmod.F90 module ? (especially now that ln_traqsr is read in namsbc namelist) 
     
    119123         sbc_tsc_b(:,:,:) = sbc_tsc(:,:,:) 
    120124      ENDIF 
     125write(numout,*) "RSRH tra_sbc tsb, tsa, tsn BB", kt, glob_sum(tsb(:,:,:,1)), glob_sum(tsa(:,:,:,1)),glob_sum(tsn(:,:,:,1)); flush(numout) 
    121126      !                             !==  Now sbc tracer content fields  ==! 
    122127      DO jj = 2, jpj 
     
    136141         IF( iom_use('emp_x_sss') )   CALL iom_put( "emp_x_sss", emp (:,:) * tsn(:,:,1,jp_sal) ) 
    137142      ENDIF 
     143write(numout,*) "RSRH tra_sbc tsb, tsa, tsn CC", kt, glob_sum(tsb(:,:,:,1)), glob_sum(tsa(:,:,:,1)),glob_sum(tsn(:,:,:,1)); flush(numout) 
     144write(numout,*) "RSRH tra_sbc qns, emp, sfx FF", kt, glob_sum(qns(:,:)), glob_sum(emp(:,:)), glob_sum(sfx(:,:)); flush(numout) 
     145write(numout,*) "RSRH tra_sbc tsn FF", kt, glob_sum(tsn(:,:,1,1)); flush(numout) 
     146write(numout,*) "RSRH tra_sbc e3t_n FF", kt, glob_sum(e3t_n(:,:,1)); flush(numout) 
     147write(numout,*) "RSRH tra_sbc sbc_tsc_b, sbc_tsc FF", kt, glob_sum(sbc_tsc_b(:,:,1)), glob_sum(sbc_tsc(:,:,1)); flush(numout) 
    138148      ! 
    139149      DO jn = 1, jpts               !==  update tracer trend  ==! 
     
    144154         END DO 
    145155      END DO 
     156write(numout,*) "RSRH tra_sbc tsb, tsa, tsn DD", kt, glob_sum(tsb(:,:,:,1)), glob_sum(tsa(:,:,:,1)),glob_sum(tsn(:,:,:,1)); flush(numout) 
    146157      !                   
    147158      IF( lrst_oce ) THEN           !==  write sbc_tsc in the ocean restart file  ==! 
     
    158169!!gm         probably taken into account in r1_hisf_tbl : to be verified 
    159170      IF( ln_isf ) THEN 
     171write(numout,*) "RSRH tra_sbc risf_tsc_b, risf_tsc FF", kt, glob_sum(risf_tsc_b(:,:,1)), glob_sum(risf_tsc(:,:,1)); flush(numout) 
    160172         zfact = 0.5_wp 
    161173         DO jj = 2, jpj 
     
    183195         END DO 
    184196      END IF 
     197write(numout,*) "RSRH tra_sbc tsb, tsa, tsn EE", kt, glob_sum(tsb(:,:,:,1)), glob_sum(tsa(:,:,:,1)),glob_sum(tsn(:,:,:,1)); flush(numout) 
    185198      ! 
    186199      !---------------------------------------- 
     
    189202      ! 
    190203      IF( ln_rnf ) THEN         ! input of heat and salt due to river runoff  
     204write(numout,*) "RSRH tra_sbc rnf_tsc_b, rnf_tsc FF", kt, glob_sum(rnf_tsc_b(:,:,1)), glob_sum(rnf_tsc(:,:,1)); flush(numout) 
    191205         zfact = 0.5_wp 
    192206         DO jj = 2, jpj  
     
    204218         END DO   
    205219      ENDIF 
     220write(numout,*) "RSRH tra_sbc tsb, tsa, tsn FF", kt, glob_sum(tsb(:,:,:,1)), glob_sum(tsa(:,:,:,1)),glob_sum(tsn(:,:,:,1)); flush(numout) 
    206221 
    207222      IF( iom_use('rnf_x_sst') )   CALL iom_put( "rnf_x_sst", rnf*tsn(:,:,1,jp_tem) )   ! runoff term on sst 
     
    226241         END DO 
    227242      ENDIF 
     243write(numout,*) "RSRH tra_sbc tsb, tsa, tsn GG", kt, glob_sum(tsb(:,:,:,1)), glob_sum(tsa(:,:,:,1)),glob_sum(tsn(:,:,:,1)); flush(numout) 
    228244 
    229245      IF( l_trdtra )   THEN                      ! save the horizontal diffusive trends for further diagnostics 
     
    239255      ! 
    240256      IF( nn_timing == 1 )  CALL timing_stop('tra_sbc') 
     257write(numout,*) "RSRH tra_sbc tsb, tsa, tsn HH", kt, glob_sum(tsb(:,:,:,1)), glob_sum(tsa(:,:,:,1)),glob_sum(tsn(:,:,:,1)); flush(numout) 
    241258      ! 
    242259   END SUBROUTINE tra_sbc 
  • branches/UKMO/dev_r8126_LIM3_couple/NEMOGCM/NEMO/OPA_SRC/step.F90

    r8879 r8892  
    3838   ! 
    3939   USE iom              ! xIOs server 
     40use in_out_manager 
     41use lib_fortran 
     42use oce 
     43 
    4044 
    4145   IMPLICIT NONE 
     
    116120      IF( ln_sto_eos ) CALL sto_pts( tsn  )          ! Random T/S fluctuations 
    117121 
     122write(numout,*) "RSRH step tsb, tsa, tsn AA", kstp, glob_sum(tsb(:,:,:,1)), glob_sum(tsa(:,:,:,1)),glob_sum(tsn(:,:,:,1)); flush(numout) 
     123 
    118124      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    119125      ! Ocean physics update                (ua, va, tsa used as workspace) 
     
    124130                         CALL bn2    ( tsb, rab_b, rn2b ) ! before Brunt-Vaisala frequency 
    125131                         CALL bn2    ( tsn, rab_n, rn2  ) ! now    Brunt-Vaisala frequency 
     132write(numout,*) "RSRH step tsb, tsa, tsn BB", kstp, glob_sum(tsb(:,:,:,1)), glob_sum(tsa(:,:,:,1)),glob_sum(tsn(:,:,:,1)); flush(numout) 
    126133 
    127134      ! 
     
    139146         avmv(:,:,:) = rn_avm0 * wvmask(:,:,:) 
    140147      ENDIF 
     148write(numout,*) "RSRH step tsb, tsa, tsn CC", kstp, glob_sum(tsb(:,:,:,1)), glob_sum(tsa(:,:,:,1)),glob_sum(tsn(:,:,:,1)); flush(numout) 
    141149 
    142150      IF( ln_rnf_mouth ) THEN                         ! increase diffusivity at rivers mouths 
     
    151159                         CALL zdf_mxl( kstp )         ! mixed layer depth 
    152160 
     161write(numout,*) "RSRH step tsb, tsa, tsn DD", kstp, glob_sum(tsb(:,:,:,1)), glob_sum(tsa(:,:,:,1)),glob_sum(tsn(:,:,:,1)); flush(numout) 
    153162                                                      ! write TKE or GLS information in the restart file 
    154163      IF( lrst_oce .AND. lk_zdftke )   CALL tke_rst( kstp, 'WRITE' ) 
     
    173182         ENDIF 
    174183      ENDIF 
     184write(numout,*) "RSRH step tsb, tsa, tsn EE", kstp, glob_sum(tsb(:,:,:,1)), glob_sum(tsa(:,:,:,1)),glob_sum(tsn(:,:,:,1)); flush(numout) 
    175185      !                                                                   ! eddy diffusivity coeff. 
    176186      IF( l_ldftra_time .OR. l_ldfeiv_time )   CALL ldf_tra( kstp )       !       and/or eiv coeff. 
     
    197207               &                                          rhd, gru , grv , grui, grvi   )  ! of t, s, rd at the first ocean level 
    198208!!jc: fs simplification 
     209write(numout,*) "RSRH step tsb, tsa, tsn FF", kstp, glob_sum(tsb(:,:,:,1)), glob_sum(tsa(:,:,:,1)),glob_sum(tsn(:,:,:,1)); flush(numout) 
    199210                             
    200211                         ua(:,:,:) = 0._wp            ! set dynamics trends to zero 
     
    220231                            CALL wzv        ( kstp )              ! now cross-level velocity  
    221232      ENDIF 
     233write(numout,*) "RSRH step tsb, tsa, tsn GG", kstp, glob_sum(tsb(:,:,:,1)), glob_sum(tsa(:,:,:,1)),glob_sum(tsn(:,:,:,1)); flush(numout) 
    222234 
    223235                         CALL dyn_bfr       ( kstp )  ! bottom friction 
     
    237249      IF( lk_diadct  )   CALL dia_dct( kstp )         ! Transports 
    238250                         CALL dia_ar5( kstp )         ! ar5 diag 
     251write(numout,*) "RSRH step tsb, tsa, tsn HH", kstp, glob_sum(tsb(:,:,:,1)), glob_sum(tsa(:,:,:,1)),glob_sum(tsn(:,:,:,1)); flush(numout) 
    239252      IF( lk_diaharm )   CALL dia_harm( kstp )        ! Tidal harmonic analysis 
     253write(numout,*) "RSRH step tsb, tsa, tsn HH.1", kstp, glob_sum(tsb(:,:,:,1)), glob_sum(tsa(:,:,:,1)),glob_sum(tsn(:,:,:,1)); flush(numout) 
    240254                         CALL dia_wri( kstp )         ! ocean model: outputs 
     255write(numout,*) "RSRH step tsb, tsa, tsn HH.2", kstp, glob_sum(tsb(:,:,:,1)), glob_sum(tsa(:,:,:,1)),glob_sum(tsn(:,:,:,1)); flush(numout) 
    241256      ! 
    242257      IF( ln_crs     )   CALL crs_fld       ( kstp )  ! ocean model: online field coarsening & output 
    243258       
     259write(numout,*) "RSRH step tsb, tsa, tsn HH.3", kstp, glob_sum(tsb(:,:,:,1)), glob_sum(tsa(:,:,:,1)),glob_sum(tsn(:,:,:,1)); flush(numout) 
    244260#if defined key_top 
    245261      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     
    247263      !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 
    248264                         CALL trc_stp       ( kstp )  ! time-stepping 
     265write(numout,*) "RSRH step tsb, tsa, tsn HH.4", kstp, glob_sum(tsb(:,:,:,1)), glob_sum(tsa(:,:,:,1)),glob_sum(tsn(:,:,:,1)); flush(numout) 
    249266#endif 
    250267 
     
    257274         & ln_trainc )   CALL tra_asm_inc   ( kstp )  ! apply tracer assimilation increment 
    258275                         CALL tra_sbc       ( kstp )  ! surface boundary condition 
     276write(numout,*) "RSRH step tsb, tsa, tsn HH.5", kstp, glob_sum(tsb(:,:,:,1)), glob_sum(tsa(:,:,:,1)),glob_sum(tsn(:,:,:,1)); flush(numout) 
    259277      IF( ln_traqsr  )   CALL tra_qsr       ( kstp )  ! penetrative solar radiation qsr 
     278write(numout,*) "RSRH step tsb, tsa, tsn HH.6", kstp, glob_sum(tsb(:,:,:,1)), glob_sum(tsa(:,:,:,1)),glob_sum(tsn(:,:,:,1)); flush(numout) 
    260279      IF( ln_trabbc  )   CALL tra_bbc       ( kstp )  ! bottom heat flux 
     280write(numout,*) "RSRH step tsb, tsa, tsn HH.7", kstp, glob_sum(tsb(:,:,:,1)), glob_sum(tsa(:,:,:,1)),glob_sum(tsn(:,:,:,1)); flush(numout) 
    261281      IF( lk_trabbl  )   CALL tra_bbl       ( kstp )  ! advective (and/or diffusive) bottom boundary layer scheme 
     282write(numout,*) "RSRH step tsb, tsa, tsn HH.8", kstp, glob_sum(tsb(:,:,:,1)), glob_sum(tsa(:,:,:,1)),glob_sum(tsn(:,:,:,1)); flush(numout) 
    262283      IF( ln_tradmp  )   CALL tra_dmp       ( kstp )  ! internal damping trends 
     284write(numout,*) "RSRH step tsb, tsa, tsn HH.9", kstp, glob_sum(tsb(:,:,:,1)), glob_sum(tsa(:,:,:,1)),glob_sum(tsn(:,:,:,1)); flush(numout) 
    263285      IF( ln_bdy     )   CALL bdy_tra_dmp   ( kstp )  ! bdy damping trends 
     286write(numout,*) "RSRH step tsb, tsa, tsn HH.10", kstp, glob_sum(tsb(:,:,:,1)), glob_sum(tsa(:,:,:,1)),glob_sum(tsn(:,:,:,1)); flush(numout) 
    264287#if defined key_agrif 
    265288      IF(.NOT. Agrif_Root())  &  
    266289               &         CALL Agrif_Sponge_tra        ! tracers sponge 
    267290#endif 
    268                          CALL tra_adv       ( kstp )  ! horizontal & vertical advection 
     291 write(numout,*) "RSRH step tsb, tsa, tsn II", kstp, glob_sum(tsb(:,:,:,1)), glob_sum(tsa(:,:,:,1)),glob_sum(tsn(:,:,:,1)); flush(numout) 
     292                        CALL tra_adv       ( kstp )  ! horizontal & vertical advection 
    269293                         CALL tra_ldf       ( kstp )  ! lateral mixing 
    270294 
     
    298322      ! 
    299323      IF( ln_diahsb        )   CALL dia_hsb( kstp )         ! - ML - global conservation diagnostics 
     324 write(numout,*) "RSRH step tsb, tsa, tsn JJ", kstp, glob_sum(tsb(:,:,:,1)), glob_sum(tsa(:,:,:,1)),glob_sum(tsn(:,:,:,1)); flush(numout) 
    300325 
    301326!!gm : This does not only concern the dynamics ==>>> add a new title 
     
    328353                         CALL ctl_stop( 'step: indic < 0' ) 
    329354                         CALL dia_wri_state( 'output.abort', kstp ) 
    330                                CALL ctl_stop('STOP','NEMO failure in stp') 
    331355      ENDIF 
    332356      IF( kstp == nit000 ) THEN 
     
    336360            &    CALL FLUSH    ( numoni )     ! flush output namelist ice (if exist) 
    337361      ENDIF 
     362 write(numout,*) "RSRH step tsb, tsa, tsn KK", kstp, glob_sum(tsb(:,:,:,1)), glob_sum(tsa(:,:,:,1)),glob_sum(tsn(:,:,:,1)); flush(numout) 
    338363 
    339364      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    340365      ! Coupled mode 
    341366      !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 
    342       !IF( lk_oasis         )   CALL sbc_cpl_snd( kstp )     ! coupled mode : field exchanges 
     367!!gm why lk_oasis and not lk_cpl ???? 
     368      !IF( lk_oasis   )   CALL sbc_cpl_snd( kstp )     ! coupled mode : field exchanges 
    343369      ! 
    344370#if defined key_iomput 
Note: See TracChangeset for help on using the changeset viewer.