Changeset 3370 for branches/2012/dev_r3337_NOCS10_ICB/NEMOGCM/NEMO/OPA_SRC/ICB/icbclv.F90

Timestamp:

2012-04-30T10:27:44+02:00 (12 years ago)

Author:

sga

Message:

NEMO branch dev_r3337_NOCS10_ICB: lots of cosmetic Gurvanistic changes (the odd space or exclamation mark!)

File:

• branches/2012/dev_r3337_NOCS10_ICB/NEMOGCM/NEMO/OPA_SRC/ICB/icbclv.F90 (modified) (11 diffs)

branches/2012/dev_r3337_NOCS10_ICB/NEMOGCM/NEMO/OPA_SRC/ICB/icbclv.F90

@@ -3 +3 @@
    !!======================================================================
    !!                       ***  MODULE  icbclv  ***
-   !! Ocean physics:  calving routines for iceberg calving
+   !! Icebergs:  calving routines for iceberg calving
    !!======================================================================
    !! History : 3.3.1 !  2010-01  (Martin&Adcroft) Original code
@@ -12 +12 @@
    !!----------------------------------------------------------------------
    !!----------------------------------------------------------------------
-   !!   accumulate_calving      :
-   !!   icb_gen       : generate test icebergs
+   !!   accumulate_calving :
+   !!   icb_gen            : generate test icebergs
    !!----------------------------------------------------------------------
    USE par_oce        ! NEMO parameters
@@ -21 +21 @@
    USE lbclnk         ! NEMO boundary exchanges for gridded data
 
-   USE icb_oce        ! define iceberg arrays
-   USE icbdia         ! iceberg utility routines
+   USE icb_oce        ! iceberg variables
+   USE icbdia         ! iceberg diagnostics
    USE icbutl         ! iceberg utility routines
 
@@ -31 +31 @@
    PUBLIC   calve_icebergs      ! routine called in icbrun.F90 module
 
+   !!----------------------------------------------------------------------
+   !! NEMO/OPA 3.3 , NEMO Consortium (2011)
+   !! $Id:$
+   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
+   !!----------------------------------------------------------------------
 CONTAINS
 
@@ -41 +46 @@
       !! ** input   : - ?
       !!----------------------------------------------------------------------
-      !
       INTEGER, INTENT(in)             :: kt
       !
@@ -63 +67 @@
          !do jn=1, nclasses
          !  where (berg_grid%calving==0.) berg_grid%stored_ice(:,:,jn)=0.
-         !enddo
+         !end do
          DO jj = 2, jpjm1
             DO ji = 2, jpim1
@@ -109 +113 @@
       !!                is considered
       !!----------------------------------------------------------------------
+      INTEGER       ::   ji, jj, jn   ! dummy loop indices
+      INTEGER       ::   icnt, icntmax
+      TYPE(iceberg) ::   newberg
+      TYPE(point)   ::   newpt
+      REAL(wp)      ::   zddt, zcalved_to_berg, zheat_to_berg
+      !!----------------------------------------------------------------------
       !
-      INTEGER                         ::   icnt,icntmax
-      TYPE(iceberg)                   ::   newberg
-      TYPE(point)                     ::   newpt
-      INTEGER                         ::   ji, jj, jn                             ! dummy loop indices
-      REAL(wp)                        ::   zddt, zcalved_to_berg, zheat_to_berg
-      !!----------------------------------------------------------------------
-
-      icntmax          = 0
-
+      icntmax = 0
+      !
       DO jn = 1, nclasses
          DO jj = nicbdj, nicbej
@@ -128 +131 @@
                DO WHILE (berg_grid%stored_ice(ji,jj,jn) >= rn_initial_mass(jn) * rn_mass_scaling(jn) )
                   !
-                  newpt%lon = glamt(ji,jj)       ! at t-point (centre of the cell)
+                  newpt%lon = glamt(ji,jj)         ! at t-point (centre of the cell)
                   newpt%lat = gphit(ji,jj)
                   newpt%xi  = REAL( nimpp+ji-1, wp )
                   newpt%yj  = REAL( njmpp+jj-1, wp )
                   !
-                  newpt%uvel = 0._wp             ! initially at rest
+                  newpt%uvel = 0._wp               ! initially at rest
                   newpt%vvel = 0._wp
                   !                                ! set berg characteristics
-                  newpt%mass         = rn_initial_mass     (jn)
-                  newpt%thickness    = rn_initial_thickness(jn)
-                  newpt%width        = first_width    (jn)
-                  newpt%length       = first_length   (jn)
+                  newpt%mass           = rn_initial_mass     (jn)
+                  newpt%thickness      = rn_initial_thickness(jn)
+                  newpt%width          = first_width         (jn)
+                  newpt%length         = first_length        (jn)
                   newberg%mass_scaling = rn_mass_scaling     (jn)
-                  newpt%mass_of_bits = 0._wp                          ! no bergy
+                  newpt%mass_of_bits   = 0._wp                          ! no bergy
                   !
                   newpt%year   = current_year
@@ -157 +160 @@
                   berg_grid%stored_heat(ji,jj) = berg_grid%stored_heat(ji,jj) - zheat_to_berg
                   !                                ! Stored mass
-                  berg_grid%stored_ice(ji,jj,jn)   = berg_grid%stored_ice(ji,jj,jn) - zcalved_to_berg
+                  berg_grid%stored_ice(ji,jj,jn) = berg_grid%stored_ice(ji,jj,jn) - zcalved_to_berg
                   !
                   zddt  = zddt + berg_dt * 2._wp / 17._wp    ! Minor offset to start day    !!gm why???
@@ -171 +174 @@
       DO jn = 1,nclasses
          CALL lbc_lnk( berg_grid%stored_ice(:,:,jn), 'T', 1._wp )
-      ENDDO
+      END DO
       CALL lbc_lnk( berg_grid%stored_heat, 'T', 1._wp )
       !
@@ -178 +181 @@
    END SUBROUTINE  calve_icebergs
 
+   !!======================================================================
 END MODULE icbclv