Changeset 5103 for LMDZ6/branches/Amaury_dev/libf/phylmd/iophy.F90

Timestamp:

Jul 23, 2024, 3:29:36 PM (8 weeks ago)

Author:

abarral

Message:

Handle CPP_INLANDSIS in lmdz_cppkeys_wrapper.F90
Remove obsolete key wrgrads_thermcell, _ADV_HALO, _ADV_HALLO, isminmax
Remove redundant uses of CPPKEY_INCA (thanks acozic)
Remove obsolete misc/write_field.F90
Remove unused ioipsl_* wrappers
Remove calls to WriteField_u with wrong signature
Convert .F -> .[fF]90
(lint) uppercase fortran operators
[note: 1d and iso still broken - working on it]

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/iophy.F90 (modified) (25 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/iophy.F90

@@ -12 +12 @@
   INTEGER, ALLOCATABLE, DIMENSION(:), SAVE :: nptabij
   INTEGER, SAVE :: itau_iophy
-  LOGICAL :: check_dim = .false.
+  LOGICAL :: check_dim = .FALSE.
 
 !$OMP THREADPRIVATE(itau_iophy)
@@ -43 +43 @@
   USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat, klon_glo, grid_type, unstructured
   USE print_control_mod, ONLY: prt_level,lunout 
-#ifdef CPP_IOIPSL
     USE ioipsl, ONLY: flio_dom_set
-#endif
   use wxios, ONLY: wxios_domain_param, wxios_domain_param_unstructured, wxios_context_init, using_xios
     IMPLICIT NONE
@@ -200 +198 @@
     
 #ifndef CPP_IOIPSL_NO_OUTPUT 
-    call flio_dom_set(mpi_size,mpi_rank,ddid,dsg,dsl,dpf,dpl,dhs,dhe, &
+    CALL flio_dom_set(mpi_size,mpi_rank,ddid,dsg,dsl,dpf,dpl,dhs,dhe, &
                       'APPLE',phys_domain_id)
 #endif
@@ -324 +322 @@
     DO i=1,pim
 
-!    print*,'CFMIP_iophy i tabij lon lat',i,tabij(i),plon(i),plat(i)
+!    PRINT*,'CFMIP_iophy i tabij lon lat',i,tabij(i),plon(i),plat(i)
 
      plon_bounds(i,1)=rlon_glo(tabij(i)-1)
@@ -346 +344 @@
      plat_bounds(i,2)=rlat_glo(tabij(i)+nbp_lon)
 
-!    print*,'CFMIP_iophy point i lon lon_bds',i,plon_bounds(i,1),rlon_glo(tabij(i)),plon_bounds(i,2) 
-!    print*,'CFMIP_iophy point i lat lat_bds',i,plat_bounds(i,1),rlat_glo(tabij(i)),plat_bounds(i,2) 
+!    PRINT*,'CFMIP_iophy point i lon lon_bds',i,plon_bounds(i,1),rlon_glo(tabij(i)),plon_bounds(i,2)
+!    PRINT*,'CFMIP_iophy point i lat lat_bds',i,plat_bounds(i,1),rlat_glo(tabij(i)),plat_bounds(i,2)
 
     ENDDO
@@ -370 +368 @@
 
     DO i=1,pim
-!    print*,'CFMIP_iophy i tabij lon lat',i,tabij(i),plon(i),plat(i)
+!    PRINT*,'CFMIP_iophy i tabij lon lat',i,tabij(i),plon(i),plat(i)
 
      plon_bounds(i,1)=zx_lon(ipt(i)-1,jpt(i))
@@ -397 +395 @@
      ENDIF
 
-!    print*,'CFMIP_iophy point i lon lon_bds',i,plon_bounds(i,1),rlon(tabij(i)),plon_bounds(i,2) 
-!    print*,'CFMIP_iophy point i lat lat_bds',i,plat_bounds(i,1),rlat(tabij(i)),plat_bounds(i,2) 
+!    PRINT*,'CFMIP_iophy point i lon lon_bds',i,plon_bounds(i,1),rlon(tabij(i)),plon_bounds(i,2)
+!    PRINT*,'CFMIP_iophy point i lat lat_bds',i,plat_bounds(i,1),rlat(tabij(i)),plat_bounds(i,2)
 
     ENDDO
@@ -417 +415 @@
       ENDIF 
      ENDDO
-!    print*,'CFMIP_iophy mpi_rank npstn',mpi_rank,npstn
+!    PRINT*,'CFMIP_iophy mpi_rank npstn',mpi_rank,npstn
      IF(.NOT. ALLOCATED(nptabij)) THEN
       ALLOCATE(nptabij(npstn))
@@ -428 +426 @@
        npproc=npproc+1
        nptabij(npproc)=tabij(ip)
-!      print*,'mpi_rank npproc ip plon plat tabij=',mpi_rank,npproc,ip, &
+!      PRINT*,'mpi_rank npproc ip plon plat tabij=',mpi_rank,npproc,ip, &
 !      plon(ip),plat(ip),tabij(ip)
        npplon(npproc)=plon(ip)
@@ -439 +437 @@
 !!! print qui sert a reordonner les points stations selon l'ordre CFMIP
 !!! ne pas enlever
-        print*,'iophy_mpi rank ip lon lat',mpi_rank,ip,plon(ip),plat(ip)
+        PRINT*,'iophy_mpi rank ip lon lat',mpi_rank,ip,plon(ip),plat(ip)
 !!!
       ENDIF
@@ -638 +636 @@
             var%description, var%unit, var%flag(iff), typeecrit)
             IF (prt_level >= 10) THEN
-              WRITE(lunout,*) 'histdef2d: call wxios_add_field_to_file var%name iff: ', &
+              WRITE(lunout,*) 'histdef2d: CALL wxios_add_field_to_file var%name iff: ', &
                               trim(var%name),iff
             ENDIF
@@ -753 +751 @@
           var%description, var%unit, var%flag(iff), typeecrit)
             IF (prt_level >= 10) THEN
