New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 16 for trunk/NEMO/OPA_SRC/DTA – NEMO

Changeset 16 for trunk/NEMO/OPA_SRC/DTA


Ignore:
Timestamp:
2004-02-17T09:06:15+01:00 (20 years ago)
Author:
opalod
Message:

CT : UPDATE001 : First major NEMO update

Location:
trunk/NEMO/OPA_SRC/DTA
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/OPA_SRC/DTA/dtasal.F90

    r3 r16  
    2323    
    2424   !! * Shared module variables 
    25    LOGICAL , PUBLIC, PARAMETER ::   lk_dtasal = .TRUE.    ! salinity data flag 
    26    REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) ::   & 
    27       s_dta       ! salinity data at given time-step 
     25   LOGICAL , PUBLIC, PARAMETER ::   lk_dtasal = .TRUE.    !: salinity data flag 
     26   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) ::   &  !: 
     27      s_dta       !: salinity data at given time-step 
    2828 
    2929   !! * Module variables 
     
    8181      INTEGER ::   ji, jj, jk, jl    ! dummy loop indicies 
    8282      INTEGER ::   & 
    83          imois, iman, ik, i15,  &  ! temporary integers 
    84          ipi, ipj, ipk, itime      !    "          " 
     83         imois, iman, ik, i15,       &  ! temporary integers 
     84         ipi, ipj, ipk, itime,       &  !    "          " 
     85         il0, il1, ii0, ii1, ij0, ij1   !    "          " 
    8586      INTEGER, DIMENSION(jpmois) ::   istep 
    8687      REAL(wp) ::   & 
     
    186187            !                                        !  ORCA_R2 configuration 
    187188            !                                        ! ======================= 
    188           
    189             DO jj = mj0(101), mj1(109)                      ! Reduced salinity in the Alboran Sea 
    190                DO ji = mi0(141), mi1(155) 
     189            ij0 = 101   ;   ij1 = 109 
     190            ii0 = 141   ;   ii1 = 155    
     191            DO jj = mj0(ij0), mj1(ij1)                      ! Reduced salinity in the Alboran Sea 
     192               DO ji = mi0(ii0), mi1(ii1) 
    191193                  DO jk = 13, 13 
    192194                     saldta(ji,jj,jk,:) = saldta(ji,jj,jk,:) - 0.15 
     
    205207            IF( n_cla == 1 ) THEN  
    206208               !                                         ! New salinity profile at Gibraltar 
    207                saldta( mi0(139):mi1(139) , mj0(101):mj1(101) , : , : ) =   & 
    208                   &                                    saldta( mi0(138):mi1(138) , mj0(101):mj1(101) , : , : ) 
    209                saldta( mi0(139):mi1(139) , mj0(102):mj1(102) , : , : ) =   & 
    210                   &                                    saldta( mi0(138):mi1(138) , mj0(102):mj1(102) , : , : ) 
    211                DO jl = mi0(138), mi1(138)                ! New temperature profile at Gibraltar 
    212                   DO jj = mj0(101), mj1(102) 
    213                      DO ji = mi0(139), mi1(139) 
     209               il0 = 138   ;   il1 = 138    
     210               ij0 = 101   ;   ij1 = 101 
     211               ii0 = 139   ;   ii1 = 139    
     212               saldta( mi0(ii0):mi1(ii1), mj0(ij0):mj1(ij1) , : , : ) =   & 
     213                  &                                    saldta( mi0(il0):mi1(il1) , mj0(ij0):mj1(ij1) , : , : ) 
     214               ij0 = 101   ;   ij1 = 101 
     215               saldta( mi0(ii0):mi1(ii1), mj0(ij0):mj1(ij1) , : , : ) =   & 
     216                  &                                    saldta( mi0(il0):mi1(il1) , mj0(ij0):mj1(ij1) , : , : ) 
     217               il0 = 138   ;   il1 = 138    
     218               ij0 = 101   ;   ij1 = 102 
     219               ii0 = 139   ;   ii1 = 139    
     220               DO jl = mi0(ii0), mi1(ii1)                ! New salinity profile at Gibraltar 
     221                  DO jj = mj0(ij0), mj1(ij1) 
     222                     DO ji = mi0(ii0), mi1(ii1) 
    214223                        saldta(ji,jj,:,:) = saldta(jl,jj,:,:) 
    215224                     END DO 
     
    217226               END DO 
    218227 
    219                DO jl = mi0(164), mi1(164)                ! New salinity profile at Bab el Mandeb 
    220                   DO jj = mj0(88), mj1(88) 
    221                      DO ji = mi0(161), mi1(163) 
     228               il0 = 164   ;   il1 = 164    
     229               ij0 =  88   ;   ij1 =  88 
     230               ii0 = 161   ;   ii1 = 163    
     231               DO jl = mi0(ii0), mi1(ii1)                ! New salinity profile at Bab el Mandeb 
     232                  DO jj = mj0(ij0), mj1(ij1) 
     233                     DO ji = mi0(ii0), mi1(ii1) 
    222234                        saldta(ji,jj,:,:) = saldta(jl,jj,:,:) 
    223235                     END DO 
    224236                  END DO 
    225                   DO jj = mj0(87), mj1(87) 
    226                      DO ji = mi0(161), mi1(163) 
     237                  ij0 =  87   ;   ij1 =  87 
     238                  DO jj = mj0(ij0), mj1(ij1) 
     239                     DO ji = mi0(ii0), mi1(ii1) 
    227240                        saldta(ji,jj,:,:) = saldta(jl,jj,:,:) 
    228241                     END DO 
     
    280293   !!   Default option:                                    NO salinity data 
    281294   !!---------------------------------------------------------------------- 
    282    LOGICAL, PUBLIC, PARAMETER ::   lk_dtasal = .FALSE.   ! salinity data flag 
     295   LOGICAL , PUBLIC, PARAMETER ::   lk_dtasal = .FALSE.   !: salinity data flag 
    283296CONTAINS 
    284297   SUBROUTINE dta_sal( kt )        ! Empty routine 
    285       WRITE(*,*) kt 
     298      WRITE(*,*) 'dta_sal: You should not have seen this print! error?', kt 
    286299   END SUBROUTINE dta_sal 
    287300#endif 
  • trunk/NEMO/OPA_SRC/DTA/dtasst.F90

    r3 r16  
    2626   !! * Shared module variables 
    2727#if defined key_dtasst 
    28    LOGICAL , PUBLIC, PARAMETER ::   lk_dtasst = .TRUE.   ! sst data flag 
     28   LOGICAL , PUBLIC, PARAMETER ::   lk_dtasst = .TRUE.   !: sst data flag 
    2929#else 
    30    LOGICAL , PUBLIC, PARAMETER ::   lk_dtasst = .TRUE.   ! sst data flag 
     30   LOGICAL , PUBLIC, PARAMETER ::   lk_dtasst = .FALSE.  !: sst data flag 
    3131#endif 
    32    REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   & 
    33       sst             ! surface temperature 
    34    REAL(wp), PUBLIC, DIMENSION(jpi,jpj,2) ::   & 
    35       rclice          ! climatological ice index (0/1) (2 months) 
     32   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   &  !: 
     33      sst             !: surface temperature 
     34   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,2) ::   &  !: 
     35      rclice          !: climatological ice index (0/1) (2 months) 
    3636   !!---------------------------------------------------------------------- 
    3737   !!   OPA 9.0 , IPSL-LODYC  (2003) 
  • trunk/NEMO/OPA_SRC/DTA/dtatem.F90

    r3 r16  
    2323 
    2424   !! * Shared module variables 
    25    LOGICAL , PUBLIC, PARAMETER ::   lk_dtatem = .TRUE.   ! temperature data flag 
    26    REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) ::   & 
    27       t_dta             ! temperature data at given time-step 
     25   LOGICAL , PUBLIC, PARAMETER ::   lk_dtatem = .TRUE.   !: temperature data flag 
     26   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) ::   &  !: 
     27      t_dta             !: temperature data at given time-step 
    2828 
    2929   !! * Module variables 
     
    9191      INTEGER ::   & 
    9292         imois, iman, itime, ik ,    &  ! temporary integers 
    93          i15, ipi, ipj, ipk             !    "          " 
     93         i15, ipi, ipj, ipk,         &  !    "          " 
     94         il0, il1, ii0, ii1, ij0, ij1   !    "          " 
    9495 
    9596      INTEGER, DIMENSION(jpmois) ::   istep 
     
    192193            !                                        ! =======================  
    193194 
    194             DO jj = mj0(101), mj1(109)                   ! Reduced temperature at Alboran Sea 
    195                DO ji = mi0(141), mi1(155) 
     195            ij0 = 101   ;   ij1 = 109 
     196            ii0 = 141   ;   ii1 = 155 
     197            DO jj = mj0(ij0), mj1(ij1)                      ! Reduced temperature in the Alboran Sea 
     198               DO ji = mi0(ii0), mi1(ii1) 
    196199                  temdta(ji,jj, 13:13 ,:) = temdta(ji,jj, 13:13 ,:) - 0.20 
    197200                  temdta(ji,jj, 14:15 ,:) = temdta(ji,jj, 14:15 ,:) - 0.35 
     
    202205            IF( n_cla == 0 ) THEN  
    203206               !                                         ! Reduced temperature at Red Sea 
    204                temdta( mi0(148):mi1(160) , mj0(87):mj1(96) ,  4:10 , : ) = 7.0  
    205                temdta( mi0(148):mi1(160) , mj0(87):mj1(96) , 11:13 , : ) = 6.5  
    206                temdta( mi0(148):mi1(160) , mj0(87):mj1(96) , 14:20 , : ) = 6.0 
     207               ij0 =  87   ;   ij1 =  96 
     208               ii0 = 148   ;   ii1 = 160 
     209               temdta( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ,  4:10 , : ) = 7.0  
     210               temdta( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 11:13 , : ) = 6.5  
     211               temdta( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 14:20 , : ) = 6.0 
    207212            ELSE 
    208                DO jl = mi0(138), mi1(138)                ! New temperature profile at Gibraltar 
    209                   DO jj = mj0(101), mj1(102) 
    210                      DO ji = mi0(139), mi1(139) 
     213               il0 = 138   ;   il1 = 138 
     214               ij0 = 101   ;   ij1 = 102 
     215               ii0 = 139   ;   ii1 = 139 
     216               DO jl = mi0(ii0), mi1(ii1)                ! New temperature profile at Gibraltar 
     217                  DO jj = mj0(ij0), mj1(ij1) 
     218                     DO ji = mi0(ii0), mi1(ii1) 
    211219                        temdta(ji,jj,:,:) = temdta(jl,jj,:,:) 
    212220                     END DO 
    213221                  END DO 
    214222               END DO 
    215                DO jl = mi0(164), mi1(164)                ! New temperature profile at Bab el Mandeb 
    216                   DO jj = mj0(88), mj1(88) 
    217                      DO ji = mi0(161), mi1(163) 
     223               il0 = 164   ;   il1 = 164 
     224               ij0 =  88   ;   ij1 =  88 
     225               ii0 = 161   ;   ii1 = 163 
     226               DO jl = mi0(ii0), mi1(ii1)                ! New temperature profile at Bab el Mandeb 
     227                  DO jj = mj0(ij0), mj1(ij1) 
     228                     DO ji = mi0(ii0), mi1(ii1) 
    218229                        temdta(ji,jj,:,:) = temdta(jl,jj,:,:) 
    219230                     END DO 
    220231                  END DO 
    221                   DO jj = mj0(87), mj1(87) 
    222                      DO ji = mi0(161), mi1(163) 
     232                  ij0 =  87   ;   ij1 =  87 
     233                  DO jj = mj0(ij0), mj1(ij1) 
     234                     DO ji = mi0(ii0), mi1(ii1) 
    223235                        temdta(ji,jj,:,:) = temdta(jl,jj,:,:) 
    224236                     END DO 
     
    274286   !!   Default case                           NO 3D temperature data field 
    275287   !!---------------------------------------------------------------------- 
    276    LOGICAL , PUBLIC, PARAMETER ::   lk_dtatem = .FALSE.   ! temperature data flag 
     288   LOGICAL , PUBLIC, PARAMETER ::   lk_dtatem = .FALSE.   !: temperature data flag 
    277289CONTAINS 
    278290   SUBROUTINE dta_tem( kt )        ! Empty routine 
    279       WRITE(*,*) kt 
     291      WRITE(*,*) 'dta_tem: You should not have seen this print! error?', kt 
    280292   END SUBROUTINE dta_tem 
    281293#endif 
Note: See TracChangeset for help on using the changeset viewer.