Changeset 2793 for LMDZ5/trunk/libf/phylmd/iophys.F90

Timestamp:

Feb 16, 2017, 8:57:30 AM (8 years ago)

Author:

fhourdin

Message:

Petite correction a la precedente commission de iophys

File:

• LMDZ5/trunk/libf/phylmd/iophys.F90 (modified) (5 diffs)

LMDZ5/trunk/libf/phylmd/iophys.F90

@@ +1 @@
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+! Interface pour ecrire en netcdf avec les routines d'enseignement
+! iotd de Frederic Hourdin
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+
       subroutine iophys_ecrit(nom,lllm,titre,unite,px)
 
@@ -62 +67 @@
       end
 
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+! Version avec reindexation pour appeler depuis les routines internes
+! à la sous surface
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+
+
+    subroutine iophys_ecrit_index(nom,lllm,titre,unite,knon,knindex,px)
+
+    USE mod_phys_lmdz_para, ONLY: klon_omp
+    USE dimphy, ONLY : klon
+    USE mod_grid_phy_lmdz, ONLY: klon_glo
+    IMPLICIT NONE
+
+! This subroutine returns the sea surface temperature already read from limit.nc
+!
+
+! Arguments on input:
+    INTEGER lllm
+    CHARACTER (len=*) :: nom,titre,unite
+    REAL px(klon_omp,lllm)
+    INTEGER, INTENT(IN)                  :: knon     ! nomber of points on compressed grid
+    INTEGER, DIMENSION(klon), INTENT(IN) :: knindex  ! grid point number for compressed grid
+    REAL, DIMENSION(klon,lllm) :: varout
+
+    INTEGER :: i,l
+
+    IF (klon/=klon_omp) THEN
+      print*,'klon, klon_omp',klon,klon_omp
+      STOP'probleme de dimension parallele'
+    ENDIF
+
+    varout(1:klon,1:lllm)=0.
+    DO l = 1, lllm
+    DO i = 1, knon
+       varout(knindex(i),l) = px(i,l)
+    END DO
+    END DO
+    CALL iophys_ecrit(nom,lllm,titre,unite,varout)
+
+  END SUBROUTINE iophys_ecrit_index
+
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
       SUBROUTINE iophys_ini
@@ -68 +114 @@
       USE regular_lonlat_mod, ONLY: lon_reg, lat_reg
       USE dimphy, ONLY: klev
+      USE mod_grid_phy_lmdz, ONLY: klon_glo
 
       IMPLICIT NONE
@@ -88 +135 @@
 
 real pi
+INTEGER nlat_eff
 
 !   Arguments:
@@ -94 +142 @@
 !$OMP MASTER
     IF (is_mpi_root) THEN       
+
+! Bidouille pour gerer le fait que lat_reg contient deux latitudes
+! en version uni-dimensionnelle (chose qui pourrait être résolue
+! par ailleurs)
+IF (klon_glo==1) THEN
+   nlat_eff=1
+ELSE
+   nlat_eff=size(lat_reg)
+ENDIF
 pi=2.*asin(1.)
 call iotd_ini('phys.nc   ', &
-!iim,jjp1,llm,rlonv(1:iim)*180./pi,rlatu*180./pi,presnivs)
-size(lon_reg),size(lat_reg),klev,lon_reg(:)*180./pi,lat_reg*180./pi,presnivs)
+size(lon_reg),nlat_eff,klev,lon_reg(:)*180./pi,lat_reg*180./pi,presnivs)
     ENDIF
 !$OMP END MASTER