Changeset 867


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

Add control prints for sea-ice

Location:
trunk/NEMO/LIM_SRC_3
Files:
6 edited

Legend:

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

    r863 r867  
    271271 
    272272      IF(ln_ctl) THEN   ! Control print 
     273         CALL prt_ctl_info(' ') 
     274         CALL prt_ctl_info(' - Cell values : ') 
     275         CALL prt_ctl_info('   ~~~~~~~~~~~~~ ') 
    273276         CALL prt_ctl(tab2d_1=tio_u     , clinfo1=' lim_dyn  : tio_u     :', tab2d_2=tio_v , clinfo2=' tio_v :') 
    274277         CALL prt_ctl(tab2d_1=ust2s     , clinfo1=' lim_dyn  : ust2s     :') 
     
    284287         CALL prt_ctl(tab2d_1=stress12_i, clinfo1=' lim_dyn  : stress12_i:') 
    285288         DO jl = 1, jpl 
     289            CALL prt_ctl_info(' ') 
    286290            CALL prt_ctl_info(' - Category : ', ivar1=jl) 
    287291            CALL prt_ctl_info('   ~~~~~~~~~~') 
     
    297301            CALL prt_ctl(tab2d_1=smv_i (:,:,jl)   , clinfo1= ' lim_dyn  : smv_i    : ') 
    298302            DO ja = 1, nlay_i 
     303               CALL prt_ctl_info(' ') 
    299304               CALL prt_ctl_info(' - Layer : ', ivar1=ja) 
    300305               CALL prt_ctl_info('   ~~~~~~~') 
  • 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      !----------------------------------------------------! 
  • trunk/NEMO/LIM_SRC_3/limitd_th.F90

    r863 r867  
    134134      d_smv_i_thd(:,:,:) = smv_i(:,:,:) - old_smv_i(:,:,:) 
    135135 
    136       !- Recover Old values 
    137       a_i(:,:,:)         = old_a_i (:,:,:) 
    138       v_s(:,:,:)         = old_v_s (:,:,:) 
    139       v_i(:,:,:)         = old_v_i (:,:,:) 
    140       e_s(:,:,:,:)       = old_e_s (:,:,:,:) 
    141       e_i(:,:,:,:)       = old_e_i (:,:,:,:) 
    142  
    143       IF ( ( num_sal .EQ. 2 ) .OR. ( num_sal .EQ. 4 ) ) & 
    144       smv_i(:,:,:)       = old_smv_i (:,:,:) 
    145  
    146136      IF(ln_ctl) THEN   ! Control print 
     137         CALL prt_ctl_info(' ') 
     138         CALL prt_ctl_info(' - Cell values : ') 
     139         CALL prt_ctl_info('   ~~~~~~~~~~~~~ ') 
    147140         CALL prt_ctl(tab2d_1=area , clinfo1=' lim_itd_th  : cell area :') 
    148141         CALL prt_ctl(tab2d_1=at_i , clinfo1=' lim_itd_th  : at_i      :') 
     
    150143         CALL prt_ctl(tab2d_1=vt_s , clinfo1=' lim_itd_th  : vt_s      :') 
    151144         DO jl = 1, jpl 
     145            CALL prt_ctl_info(' ') 
    152146            CALL prt_ctl_info(' - Category : ', ivar1=jl) 
    153147            CALL prt_ctl_info('   ~~~~~~~~~~') 
     
    163157            CALL prt_ctl(tab2d_1=smv_i (:,:,jl)   , clinfo1= ' lim_itd_th  : smv_i    : ') 
    164158            DO ja = 1, nlay_i 
     159               CALL prt_ctl_info(' ') 
    165160               CALL prt_ctl_info(' - Layer : ', ivar1=ja) 
    166161               CALL prt_ctl_info('   ~~~~~~~') 
     
    171166      ENDIF 
    172167       
     168      !- Recover Old values 
     169      a_i(:,:,:)         = old_a_i (:,:,:) 
     170      v_s(:,:,:)         = old_v_s (:,:,:) 
     171      v_i(:,:,:)         = old_v_i (:,:,:) 
     172      e_s(:,:,:,:)       = old_e_s (:,:,:,:) 
     173      e_i(:,:,:,:)       = old_e_i (:,:,:,:) 
     174 
     175      IF ( ( num_sal .EQ. 2 ) .OR. ( num_sal .EQ. 4 ) ) & 
     176      smv_i(:,:,:)       = old_smv_i (:,:,:) 
     177 
    173178 
    174179      END SUBROUTINE lim_itd_th 
  • trunk/NEMO/LIM_SRC_3/limthd.F90

    r863 r867  
    476476 
    477477      IF(ln_ctl) THEN   ! Control print 
     478         CALL prt_ctl_info(' ') 
     479         CALL prt_ctl_info(' - Cell values : ') 
     480         CALL prt_ctl_info('   ~~~~~~~~~~~~~ ') 
    478481         CALL prt_ctl(tab2d_1=area , clinfo1=' lim_thd  : cell area :') 
    479482         CALL prt_ctl(tab2d_1=at_i , clinfo1=' lim_thd  : at_i      :') 
     
    481484         CALL prt_ctl(tab2d_1=vt_s , clinfo1=' lim_thd  : vt_s      :') 
    482485         DO jl = 1, jpl 
     486            CALL prt_ctl_info(' ') 
    483487            CALL prt_ctl_info(' - Category : ', ivar1=jl) 
    484488            CALL prt_ctl_info('   ~~~~~~~~~~') 
     
    494498            CALL prt_ctl(tab2d_1=smv_i (:,:,jl)   , clinfo1= ' lim_thd  : smv_i    : ') 
    495499            DO jk = 1, nlay_i 
     500               CALL prt_ctl_info(' ') 
    496501               CALL prt_ctl_info(' - Layer : ', ivar1=jk) 
    497502               CALL prt_ctl_info('   ~~~~~~~') 
  • trunk/NEMO/LIM_SRC_3/limtrp.F90

    r863 r867  
    554554 
    555555      IF(ln_ctl) THEN   ! Control print 
     556         CALL prt_ctl_info(' ') 
     557         CALL prt_ctl_info(' - Cell values : ') 
     558         CALL prt_ctl_info('   ~~~~~~~~~~~~~ ') 
    556559         CALL prt_ctl(tab2d_1=area , clinfo1=' lim_trp  : cell area :') 
    557560         CALL prt_ctl(tab2d_1=at_i , clinfo1=' lim_trp  : at_i      :') 
     
    559562         CALL prt_ctl(tab2d_1=vt_s , clinfo1=' lim_trp  : vt_s      :') 
    560563         DO jl = 1, jpl 
     564            CALL prt_ctl_info(' ') 
    561565            CALL prt_ctl_info(' - Category : ', ivar1=jl) 
    562566            CALL prt_ctl_info('   ~~~~~~~~~~') 
     
    572576            CALL prt_ctl(tab2d_1=smv_i (:,:,jl)   , clinfo1= ' lim_trp  : smv_i    : ') 
    573577            DO jk = 1, nlay_i 
     578               CALL prt_ctl_info(' ') 
    574579               CALL prt_ctl_info(' - Layer : ', ivar1=jk) 
    575580               CALL prt_ctl_info('   ~~~~~~~') 
  • trunk/NEMO/LIM_SRC_3/limupdate.F90

    r863 r867  
    998998 
    999999      IF(ln_ctl) THEN   ! Control print 
     1000         CALL prt_ctl_info(' ') 
    10001001         CALL prt_ctl_info(' - Cell values : ') 
    10011002         CALL prt_ctl_info('   ~~~~~~~~~~~~~ ') 
     
    10101011 
    10111012         DO jl = 1, jpl 
     1013            CALL prt_ctl_info(' ') 
    10121014            CALL prt_ctl_info(' - Category : ', ivar1=jl) 
    10131015            CALL prt_ctl_info('   ~~~~~~~~~~') 
     
    10541056            DO jk = 1, nlay_i 
    10551057               CALL prt_ctl_info(' - Layer : ', ivar1=jk) 
    1056                CALL prt_ctl(tab2d_1=t_i(:,:,jk,jl) , clinfo1= ' lim_dyn  : t_i       : ') 
     1058               CALL prt_ctl(tab2d_1=t_i(:,:,jk,jl) , clinfo1= ' lim_update  : t_i       : ') 
    10571059            END DO 
    10581060         END DO 
    10591061 
     1062         CALL prt_ctl_info(' ') 
    10601063         CALL prt_ctl_info(' - Heat / FW fluxes : ') 
    10611064         CALL prt_ctl_info('   ~~~~~~~~~~~~~~~~~~ ') 
     
    10641067         CALL prt_ctl(tab2d_1=fhbri  , clinfo1= ' lim_update : fhbri : ', tab2d_2=fheat_rpo , clinfo2= ' fheat_rpo : ') 
    10651068 
     1069         CALL prt_ctl_info(' ') 
    10661070         CALL prt_ctl_info(' - Stresses : ') 
    10671071         CALL prt_ctl_info('   ~~~~~~~~~~ ') 
Note: See TracChangeset for help on using the changeset viewer.