Changeset 494 for trunk/NEMO/TOP_SRC/TRP/trcbbc.F90

Timestamp:

2006-09-01T16:03:49+02:00 (18 years ago)

Author:

opalod

Message:

nemo_v1_update_062:CE+RB: use IOM for passive tracers

File:

• trunk/NEMO/TOP_SRC/TRP/trcbbc.F90 (modified) (3 diffs)

trunk/NEMO/TOP_SRC/TRP/trcbbc.F90

@@ -142 +142 @@
       !!----------------------------------------------------------------------
       !! * Modules used
-      USE ioipsl
+      USE iom
 
       !! * local declarations
@@ -148 +148 @@
       INTEGER  ::   ji, jj              ! dummy loop indices
       INTEGER  ::   inum = 11           ! temporary logical unit
-      INTEGER  ::   itime               ! temporary integers
-      REAL(wp) ::   zdate0, zdt         ! temporary scalars
-      REAL(wp), DIMENSION(1) :: zdept   ! temporary workspace
-      REAL(wp), DIMENSION(jpidta,jpjdta) ::   &
-         zlamt, zphit, zdta   ! temporary workspace
 
       NAMELIST/namtrcbbc/ngeo_trc_flux, ngeo_trc_flux_const 
@@ -193 +188 @@
       CASE ( 2 )                ! variable geothermal heat flux
          ! read the geothermal fluxes in mW/m2
-         clname = 'passivetrc_geothermal_heating'
-         itime = 1
-         zlamt(:,:) = 0.
-         zphit(:,:) = 0.
-         IF(lwp) WRITE(numout,*) '             *** variable geothermal heat flux read in ', clname, ' file'
-         CALL restini( clname, jpidta, jpjdta, zlamt, zphit, 1, zdept , clname,   &
-                       itime, zdate0, zdt, inum , domain_id=nidom )
-         CALL restget( inum, 'heatflow', jpidta, jpjdta, 1, 0, .FALSE., zdta )
-         DO jj = 1, nlcj
-            DO ji = 1, nlci
-              qgh_trd(ji,jj) = zdta(mig(ji),mjg(jj))
-            END DO
-         END DO
-
-         CALL restclo( inum )
+         CALL iom_open ( 'geothermal_heating_trc.nc', inum )
+         CALL iom_get ( inum, jpdom_data, 'heatflow', qgh_trd )
+         CALL iom_close (inum)
+
          qgh_trd(:,:) = qgh_trd(:,:) * 1.e-3 ! conversion in W/m2
 
       CASE DEFAULT
-         IF(lwp) WRITE(numout,cform_err)
-         IF(lwp) WRITE(numout,*) '     bad flag value for ngeo_trc_flux = ', ngeo_trc_flux
-         nstop = nstop + 1
+         WRITE(ctmp1,*) '     bad flag value for ngeo_flux = ', ngeo_flux
+         CALL ctl_stop( ctmp1 )
 
       END SELECT