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 497 for trunk/NEMO/OFF_SRC/DOM/domrea_dimg.h90 – NEMO

Ignore:
Timestamp:
2006-09-12T13:03:53+02:00 (18 years ago)
Author:
opalod
Message:

nemo_v1_update_064:CE:re-organization of coordinate definition and scale factors

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/NEMO/OFF_SRC/DOM/domrea_dimg.h90

    r343 r497  
    7171       irecl8 = jpi*jpj*wp 
    7272       ios1 = 0 ; ios2 = 0 
    73        IF (lk_zps ) ios1= 1  
    74        IF (lk_sco ) ios2= 1 
     73       IF (ln_zps ) ios1= 1  
     74       IF (ln_sco ) ios2= 1 
    7575    
    7676       SELECT CASE (nmsh ) 
    … …  
    262262          mbathy=zprt*tmask(:,:,1)+1 
    263263 
    264 #if defined key_s_coord 
    265          ! 
    266          ! hbat 
    267          irec = irecv(inum4) + (narea - 1 ) 
    268          READ(inum4,REC=irec) hbatt(:,:) 
    269          irecv(inum4) = irecv(inum4) + jpnij  
    270         
    271          irec = irecv(inum4) + (narea - 1 ) 
    272          READ(inum4,REC=irec) hbatu(:,:) 
    273          irecv(inum4) = irecv(inum4) + jpnij  
    274  
    275          irec = irecv(inum4) + (narea - 1 ) 
    276          READ(inum4,REC=irec) hbatv(:,:) 
    277          irecv(inum4) = irecv(inum4) + jpnij  
    278  
    279          irec = irecv(inum4) + (narea - 1 ) 
    280          READ(inum4,REC=irec) hbatf(:,:) 
    281          irecv(inum4) = irecv(inum4) + jpnij  
    282          ! 
    283          ! gsig and esig ( as vectors of jpk per record ) 
    284          irec =  irecv(inum4) + (narea - 1 ) 
    285          READ(inum4,REC=irec) gsigt(:) 
    286          irecv(inum4) = irecv(inum4) + jpnij  
    287  
    288          irec =  irecv(inum4) + (narea - 1 ) 
    289          READ(inum4,REC=irec) gsigw(:) 
    290          irecv(inum4) = irecv(inum4) + jpnij  
    291  
    292          irec =  irecv(inum4) + (narea - 1 ) 
    293          READ(inum4,REC=irec) gsi3w(:) 
    294          irecv(inum4) = irecv(inum4) + jpnij  
    295  
    296          irec =  irecv(inum4) + (narea - 1 ) 
    297          READ(inum4,REC=irec) esigt(:) 
    298          irecv(inum4) = irecv(inum4) + jpnij  
    299  
    300          irec =  irecv(inum4) + (narea - 1 ) 
    301          READ(inum4,REC=irec) esigw(:) 
    302          irecv(inum4) = irecv(inum4) + jpnij  
    303  
    304 # elif defined key_partial_steps 
    305          ! 
    306          ! hdep 
    307          irec = irecv(inum4) + (narea - 1 ) 
    308          READ(inum4,REC=irec) hdept(:,:) 
    309          irecv(inum4) = irecv(inum4) + jpnij  
    310         
    311          irec = irecv(inum4) + (narea - 1 ) 
    312          READ(inum4,REC=irec) hdepw(:,:) 
    313          irecv(inum4) = irecv(inum4) + jpnij  
    314          ! 
    315          ! e3t_ps (3D) 
    316          DO jk=1,jpk 
    317              irec=irecv(inum4) + (jk -1 ) + jpk * (narea - 1 ) 
    318              READ(inum4,REC=irec) e3t_ps(:,:,jk) 
    319          END DO 
    320          irecv(inum4) = irecv(inum4) + jpk * jpnij  
    321  
    322          ! e3u_ps e3v_ps e3w_ps (3D) 
    323          DO jk=1,jpk 
    324              irec=irecv(inum4) + (jk -1 ) + jpk * (narea - 1 ) 
    325              READ(inum4,REC=irec) e3u_ps(:,:,jk) 
    326          END DO 
    327          irecv(inum4) = irecv(inum4) + jpk * jpnij  
    328           
    329          DO jk=1,jpk 
    330              irec=irecv(inum4) + (jk -1 ) + jpk * (narea - 1 ) 
    331              READ(inum4,REC=irec) e3v_ps(:,:,jk) 
    332          END DO 
    333          irecv(inum4) = irecv(inum4) + jpk * jpnij  
    334  
    335          DO jk=1,jpk 
    336              irec=irecv(inum4) + (jk -1 ) + jpk * (narea - 1 ) 
    337              READ(inum4,REC=irec) e3w_ps(:,:,jk) 
    338          END DO 
    339          irecv(inum4) = irecv(inum4) + jpk * jpnij  
    340          !  
    341          !  
    342          ! gdep 
    343          irec =  irecv(inum4) + (narea - 1 ) 
    344          READ(inum4,REC=irec) gdept(:) 
    345          irecv(inum4) = irecv(inum4) + jpnij  
    346  
    347          irec =  irecv(inum4) + (narea - 1 ) 
    348          READ(inum4,REC=irec) gdepw(:) 
    349          irecv(inum4) = irecv(inum4) + jpnij  
    350          ! 
    351          ! e3 
    352          irec =  irecv(inum4) + (narea - 1 ) 
    353          READ(inum4,REC=irec) e3t(:) 
    354          irecv(inum4) = irecv(inum4) + jpnij  
    355  
    356          irec =  irecv(inum4) + (narea - 1 ) 
    357          READ(inum4,REC=irec) e3w(:) 
     264#if ! defined key_zco 
     265          IF( ln_sco ) THEN                                         ! s-coordinate 
     266             ! 
     267             ! hbat 
     268             irec = irecv(inum4) + (narea - 1 ) 
     269             READ(inum4,REC=irec) hbatt(:,:) 
     270             irecv(inum4) = irecv(inum4) + jpnij  
     271              
     272             irec = irecv(inum4) + (narea - 1 ) 
     273             READ(inum4,REC=irec) hbatu(:,:) 
     274             irecv(inum4) = irecv(inum4) + jpnij  
     275              
     276             irec = irecv(inum4) + (narea - 1 ) 
     277             READ(inum4,REC=irec) hbatv(:,:) 
     278             irecv(inum4) = irecv(inum4) + jpnij  
     279              
     280             irec = irecv(inum4) + (narea - 1 ) 
     281             READ(inum4,REC=irec) hbatf(:,:) 
     282             irecv(inum4) = irecv(inum4) + jpnij  
     283             ! 
     284             ! gsig and esig ( as vectors of jpk per record ) 
     285             irec =  irecv(inum4) + (narea - 1 ) 
     286             READ(inum4,REC=irec) gsigt(:) 
     287             irecv(inum4) = irecv(inum4) + jpnij  
     288              
     289             irec =  irecv(inum4) + (narea - 1 ) 
     290             READ(inum4,REC=irec) gsigw(:) 
     291             irecv(inum4) = irecv(inum4) + jpnij  
     292              
     293             irec =  irecv(inum4) + (narea - 1 ) 
     294             READ(inum4,REC=irec) gsi3w(:) 
     295             irecv(inum4) = irecv(inum4) + jpnij  
     296              
     297             irec =  irecv(inum4) + (narea - 1 ) 
     298             READ(inum4,REC=irec) esigt(:) 
     299             irecv(inum4) = irecv(inum4) + jpnij  
     300              
     301             irec =  irecv(inum4) + (narea - 1 ) 
     302             READ(inum4,REC=irec) esigw(:) 
     303             irecv(inum4) = irecv(inum4) + jpnij  
     304 
     305             ! e3 (3D) 
     306             DO jk=1,jpk 
     307                irec =  irecv(inum4) + (jk -1 ) + jpk * (narea - 1 ) 
     308                READ(inum4,REC=irec) e3t(:,:,jk) 
     309             ENDIF 
     310             irecv(inum4) = irecv(inum4) + jpk * jpnij  
     311              
     312             DO jk=1,jpk 
     313                irec =  irecv(inum4) + (jk -1 ) + jpk * (narea - 1 ) 
     314                READ(inum4,REC=irec) e3u(:,:,jk) 
     315             ENDIF 
     316             irecv(inum4) = irecv(inum4) + jpk * jpnij  
     317              
     318             DO jk=1,jpk 
     319                irec =  irecv(inum4) + (jk -1 ) + jpk * (narea - 1 ) 
     320                READ(inum4,REC=irec) e3v(:,:,jk) 
     321             ENDIF 
     322             irecv(inum4) = irecv(inum4) + jpk * jpnij  
     323              
     324             DO jk=1,jpk 
     325                irec =  irecv(inum4) + (jk -1 ) + jpk * (narea - 1 ) 
     326                READ(inum4,REC=irec) e3w(:,:,jk) 
     327             ENDIF 
     328             irecv(inum4) = irecv(inum4) + jpk * jpnij  
     329 
     330             ! gdep 
     331             irec =  irecv(inum4) + (narea - 1 ) 
     332             READ(inum4,REC=irec) gdept_0(:) 
     333             irecv(inum4) = irecv(inum4) + jpnij  
     334              
     335             irec =  irecv(inum4) + (narea - 1 ) 
     336             READ(inum4,REC=irec) gdepw_0(:) 
     337             irecv(inum4) = irecv(inum4) + jpnij  
     338             ! 
     339          ENDIF 
     340              
     341          IF( ln_zps ) THEN 
     342             ! 
     343             ! hdep 
     344             irec = irecv(inum4) + (narea - 1 ) 
     345             READ(inum4,REC=irec) hdept(:,:) 
     346             irecv(inum4) = irecv(inum4) + jpnij  
     347              
     348             irec = irecv(inum4) + (narea - 1 ) 
     349             READ(inum4,REC=irec) hdepw(:,:) 
     350             irecv(inum4) = irecv(inum4) + jpnij  
     351             ! 
     352             ! e3t_ps (3D) 
     353             DO jk=1,jpk 
     354                irec=irecv(inum4) + (jk -1 ) + jpk * (narea - 1 ) 
     355                READ(inum4,REC=irec) e3t(:,:,jk) 
     356             END DO 
     357             irecv(inum4) = irecv(inum4) + jpk * jpnij  
     358              
     359             ! e3u_ps e3v_ps e3w_ps (3D) 
     360             DO jk=1,jpk 
     361                irec=irecv(inum4) + (jk -1 ) + jpk * (narea - 1 ) 
     362                READ(inum4,REC=irec) e3u(:,:,jk) 
     363             END DO 
     364             irecv(inum4) = irecv(inum4) + jpk * jpnij  
     365              
     366             DO jk=1,jpk 
     367                irec=irecv(inum4) + (jk -1 ) + jpk * (narea - 1 ) 
     368                READ(inum4,REC=irec) e3v(:,:,jk) 
     369             END DO 
     370             irecv(inum4) = irecv(inum4) + jpk * jpnij  
     371              
     372             DO jk=1,jpk 
     373                irec=irecv(inum4) + (jk -1 ) + jpk * (narea - 1 ) 
     374                READ(inum4,REC=irec) e3w(:,:,jk) 
     375             END DO 
     376             irecv(inum4) = irecv(inum4) + jpk * jpnij  
     377             !  
     378             !  
     379             ! gdep 
     380             irec =  irecv(inum4) + (narea - 1 ) 
     381             READ(inum4,REC=irec) gdept_0(:) 
     382             irecv(inum4) = irecv(inum4) + jpnij  
     383              
     384             irec =  irecv(inum4) + (narea - 1 ) 
     385             READ(inum4,REC=irec) gdepw_0(:) 
     386             irecv(inum4) = irecv(inum4) + jpnij  
     387             ! 
     388             ! e3 
     389             irec =  irecv(inum4) + (narea - 1 ) 
     390             READ(inum4,REC=irec) e3t_0(:) 
     391             irecv(inum4) = irecv(inum4) + jpnij  
     392              
     393             irec =  irecv(inum4) + (narea - 1 ) 
     394             READ(inum4,REC=irec) e3w_0(:) 
     395          ENDIF 
    358396#else 
    359          ! 
    360          ! gdep 
    361          irec =  irecv(inum4) + (narea - 1 ) 
    362          READ(inum4,REC=irec) gdept(:) 
    363          irecv(inum4) = irecv(inum4) + jpnij  
    364  
    365          irec =  irecv(inum4) + (narea - 1 ) 
    366          READ(inum4,REC=irec) gdepw(:) 
    367          irecv(inum4) = irecv(inum4) + jpnij  
    368          ! 
    369          ! e3 
    370          irec =  irecv(inum4) + (narea - 1 ) 
    371          READ(inum4,REC=irec) e3t(:) 
    372          irecv(inum4) = irecv(inum4) + jpnij  
    373  
    374          irec =  irecv(inum4) + (narea - 1 ) 
    375          READ(inum4,REC=irec) e3w(:) 
    376          irecv(inum4) = irecv(inum4) + jpnij  
     397          ! 
     398          ! gdep 
     399          irec =  irecv(inum4) + (narea - 1 ) 
     400          READ(inum4,REC=irec) gdept_0(:) 
     401          irecv(inum4) = irecv(inum4) + jpnij  
     402           
     403          irec =  irecv(inum4) + (narea - 1 ) 
     404          READ(inum4,REC=irec) gdepw_0(:) 
     405          irecv(inum4) = irecv(inum4) + jpnij  
     406          ! 
     407          ! e3 
     408          irec =  irecv(inum4) + (narea - 1 ) 
     409          READ(inum4,REC=irec) e3t_0(:) 
     410          irecv(inum4) = irecv(inum4) + jpnij  
     411           
     412          irec =  irecv(inum4) + (narea - 1 ) 
     413          READ(inum4,REC=irec) e3w_0(:) 
     414          irecv(inum4) = irecv(inum4) + jpnij  
    377415         ! 
    378416#endif 
Note: See TracChangeset for help on using the changeset viewer.