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 867 for trunk/NEMO/LIM_SRC_3/limitd_me.F90 – NEMO

Ignore:
Timestamp:
2008-03-14T15:57:19+01:00 (16 years ago)
Author:
ctlod
Message:

Add control prints for sea-ice

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/LIM_SRC_3/limitd_me.F90

    r863 r867  
    462462      d_smv_i_trp(:,:,:)  = smv_i(:,:,:) - old_smv_i(:,:,:) 
    463463 
    464       !-------------------------! 
    465       ! Back to initial values 
    466       !-------------------------! 
    467  
    468       ! update of fields will be made later in lim update 
    469       u_ice(:,:)    = old_u_ice(:,:) 
    470       v_ice(:,:)    = old_v_ice(:,:) 
    471       a_i(:,:,:)    = old_a_i(:,:,:) 
    472       v_s(:,:,:)    = old_v_s(:,:,:) 
    473       v_i(:,:,:)    = old_v_i(:,:,:) 
    474       e_s(:,:,:,:)  = old_e_s(:,:,:,:) 
    475       e_i(:,:,:,:)  = old_e_i(:,:,:,:) 
    476       oa_i(:,:,:)   = old_oa_i(:,:,:) 
    477       IF ( ( num_sal .EQ. 2 ) .OR. ( num_sal .EQ. 4 ) ) &  
    478       smv_i(:,:,:)  = old_smv_i(:,:,:) 
    479  
    480464      IF(ln_ctl) THEN     ! Control print 
     465         CALL prt_ctl_info(' ') 
     466         CALL prt_ctl_info(' - Cell values : ') 
     467         CALL prt_ctl_info('   ~~~~~~~~~~~~~ ') 
    481468         CALL prt_ctl(tab2d_1=area , clinfo1=' lim_itd_me  : cell area :') 
    482469         CALL prt_ctl(tab2d_1=at_i , clinfo1=' lim_itd_me  : at_i      :') 
     
    484471         CALL prt_ctl(tab2d_1=vt_s , clinfo1=' lim_itd_me  : vt_s      :') 
    485472         DO jl = 1, jpl 
     473            CALL prt_ctl_info(' ') 
    486474            CALL prt_ctl_info(' - Category : ', ivar1=jl) 
    487475            CALL prt_ctl_info('   ~~~~~~~~~~') 
     
    497485            CALL prt_ctl(tab2d_1=smv_i (:,:,jl)   , clinfo1= ' lim_itd_me  : smv_i    : ') 
    498486            DO jk = 1, nlay_i 
     487               CALL prt_ctl_info(' ') 
    499488               CALL prt_ctl_info(' - Layer : ', ivar1=jk) 
    500489               CALL prt_ctl_info('   ~~~~~~~') 
     
    505494      ENDIF 
    506495 
     496      !-------------------------! 
     497      ! Back to initial values 
     498      !-------------------------! 
     499 
     500      ! update of fields will be made later in lim update 
     501      u_ice(:,:)    = old_u_ice(:,:) 
     502      v_ice(:,:)    = old_v_ice(:,:) 
     503      a_i(:,:,:)    = old_a_i(:,:,:) 
     504      v_s(:,:,:)    = old_v_s(:,:,:) 
     505      v_i(:,:,:)    = old_v_i(:,:,:) 
     506      e_s(:,:,:,:)  = old_e_s(:,:,:,:) 
     507      e_i(:,:,:,:)  = old_e_i(:,:,:,:) 
     508      oa_i(:,:,:)   = old_oa_i(:,:,:) 
     509      IF ( ( num_sal .EQ. 2 ) .OR. ( num_sal .EQ. 4 ) ) &  
     510      smv_i(:,:,:)  = old_smv_i(:,:,:) 
    507511 
    508512      !----------------------------------------------------! 
Note: See TracChangeset for help on using the changeset viewer.