Changeset 5160 for LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/dustemission_mod.F90

Timestamp:

Aug 3, 2024, 2:56:58 PM (7 weeks ago)

Author:

abarral

Message:

Put .h into modules

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/dustemission_mod.F90 (modified) (17 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/dustemission_mod.F90

@@ -350 +350 @@
    CALL abort_gcm('adaptdustemission', 'Dust range problem',1)
   ENDIF
-  print *,'FINE DUST BIN: tuning EMISSION factor= ',tuningfactorfine
-  print *,'COA DUST BIN: tuning EMISSION factor= ',tuningfactorcoa
-  print *,'SCO DUST BIN: tuning EMISSION factor= ',tuningfactorsco
-  print *,'ALL DUST BIN: SUM to the emissions (mg/m2/s) = ',basesumemission
+  PRINT *,'FINE DUST BIN: tuning EMISSION factor= ',tuningfactorfine
+  PRINT *,'COA DUST BIN: tuning EMISSION factor= ',tuningfactorcoa
+  PRINT *,'SCO DUST BIN: tuning EMISSION factor= ',tuningfactorsco
+  PRINT *,'ALL DUST BIN: SUM to the emissions (mg/m2/s) = ',basesumemission
   auxr1=9999.
   auxr2=9999.
@@ -418 +418 @@
 
 
-!print *,'JE'
-!print *,iminacclow,iminacchigh,imincoalow,imincoahigh
+!PRINT *,'JE'
+!PRINT *,iminacclow,iminacchigh,imincoalow,imincoahigh
 
 ! estimate and integrate bins into only accumulation and coarse
@@ -649 +649 @@
   varname='A'
   CALL read_surface(varname,Aini)
-print *,'beforewritephy',mpi_rank,omp_rank
+PRINT *,'beforewritephy',mpi_rank,omp_rank
   CALL writefield_phy("SOLinit",solini,5)
   CALL writefield_phy("Pinit",Pini,5)
@@ -657 +657 @@
   CALL writefield_phy("Dinit",Dini,5)
   CALL writefield_phy("Ainit",Aini,5)
-print *,'afterwritephy',mpi_rank,omp_rank
+PRINT *,'afterwritephy',mpi_rank,omp_rank
 
   DO i=1,klon
@@ -667 +667 @@
     ENDDO
   ENDDO
-!print *,'JEOK1',mpi_rank,omp_rank
+!PRINT *,'JEOK1',mpi_rank,omp_rank
   DO i=1,klon
-    !print *,Pini(i,1),Pini(i,2),Pini(i,3),Pini(i,4),Pini(i,5)
+    !PRINT *,Pini(i,1),Pini(i,2),Pini(i,3),Pini(i,4),Pini(i,5)
     DO nts=1,ntyp
       !IF(xlon(i).ge.longmin.AND.xlon(i).le.longmax.AND. &
@@ -698 +698 @@
   ENDDO
 
-! print *,'JEOK2',mpi_rank,omp_rank
+! PRINT *,'JEOK2',mpi_rank,omp_rank
 IF ( 1==1 ) THEN
-! print *,'JEOK4',mpi_rank,omp_rank
+! PRINT *,'JEOK4',mpi_rank,omp_rank
    CALL writefield_phy("SOL",sol,5)
    CALL writefield_phy("P"  ,P  ,5)
@@ -712 +712 @@
    CALL writefield_phy("xlon",xlon,1)
    CALL writefield_phy("xlat",xlat,1)
-!print *,'JEOK5',mpi_rank,omp_rank
-!print *,'JEOK6',mpi_rank,omp_rank
+!PRINT *,'JEOK5',mpi_rank,omp_rank
+!PRINT *,'JEOK6',mpi_rank,omp_rank
 
 END IF
@@ -724 +724 @@
 ! data format in columns
 ! mmd1  sigma1  p1  mmd2  sigma2  p2 mmd3 ... alpha
-!print *,'JEOK7',mpi_rank,omp_rank
+!PRINT *,'JEOK7',mpi_rank,omp_rank
 !$OMP MASTER
 IF (is_mpi_root .AND. is_omp_root) THEN
-!print *,'JEOK9',mpi_rank,omp_rank
+!PRINT *,'JEOK9',mpi_rank,omp_rank
  fnsolspe='SOILSPEC.data'
   PRINT*,'  o Reading ',fnsolspe(1:40)
@@ -744 +744 @@
 !$OMP END MASTER
 !$OMP BARRIER
-!print *,'JEOK10',mpi_rank,omp_rank
+!PRINT *,'JEOK10',mpi_rank,omp_rank
   CALL bcast(solspe)
 ! Calcul de la distribution en taille des particules de Dust
@@ -769 +769 @@
 !if (.TRUE.) THEN
 !c 0: Iversen and White 1982
-       print *,'Using  Iversen and White 1982 Uth'
+       PRINT *,'Using  Iversen and White 1982 Uth'
          DO i=1,ncl
             bb=adust*(sizeclass(i)**xdust)+bdust
@@ -786 +786 @@
 IF(.TRUE.) THEN
 ! 1: Shao and Lu 2000
-       print *,'Using  Shao and Lu 2000 Uth'
+       PRINT *,'Using  Shao and Lu 2000 Uth'
             an=0.0123
             gam=0.3
@@ -886 +886 @@
 ! JE20150120<<
   IF (flag_feff == 0) THEN
-    print *,'JE_dbg FORCED deactivated feff'
+    PRINT *,'JE_dbg FORCED deactivated feff'
     DO i=1,klon
       DO k=1,ntyp
@@ -1072 +1072 @@
             numax=log(d2max/diam(k))
             massfrac(k,nb)=0.5*(erf(numax/exden)-erf(numin/exden))
-            !print *,k,nb,massfrac(k,nb)
+            !PRINT *,k,nb,massfrac(k,nb)
          ENDDO
       ENDDO
@@ -1209 +1209 @@
                    pdfu=(kref/weilambda)*(wind10ms(kwb)/weilambda)**(kref-1) &
                       *exp(-(wind10ms(kwb)/weilambda)**kref)
-!                   !print *,'JEdbg  U10mMOD weilambda  ',U10mMOD,weilambda
+!                   !PRINT *,'JEdbg  U10mMOD weilambda  ',U10mMOD,weilambda
 !JE20141205>>
 
@@ -1228 +1228 @@
              IF(.FALSE.) THEN
 !                  nat=int(sol(i,n))
-!                    print *,i,n
+!                    PRINT *,i,n
                     IF(sol(i,n)>1..and.sol(i,n)<15.) nat=int(sol(i,n))
 !JE20140526<<
-!                    print *,'JE: WARNING: nat=0 forced to nat=99!! and doing nothing'
+!                    PRINT *,'JE: WARNING: nat=0 forced to nat=99!! and doing nothing'
                    IF(sol(i,n)<0.5) THEN
                       nat=99
@@ -1317 +1317 @@
         flux2=flux2+ft2*cpcent*cerod
         flux3=flux3+ft3*cpcent*cerod
-!       print *,'JEflux :',kwb,n,flux1,flux2,flux3
+!       PRINT *,'JEflux :',kwb,n,flux1,flux2,flux3
 80 CONTINUE
              ENDDO !n=1,ntyp
@@ -1348 +1348 @@
 
 
-!print *,' JE fluxdust in calcdust'
+!PRINT *,' JE fluxdust in calcdust'
 ! DO k=1,klon
 !   DO i=1,ndistb
-!!print *,k,i,fluxdust(k,i)
+!!PRINT *,k,i,fluxdust(k,i)
 !enddo
 !enddo
-!print *,' JE emisbin in calcdust'
+!PRINT *,' JE emisbin in calcdust'
 !do k=1,klon
 !do j=1,nbins
-!!print *,k,j,emisbin(k,j)
+!!PRINT *,k,j,emisbin(k,j)
 !enddo
 !enddo