Changeset 358


Ignore:
Timestamp:
09/07/15 18:28:16 (9 years ago)
Author:
llfita
Message:

Modifying code to be 'gfortran' compatible:

  • Change WRITE(' ',)... by PRINT *,
  • Shorten some code lines
Location:
codes/icosagcm/trunk/src
Files:
5 edited

Legend:

Unmodified
Added
Removed

  • codes/icosagcm/trunk/src/dissip_gcm.f90

    r295 r358  
    111111      IF (is_master) PRINT *, 'Rayleigh friction : Schaer-like mountain with shear DCMIP2.2' 
    112112   CASE DEFAULT 
    113       IF (is_master) PRINT *, 'Bad selector : rayleigh_friction_type =', TRIM(rayleigh_friction_key), ' in dissip_gcm.f90/init_dissip' 
     113      IF (is_master) PRINT *, 'Bad selector : rayleigh_friction_type =', TRIM(rayleigh_friction_key), & 
     114        ' in dissip_gcm.f90/init_dissip' 
    114115      STOP 
    115116   END SELECT 

  • codes/icosagcm/trunk/src/domain.f90

    r327 r358  
    472472    
    473473      ind_glo=0 
    474       WRITE(*,'') 
     474      PRINT *,'' 
    475475      PRINT*, '      MPI PROCESS DISTRIBUTION' 
    476       WRITE(*,'') 
     476      PRINT *,'' 
    477477       
    478478      WRITE(*,"(' ')", ADVANCE='NO') 
    … …  
    480480        WRITE(*,"('=')", ADVANCE='NO') 
    481481      ENDDO 
    482       WRITE(*,'') 
     482      PRINT *,'' 
    483483 
    484484      DO nf=1,nb_face 
    … …  
    489489              WRITE(*,"('-')", ADVANCE='NO') 
    490490            ENDDO 
    491             WRITE(*,'') 
     491            PRINT *,'' 
    492492          ENDIF 
    493493 
    … …  
    496496            WRITE(*,"(' ','    ',' |')",ADVANCE='NO')          
    497497          ENDDO 
    498           WRITE(*,'') 
     498          PRINT *,'' 
    499499 
    500500          WRITE(*,"('|')", ADVANCE='NO') 
    … …  
    503503            WRITE(*,"(' ',i4.4  ,' |')",ADVANCE='NO'),domglo_rank(ind_glo)            
    504504          END DO 
    505           WRITE(*,'') 
     505          PRINT *,'' 
    506506 
    507507          WRITE(*,"('|')", ADVANCE='NO') 
    … …  
    509509            WRITE(*,"(' ','    ',' |')",ADVANCE='NO')          
    510510          ENDDO 
    511           WRITE(*,'') 
     511          PRINT *,'' 
    512512 
    513513        ENDDO 
    … …  
    517517          WRITE(*,"('=')", ADVANCE='NO') 
    518518        ENDDO 
    519         WRITE(*,'') 
     519        PRINT *,'' 
    520520      ENDDO 
    521521    ENDIF 

  • codes/icosagcm/trunk/src/restart.f90

    r327 r358  
    269269                DO i=d%ii_begin,d%ii_end 
    270270                  DO k=0,5 
    271                    IF (d%edge_assign_domain(k,i,j)==ind .AND. d%edge_assign_i(k,i,j)==i .AND. d%edge_assign_j(k,i,j)==j .AND. d%edge_assign_pos(k,i,j)==k) THEN 
     271                   IF (d%edge_assign_domain(k,i,j)==ind .AND. d%edge_assign_i(k,i,j)==i .AND. d%edge_assign_j(k,i,j)==j & 
     272                     .AND. d%edge_assign_pos(k,i,j)==k) THEN 
    272273                       ij=(j-1)*d%iim+i 
    273274                      ind_glo=d%assign_cell_glo(i,j) 
    … …  
    287288                DO i=d%ii_begin,d%ii_end 
    288289                  DO k=0,5 
    289                    IF (d%edge_assign_domain(k,i,j)==ind .AND. d%edge_assign_i(k,i,j)==i .AND. d%edge_assign_j(k,i,j)==j .AND. d%edge_assign_pos(k,i,j)==k) THEN 
     290                   IF (d%edge_assign_domain(k,i,j)==ind .AND. d%edge_assign_i(k,i,j)==i .AND. d%edge_assign_j(k,i,j)==j & 
     291                     .AND. d%edge_assign_pos(k,i,j)==k) THEN 
    290292                       ij=(j-1)*d%iim+i 
    291293                      ind_glo=d%assign_cell_glo(i,j) 
    … …  
    305307                DO i=d%ii_begin,d%ii_end 
    306308                  DO k=0,5 
    307                     IF (d%edge_assign_domain(k,i,j)==ind .AND. d%edge_assign_i(k,i,j)==i .AND. d%edge_assign_j(k,i,j)==j .AND. d%edge_assign_pos(k,i,j)==k) THEN 
     309                    IF (d%edge_assign_domain(k,i,j)==ind .AND. d%edge_assign_i(k,i,j)==i .AND. d%edge_assign_j(k,i,j)==j& 
     310                      .AND. d%edge_assign_pos(k,i,j)==k) THEN 
    308311                      ij=(j-1)*d%iim+i 
    309312                      ind_glo=d%assign_cell_glo(i,j) 

  • codes/icosagcm/trunk/src/transfert.F90

    r266 r358  
    33#ifdef CPP_USING_MPI 
    44  USE transfert_mpi_mod, ONLY : init_transfert, transfert_request=>transfert_request_mpi, req_i1,req_e1_vect, & 
    5                                 req_e1_scal, req_i0, req_e0_vect, req_e0_scal, request_add_point, create_request, gather_field, scatter_field, & 
    6                                 t_message,init_message=>init_message_mpi,transfert_message=>transfert_message_mpi,  & 
    7                                 send_message=>send_message_mpi,test_message=>test_message_mpi,wait_message=>wait_message_mpi,barrier, & 
    8                                 bcast_mpi 
     5                                req_e1_scal, req_i0, req_e0_vect, req_e0_scal, request_add_point,             & 
     6                                create_request, gather_field, scatter_field,         & 
     7                                t_message, init_message=>init_message_mpi,           & 
     8                                transfert_message=>transfert_message_mpi,            & 
     9                                send_message=>send_message_mpi,                      & 
     10                                test_message=>test_message_mpi,                      & 
     11                                wait_message=>wait_message_mpi,barrier,bcast_mpi 
    912#else  
    1013  USE transfert_mpi_mod, ONLY : init_transfert, transfert_request=>transfert_request_seq, req_i1,req_e1_vect, & 
    11                                 req_e1_scal,req_i0, req_e0_vect, req_e0_scal, request_add_point, create_request, gather_field, scatter_field,& 
    12                                 t_message,init_message=>init_message_seq,transfert_message=>transfert_message_seq,  & 
    13                                 send_message=>send_message_seq,test_message=>test_message_seq,wait_message=>wait_message_seq,barrier, & 
    14                                 bcast_mpi 
     14                                req_e1_scal,req_i0, req_e0_vect, req_e0_scal,        & 
     15                                request_add_point, create_request, gather_field,     & 
     16                                scatter_field, t_message,                            & 
     17                                init_message=>init_message_seq,                      & 
     18                                transfert_message=>transfert_message_seq,            & 
     19                                send_message=>send_message_seq,                      & 
     20                                test_message=>test_message_seq,                      & 
     21                                wait_message=>wait_message_seq,barrier, bcas_mpi 
    1522#endif 
    1623 

  • codes/icosagcm/trunk/src/transfert_mpi.f90

    r327 r358  
    11501150              IF (mpi_threading_mode==MPI_THREAD_SERIALIZED) THEN 
    11511151                !$OMP CRITICAL             
    1152                 CALL MPI_ISEND(buffer_r,send%size,MPI_REAL8,send%rank,send%tag+1000*message%number,comm_icosa, message%mpi_req(ireq),ierr) 
     1152                CALL MPI_ISEND(buffer_r,send%size,MPI_REAL8,send%rank,               & 
     1153                  send%tag+1000*message%number,comm_icosa, message%mpi_req(ireq),ierr) 
    11531154                !$OMP END CRITICAL 
    11541155              ELSE IF (mpi_threading_mode==MPI_THREAD_MULTIPLE) THEN 
    1155                 CALL MPI_ISEND(buffer_r,send%size,MPI_REAL8,send%rank,send%tag+1000*message%number,comm_icosa, message%mpi_req(ireq),ierr) 
     1156                CALL MPI_ISEND(buffer_r,send%size,MPI_REAL8,send%rank,               & 
     1157                  send%tag+1000*message%number,comm_icosa, message%mpi_req(ireq),ierr) 
    11561158              ENDIF 
    11571159              
    … …  
    11861188              IF (mpi_threading_mode==MPI_THREAD_SERIALIZED) THEN 
    11871189               !$OMP CRITICAL             
    1188                 CALL MPI_IRECV(buffer_r,recv%size,MPI_REAL8,recv%rank,recv%tag+1000*message%number,comm_icosa, message%mpi_req(ireq),ierr) 
     1190                CALL MPI_IRECV(buffer_r,recv%size,MPI_REAL8,recv%rank,               & 
     1191                  recv%tag+1000*message%number,comm_icosa, message%mpi_req(ireq),ierr) 
    11891192               !$OMP END CRITICAL 
    11901193              ELSE IF (mpi_threading_mode==MPI_THREAD_MULTIPLE) THEN 
    1191                  CALL MPI_IRECV(buffer_r,recv%size,MPI_REAL8,recv%rank,recv%tag+1000*message%number,comm_icosa, message%mpi_req(ireq),ierr) 
     1194                 CALL MPI_IRECV(buffer_r,recv%size,MPI_REAL8,recv%rank,              & 
     1195                   recv%tag+1000*message%number,comm_icosa, message%mpi_req(ireq),ierr) 
    11921196              ENDIF 
    11931197             
    … …  
    12321236                IF (mpi_threading_mode==MPI_THREAD_SERIALIZED) THEN 
    12331237                  !$OMP CRITICAL    
    1234                   CALL MPI_ISEND(buffer_r,send%size*dim3,MPI_REAL8,send%rank,send%tag+1000*message%number,comm_icosa, message%mpi_req(ireq),ierr) 
     1238                  CALL MPI_ISEND(buffer_r,send%size*dim3,MPI_REAL8,send%rank,        & 
     1239                    send%tag+1000*message%number,comm_icosa, message%mpi_req(ireq),ierr) 
    12351240                  !$OMP END CRITICAL 
    12361241                ELSE IF (mpi_threading_mode==MPI_THREAD_MULTIPLE) THEN 
    1237                   CALL MPI_ISEND(buffer_r,send%size*dim3,MPI_REAL8,send%rank,send%tag+1000*message%number,comm_icosa, message%mpi_req(ireq),ierr) 
     1242                  CALL MPI_ISEND(buffer_r,send%size*dim3,MPI_REAL8,send%rank,        & 
     1243                    send%tag+1000*message%number,comm_icosa, message%mpi_req(ireq),ierr) 
    12381244                ENDIF 
    12391245              ENDIF 
    … …  
    12741280                IF (mpi_threading_mode==MPI_THREAD_SERIALIZED) THEN 
    12751281                  !$OMP CRITICAL 
    1276                   CALL MPI_IRECV(buffer_r,recv%size*dim3,MPI_REAL8,recv%rank,recv%tag+1000*message%number,comm_icosa, message%mpi_req(ireq),ierr) 
     1282                  CALL MPI_IRECV(buffer_r,recv%size*dim3,MPI_REAL8,recv%rank,        & 
     1283                    recv%tag+1000*message%number,comm_icosa, message%mpi_req(ireq),ierr) 
    12771284                  !$OMP END CRITICAL 
    12781285                ELSE IF (mpi_threading_mode==MPI_THREAD_MULTIPLE) THEN 
    1279                   CALL MPI_IRECV(buffer_r,recv%size*dim3,MPI_REAL8,recv%rank,recv%tag+1000*message%number,comm_icosa, message%mpi_req(ireq),ierr) 
     1286                  CALL MPI_IRECV(buffer_r,recv%size*dim3,MPI_REAL8,recv%rank,        & 
     1287                    recv%tag+1000*message%number,comm_icosa, message%mpi_req(ireq),ierr) 
    12801288                ENDIF 
    12811289              ENDIF 
    … …  
    13071315              CALL trace_start("copy_to_buffer") 
    13081316              DO d4=1,dim4 
    1309                 offset=send%offset*dim3*dim4 + dim3*send%size*(d4-1) + (lbegin-1)*send%size 
     1317                offset=send%offset*dim3*dim4 + dim3*send%size*(d4-1) +               & 
     1318                  (lbegin-1)*send%size 
    13101319 
    13111320                DO d3=lbegin,lend 
    … …  
    13211330                IF (mpi_threading_mode==MPI_THREAD_SERIALIZED) THEN 
    13221331                  !$OMP CRITICAL 
    1323                   CALL MPI_ISEND(buffer_r,send%size*dim3*dim4,MPI_REAL8,send%rank,send%tag+1000*message%number,comm_icosa, message%mpi_req(ireq),ierr) 
     1332                  CALL MPI_ISEND(buffer_r,send%size*dim3*dim4,MPI_REAL8,send%rank,   & 
     1333                    send%tag+1000*message%number,comm_icosa, message%mpi_req(ireq),ierr) 
    13241334                  !$OMP END CRITICAL 
    13251335                ELSE IF (mpi_threading_mode==MPI_THREAD_MULTIPLE) THEN 
    1326                   CALL MPI_ISEND(buffer_r,send%size*dim3*dim4,MPI_REAL8,send%rank,send%tag+1000*message%number,comm_icosa, message%mpi_req(ireq),ierr) 
     1336                  CALL MPI_ISEND(buffer_r,send%size*dim3*dim4,MPI_REAL8,send%rank,   & 
     1337                    send%tag+1000*message%number,comm_icosa, message%mpi_req(ireq),ierr) 
    13271338                ENDIF 
    13281339              ENDIF 
    … …  
    13671378                IF (mpi_threading_mode==MPI_THREAD_SERIALIZED) THEN 
    13681379                 !$OMP CRITICAL            
    1369                   CALL MPI_IRECV(buffer_r,recv%size*dim3*dim4,MPI_REAL8,recv%rank,recv%tag+1000*message%number,comm_icosa, message%mpi_req(ireq),ierr)  
     1380                  CALL MPI_IRECV(buffer_r,recv%size*dim3*dim4,MPI_REAL8,recv%rank,   & 
     1381                    recv%tag+1000*message%number,comm_icosa, message%mpi_req(ireq),ierr)  
    13701382                  !$OMP END CRITICAL 
    13711383                ELSE IF (mpi_threading_mode==MPI_THREAD_MULTIPLE) THEN 
    1372                   CALL MPI_IRECV(buffer_r,recv%size*dim3*dim4,MPI_REAL8,recv%rank,recv%tag+1000*message%number,comm_icosa, message%mpi_req(ireq),ierr)  
     1384                  CALL MPI_IRECV(buffer_r,recv%size*dim3*dim4,MPI_REAL8,recv%rank,   & 
     1385                    recv%tag+1000*message%number,comm_icosa, message%mpi_req(ireq),ierr)  
    13731386                ENDIF 
    13741387              ENDIF 
    … …  
    13871400          msize=message%buffers(ireq)%size 
    13881401          rank=message%buffers(ireq)%rank 
    1389           CALL MPI_ISEND(buffer_r,msize,MPI_REAL8,rank,1000*message%number,comm_icosa, message%mpi_req(ireq),ierr) 
     1402          CALL MPI_ISEND(buffer_r,msize,MPI_REAL8,rank,1000*message%number,          & 
     1403            comm_icosa, message%mpi_req(ireq),ierr) 
    13901404        ENDDO 
    13911405 
    … …  
    13941408          msize=message%buffers(ireq)%size 
    13951409          rank=message%buffers(ireq)%rank 
    1396           CALL MPI_IRECV(buffer_r,msize,MPI_REAL8,rank,1000*message%number,comm_icosa, message%mpi_req(ireq),ierr) 
     1410          CALL MPI_IRECV(buffer_r,msize,MPI_REAL8,rank,1000*message%number,          & 
     1411            comm_icosa, message%mpi_req(ireq),ierr) 
    13971412        ENDDO 
    13981413 
    … …  
    14131428 
    14141429!$OMP MASTER 
    1415     IF (message%pending .AND. .NOT. message%completed) CALL MPI_TESTALL(message%nreq,message%mpi_req,message%completed,message%status,ierr) 
     1430    IF (message%pending .AND. .NOT. message%completed) CALL MPI_TESTALL(message%nreq,& 
     1431      message%mpi_req,message%completed,message%status,ierr) 
    14161432!$OMP END MASTER 
    14171433  END SUBROUTINE  test_message_mpi 
    … …  
    14561472 
    14571473!$OMP MASTER 
    1458         IF (.NOT. message%completed) CALL MPI_WAITALL(nreq,message%mpi_req,message%status,ierr) 
     1474        IF (.NOT. message%completed) CALL MPI_WAITALL(nreq,message%mpi_req,          & 
     1475          message%status,ierr) 
    14591476!$OMP END MASTER 
    14601477!$OMP BARRIER 
    … …  
    14861503 
    14871504!$OMP MASTER 
    1488         IF (.NOT. message%completed) CALL MPI_WAITALL(nreq,message%mpi_req,message%status,ierr) 
     1505        IF (.NOT. message%completed) CALL MPI_WAITALL(nreq,message%mpi_req,          & 
     1506          message%status,ierr) 
    14891507!$OMP END MASTER 
    14901508!$OMP BARRIER 
    … …  
    15341552      ELSE  IF (field(1)%ndim==4) THEN 
    15351553!$OMP MASTER 
    1536         IF (.NOT. message%completed) CALL MPI_WAITALL(nreq,message%mpi_req,message%status,ierr) 
     1554        IF (.NOT. message%completed) CALL MPI_WAITALL(nreq,message%mpi_req,          & 
     1555          message%status,ierr) 
    15371556!$OMP END MASTER 
    15381557!$OMP BARRIER 
    … …  
    15571576              CALL trace_start("copy_from_buffer") 
    15581577              DO d4=1,dim4 
    1559                 offset=recv%offset*dim3*dim4 + dim3*recv%size*(d4-1) + (lbegin-1)*recv%size 
     1578                offset=recv%offset*dim3*dim4 + dim3*recv%size*(d4-1) +               & 
     1579                  (lbegin-1)*recv%size 
    15601580                DO d3=lbegin,lend 
    15611581                  DO n=1,recv%size 
    … …  
    16241644        IF (field(ind)%ndim==2) THEN 
    16251645          DO n=1,req%size 
    1626             rval2d(req%target_ind(n))=field(req%src_domain(n))%rval2d(req%src_ind(n))*req%target_sign(n) 
     1646            rval2d(req%target_ind(n))=field(req%src_domain(n))%rval2d(req%src_ind(n))*& 
     1647              req%target_sign(n) 
    16271648          ENDDO 
    16281649        ELSE IF (field(ind)%ndim==3) THEN 
    16291650          DO n=1,req%size 
    1630             rval3d(req%target_ind(n),:)=field(req%src_domain(n))%rval3d(req%src_ind(n),:)*req%target_sign(n) 
     1651            rval3d(req%target_ind(n),:)=field(req%src_domain(n))%rval3d(req%src_ind(n),:)*& 
     1652              req%target_sign(n) 
    16311653          ENDDO 
    16321654        ELSE IF (field(ind)%ndim==4) THEN 
    16331655          DO n=1,req%size 
    1634             rval4d(req%target_ind(n),:,:)=field(req%src_domain(n))%rval4d(req%src_ind(n),:,:)*req%target_sign(n) 
     1656            rval4d(req%target_ind(n),:,:)=field(req%src_domain(n))%rval4d(req%src_ind(n),:,:)*& 
     1657              req%target_sign(n) 
    16351658          ENDDO 
    16361659        ENDIF 
Note: See TracChangeset for help on using the changeset viewer.