Changeset 16 for trunk/NEMO/OPA_SRC/DTA/dtatem.F90

Timestamp:

2004-02-17T09:06:15+01:00 (20 years ago)

Author:

opalod

Message:

CT : UPDATE001 : First major NEMO update

File:

• trunk/NEMO/OPA_SRC/DTA/dtatem.F90 (modified) (5 diffs)

trunk/NEMO/OPA_SRC/DTA/dtatem.F90

@@ -23 +23 @@
 
    !! * Shared module variables
-   LOGICAL , PUBLIC, PARAMETER ::   lk_dtatem = .TRUE.   ! temperature data flag
-   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) ::   &
-      t_dta             ! temperature data at given time-step
+   LOGICAL , PUBLIC, PARAMETER ::   lk_dtatem = .TRUE.   !: temperature data flag
+   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) ::   &  !:
+      t_dta             !: temperature data at given time-step
 
    !! * Module variables
@@ -91 +91 @@
       INTEGER ::   &
          imois, iman, itime, ik ,    &  ! temporary integers
-         i15, ipi, ipj, ipk             !    "          "
+         i15, ipi, ipj, ipk,         &  !    "          "
+         il0, il1, ii0, ii1, ij0, ij1   !    "          "
 
       INTEGER, DIMENSION(jpmois) ::   istep
@@ -192 +193 @@
             !                                        ! ======================= 
 
-            DO jj = mj0(101), mj1(109)                   ! Reduced temperature at Alboran Sea
-               DO ji = mi0(141), mi1(155)
+            ij0 = 101   ;   ij1 = 109
+            ii0 = 141   ;   ii1 = 155
+            DO jj = mj0(ij0), mj1(ij1)                      ! Reduced temperature in the Alboran Sea
+               DO ji = mi0(ii0), mi1(ii1)
                   temdta(ji,jj, 13:13 ,:) = temdta(ji,jj, 13:13 ,:) - 0.20
                   temdta(ji,jj, 14:15 ,:) = temdta(ji,jj, 14:15 ,:) - 0.35
@@ -202 +205 @@
             IF( n_cla == 0 ) THEN 
                !                                         ! Reduced temperature at Red Sea
-               temdta( mi0(148):mi1(160) , mj0(87):mj1(96) ,  4:10 , : ) = 7.0 
-               temdta( mi0(148):mi1(160) , mj0(87):mj1(96) , 11:13 , : ) = 6.5 
-               temdta( mi0(148):mi1(160) , mj0(87):mj1(96) , 14:20 , : ) = 6.0
+               ij0 =  87   ;   ij1 =  96
+               ii0 = 148   ;   ii1 = 160
+               temdta( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ,  4:10 , : ) = 7.0 
+               temdta( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 11:13 , : ) = 6.5 
+               temdta( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 14:20 , : ) = 6.0
             ELSE
-               DO jl = mi0(138), mi1(138)                ! New temperature profile at Gibraltar
-                  DO jj = mj0(101), mj1(102)
-                     DO ji = mi0(139), mi1(139)
+               il0 = 138   ;   il1 = 138
+               ij0 = 101   ;   ij1 = 102
+               ii0 = 139   ;   ii1 = 139
+               DO jl = mi0(ii0), mi1(ii1)                ! New temperature profile at Gibraltar
+                  DO jj = mj0(ij0), mj1(ij1)
+                     DO ji = mi0(ii0), mi1(ii1)
                         temdta(ji,jj,:,:) = temdta(jl,jj,:,:)
                      END DO
                   END DO
                END DO
-               DO jl = mi0(164), mi1(164)                ! New temperature profile at Bab el Mandeb
-                  DO jj = mj0(88), mj1(88)
-                     DO ji = mi0(161), mi1(163)
+               il0 = 164   ;   il1 = 164
+               ij0 =  88   ;   ij1 =  88
+               ii0 = 161   ;   ii1 = 163
+               DO jl = mi0(ii0), mi1(ii1)                ! New temperature profile at Bab el Mandeb
+                  DO jj = mj0(ij0), mj1(ij1)
+                     DO ji = mi0(ii0), mi1(ii1)
                         temdta(ji,jj,:,:) = temdta(jl,jj,:,:)
                      END DO
                   END DO
-                  DO jj = mj0(87), mj1(87)
-                     DO ji = mi0(161), mi1(163)
+                  ij0 =  87   ;   ij1 =  87
+                  DO jj = mj0(ij0), mj1(ij1)
+                     DO ji = mi0(ii0), mi1(ii1)
                         temdta(ji,jj,:,:) = temdta(jl,jj,:,:)
                      END DO
@@ -274 +286 @@
    !!   Default case                           NO 3D temperature data field
    !!----------------------------------------------------------------------
-   LOGICAL , PUBLIC, PARAMETER ::   lk_dtatem = .FALSE.   ! temperature data flag
+   LOGICAL , PUBLIC, PARAMETER ::   lk_dtatem = .FALSE.   !: temperature data flag
 CONTAINS
    SUBROUTINE dta_tem( kt )        ! Empty routine
-      WRITE(*,*) kt
+      WRITE(*,*) 'dta_tem: You should not have seen this print! error?', kt
    END SUBROUTINE dta_tem
 #endif