Changeset 3372 for branches/2012/dev_r3337_NOCS10_ICB/NEMOGCM/NEMO/OPA_SRC/ICB/icbthm.F90

Timestamp:

2012-04-30T12:50:36+02:00 (12 years ago)

Author:

sga

Message:

NEMO branch dev_r3337_NOCS10_ICB: change all routine names and create more Gurvanistic havoc

File:

• branches/2012/dev_r3337_NOCS10_ICB/NEMOGCM/NEMO/OPA_SRC/ICB/icbthm.F90 (modified) (7 diffs)

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

@@ -12 +12 @@
    !!----------------------------------------------------------------------
    !!----------------------------------------------------------------------
-   !!   thermodynamics : initialise
-   !!                    reference for equations - M = Martin + Adcroft, OM 34, 2010
+   !!   icb_thm : initialise
+   !!             reference for equations - M = Martin + Adcroft, OM 34, 2010
    !!----------------------------------------------------------------------
    USE par_oce        ! NEMO parameters
@@ -29 +29 @@
    PRIVATE
 
-   PUBLIC   thermodynamics ! routine called in icbrun.F90 module
+   PUBLIC   icb_thm ! routine called in icbrun.F90 module
 
 CONTAINS
 
-   SUBROUTINE thermodynamics( kt )
+   SUBROUTINE icb_thm( kt )
       !!----------------------------------------------------------------------
-      !!                  ***  ROUTINE thermodynamics  ***
+      !!                  ***  ROUTINE icb_thm  ***
       !!
       !! ** Purpose :   compute the iceberg thermodynamics.
@@ -41 +41 @@
       !! ** Method  : - blah blah
       !!----------------------------------------------------------------------
-      INTEGER, INTENT(in) ::   kt   ! timestep number, just passed to print_berg
+      INTEGER, INTENT(in) ::   kt   ! timestep number, just passed to icb_utl_print_berg
       !
       INTEGER  ::   ii, ij
@@ -68 +68 @@
          pt => this%current_point
          nknberg = this%number(1)
-         CALL interp_flds( pt%xi, pt%e1, pt%uo, pt%ui, pt%ua, pt%ssh_x, &
-         &                 pt%yj, pt%e2, pt%vo, pt%vi, pt%va, pt%ssh_y, &
+         CALL icb_utl_interp( pt%xi, pt%e1, pt%uo, pt%ui, pt%ua, pt%ssh_x, &
+         &                    pt%yj, pt%e2, pt%vo, pt%vi, pt%va, pt%ssh_y, &
          &                 pt%sst, pt%cn, pt%hi, zff )
          !
@@ -159 +159 @@
             zheat = zmelt * pt%heat_density              ! kg/s x J/kg = J/s
             berg_grid%calving_hflx (ii,ij) = berg_grid%calving_hflx (ii,ij) + zheat    * z1_e1e2    ! W/m2
-            CALL melt_budget( ii, ij, zMnew, zheat, this%mass_scaling,       &
-            &                         zdM, zdMbitsE, zdMbitsM, zdMb, zdMe,   &
-            &                         zdMv, z1_dt_e1e2 )
+            CALL icb_dia_melt( ii, ij, zMnew, zheat, this%mass_scaling,       &
+            &                          zdM, zdMbitsE, zdMbitsM, zdMb, zdMe,   &
+            &                          zdMv, z1_dt_e1e2 )
          ELSE
-            WRITE(numout,*) 'thermodynamics: berg ',this%number(:),' appears to have grounded  at ',narea,ii,ij
-            CALL print_berg( this, kt )
+            WRITE(numout,*) 'icb_thm: berg ',this%number(:),' appears to have grounded  at ',narea,ii,ij
+            CALL icb_utl_print_berg( this, kt )
             WRITE(numout,*) 'msk=',tmask(ii,ij,1), e1e2t(ii,ij)
-            CALL ctl_stop('thermodynamics', 'berg appears to have grounded!')
+            CALL ctl_stop('icb_thm', 'berg appears to have grounded!')
          ENDIF
 
@@ -189 +189 @@
 
          IF( zMnew <= 0._wp ) THEN       ! Delete the berg if completely melted
-            CALL delete_iceberg_from_list( first_berg, this )
+            CALL icb_utl_delete( first_berg, this )
             !
          ELSE                            ! Diagnose mass distribution on grid
             z1_e1e2 = 1._wp / e1e2t(ii,ij) * this%mass_scaling
-            CALL size_budget( ii, ij, zWn, zLn, zAbits,   &
-            &                 this%mass_scaling, zMnew, znMbits, z1_e1e2)
+            CALL icb_dia_size( ii, ij, zWn, zLn, zAbits,   &
+            &                  this%mass_scaling, zMnew, znMbits, z1_e1e2)
          ENDIF
          !
@@ -209 +209 @@
       ENDIF
       !
-   END SUBROUTINE thermodynamics
+   END SUBROUTINE icb_thm
 
    !!======================================================================