Changeset 2667 for LMDZ5/trunk/libf/phylmd/limit_read_mod.F90

Timestamp:

Oct 13, 2016, 11:41:08 AM (8 years ago)

Author:

fhourdin

Message:

Suite de la correction precedente

File:

• LMDZ5/trunk/libf/phylmd/limit_read_mod.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/limit_read_mod.F90

@@ -237 +237 @@
           ierr=NF90_INQUIRE(nid, UnlimitedDimID=ndimid)
           ierr=NF90_INQUIRE_DIMENSION(nid, ndimid, len=nn)
-          str=""; DO WHILE(nn>0); str=TRIM(str)//CHAR(nn-10*(nn/10)-48); nn=nn/10; END DO
+          CALL num2str(nn,str)
           abort_message='limit.nc records number ('//TRIM(str)//') does'//&
             ' not match year length ('
-          nn=year_len
-          str=""; DO WHILE(nn>0); str=TRIM(str)//CHAR(nn-10*(nn/10)-48); nn=nn/10; END DO
+          CALL num2str(year_len,str)
           abort_message=TRIM(abort_message)//TRIM(str)//')'
           IF(nn/=year_len) CALL abort_physic(modname,abort_message,1)
@@ -248 +247 @@
           ierr=NF90_INQ_DIMID(nid, 'points_physiques', ndimid)
           ierr=NF90_INQUIRE_DIMENSION(nid, ndimid, len=nn)
-          str=""; DO WHILE(nn>0); str=TRIM(str)//CHAR(nn-10*(nn/10)-48); nn=nn/10; END DO
+          CALL num2str(nn,str)
           abort_message='limit.nc horizontal number of cells ('//TRIM(str)//') does'//&
             ' not match LMDZ klon_glo ('
-          nn=klon_glo
-          str=""; DO WHILE(nn>0); str=TRIM(str)//CHAR(nn-10*(nn/10)-48); nn=nn/10; END DO
+          CALL num2str(klon_glo,str)
           abort_message=TRIM(abort_message)//TRIM(str)//')'
           IF(nn/=klon_glo) CALL abort_physic(modname,abort_message,1)
@@ -372 +370 @@
   END SUBROUTINE limit_read_tot
 
+  !--------------------------------------------------------------------------------------
+  SUBROUTINE num2str(n,str)
+  !--------------------------------------------------------------------------------------
+  ! Arguments:
+    INTEGER,           INTENT(IN)  :: n
+    CHARACTER(LEN=99), INTENT(OUT) :: str
+  !--------------------------------------------------------------------------------------
+  ! Local variables:
+    INTEGER :: nn
+  !--------------------------------------------------------------------------------------
+    nn=n; str=''; DO WHILE(nn>0); str=CHAR(nn-10*(nn/10)-48)//TRIM(str); nn=nn/10; END DO
+  END SUBROUTINE num2str
+  !--------------------------------------------------------------------------------------
 
 END MODULE limit_read_mod