-              WRITE(lunout,*) 'histdef3d: call wxios_add_field_to_file var%name iff: ', &
+              WRITE(lunout,*) 'histdef3d: CALL wxios_add_field_to_file var%name iff: ', &
                               trim(var%name),iff
             ENDIF
@@ -851 +849 @@
      ELSE
       DO ip=1, npstn
-!     print*,'histwrite2d is_sequential npstn ip name nptabij',npstn,ip,name,nptabij(ip)
+!     PRINT*,'histwrite2d is_sequential npstn ip name nptabij',npstn,ip,name,nptabij(ip)
        IF(nptabij(ip)>=klon_mpi_begin.AND. &
           nptabij(ip)<=klon_mpi_end) THEN
@@ -975 +973 @@
   INTEGER :: ip
   REAL, ALLOCATABLE, DIMENSION(:) :: fieldok
-  logical, save :: is_active = .true.
+  logical, save :: is_active = .TRUE.
 
   IF (check_dim .AND. is_master) WRITE(lunout,*)'histwrite2d_phy for ',trim(var%name)
@@ -1019 +1017 @@
         !$omp barrier
         !$omp master
-        is_active = xios_field_is_active(var%name, at_current_timestep_arg=.false.)
+        is_active = xios_field_is_active(var%name, at_current_timestep_arg=.FALSE.)
         !$omp end master
         !$omp barrier
@@ -1044 +1042 @@
 
 ! La boucle sur les fichiers:
-      firstx=.true.
+      firstx=.TRUE.
 
       IF (ok_all_xml) THEN
@@ -1082 +1080 @@
                       write(lunout,*)'Dans iophy histwrite2D,iff,var%name ',&
                                      iff,trim(var%name)                       
-                      write(lunout,*)"histwrite2d_phy:.NOT.clef_stations(iff)and iff==iff_beg, call xios_send_field"
+                      write(lunout,*)"histwrite2d_phy:.NOT.clef_stations(iff)and iff==iff_beg, CALL xios_send_field"
                    ENDIF
                    IF (grid_type==regular_lonlat) THEN
@@ -1098 +1096 @@
                    ENDIF
 
-                   firstx=.false.
+                   firstx=.FALSE.
                 ENDIF
               ENDIF
@@ -1111 +1109 @@
 !                        IF (iff == iff_beg) THEN
 !                          IF (prt_level >= 10) THEN
-!                            write(lunout,*)"histwrite2d_phy: .NOT.clef_stations(iff) and iff==iff_beg, call xios_send_field"
+!                            write(lunout,*)"histwrite2d_phy: .NOT.clef_stations(iff) and iff==iff_beg, CALL xios_send_field"
 !                          ENDIF
 !                          CALL xios_send_field(var%name, Field2d)
@@ -1135 +1133 @@
 #ifndef CPP_IOIPSL_NO_OUTPUT 
                        IF (prt_level >= 10) THEN
-                         write(lunout,*)"histwrite2d_phy: clef_stations(iff) and iff==iff_beg, call wxios_write_2D"
+                         write(lunout,*)"histwrite2d_phy: clef_stations(iff) and iff==iff_beg, CALL wxios_write_2D"
                        ENDIF
                        CALL histwrite(nid_files(iff),var%name,itau_iophy,fieldok,npstn,index2d)
@@ -1185 +1183 @@
   INTEGER, ALLOCATABLE, DIMENSION(:) :: index3d
   REAL,ALLOCATABLE, DIMENSION(:,:) :: fieldok
-  logical, save :: is_active = .true.
+  logical, save :: is_active = .TRUE.
 
   IF (check_dim .AND. is_master) WRITE(lunout,*)'histwrite3d_phy for ', trim(var%name)
@@ -1217 +1215 @@
         !$omp barrier
         !$omp master
-        is_active = xios_field_is_active(var%name, at_current_timestep_arg=.false.)
+        is_active = xios_field_is_active(var%name, at_current_timestep_arg=.FALSE.)
         !$omp end master
         !$omp barrier
@@ -1247 +1245 @@
 
 ! BOUCLE SUR LES FICHIERS
-    firstx=.true.
+    firstx=.TRUE.
 
     IF (ok_all_xml) THEN
@@ -1281 +1279 @@
                   write(lunout,*)'Dans iophy, histwrite3D iff nlev klev firstx', &
                                   iff,nlev,klev, firstx                       
-                  write(lunout,*)'histwrite3d_phy: call xios_send_field for ', &
+                  write(lunout,*)'histwrite3d_phy: CALL xios_send_field for ', &
                                   trim(var%name), ' with iim jjm nlevx = ', &
                                   nbp_lon,jj_nb,nlevx
@@ -1299 +1297 @@
                 ENDIF
 
-                firstx=.false.
+                firstx=.FALSE.
               ENDIF
             ENDIF