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

Timestamp:

Jul 22, 2024, 9:29:09 PM (2 months ago)

Author:

abarral

Message:

Replace most uses of CPP_DUST by the corresponding logical defined in lmdz_cppkeys_wrapper.F90
Convert several files from .F to .f90 to allow Dust to compile w/o rrtm/ecrad
Create lmdz_yoerad.f90
(lint) Remove "!" on otherwise empty line

File:

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

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

@@ -466 +466 @@
 !  emdustsco(k)=auxr3*tuningfactor
 !enddo
-!
-
-
-
 
 !JEdbg<<
@@ -904 +900 @@
 if (1==1) then
 !  !  CALL writefield_phy("AA",tmp1(1:klon,1:5),5)
-!
+
     CALL writefield_phy("REPART5",feff(1:klon,1:5),5)
     CALL writefield_phy("REPART5dbg",feffdbg(1:klon,1:5),5)
@@ -925 +921 @@
 !  if (.false.) then
 !!**************L718
-!
+
 !!c------------------------------------------------------------------------
 !! isolog distrib and massfrac calculations.
-!
+
 !      nbinsout=nbins+1
 !      b1=log(sizedustmin)
 !      b2=log(sizedustmax)
 !! restricted ISOLOG bins distributions
-!
+
 !!      step=(b2-b1)/(nbinsout-1)
 !!      DO ni=1,nbinsout
@@ -945 +941 @@
 !     binsHR(nb)=exp(b1+(nb-1)*stepbin)
 !  END DO
-!
+
 !  DO nb=1,nbinsHR
 !     binsHRcm(nb)=binsHR(nb)*1.e-4
@@ -951 +947 @@
 !! Making HIGH RESOLUTION dry deposition velocity
 !     CALL calcvd(vdout)
-!
-!
+
+
 ! DO nb=1,nbinsHR
 !     vdHR(nb)=vdout(nb)
 !!  WRITE(18,*),binsHR(nb),vdHR(nb)
 !  END DO
-!
+
 !   !searching for minimum value of dry deposition velocity
 !  minisograd=1.e20
@@ -966 +962 @@
 !     END IF
 !  END DO
-!
+
 !! searching for optimal number of bins in positive slope Vd part
-!
+
 !  nbins1=1
 !  nbins2=nbinsout-1
@@ -976 +972 @@
 !        IF(delta2.GE.delta1)THEN
 !        GOTO 50
-!
+
 !        ELSE
 !           nbins2=nbins2-1
@@ -994 +990 @@
 !     logvdISOGRAD(k)=logvdISOGRAD(1)-(k-1)*delta1
 !  END DO
-!
+
 !  logvdISOGRAD(nbins1+1)=log(minisograd)
-!
+
 !  DO k=1,nbins2
 !     logvdISOGRAD(nbins1+1+k)=logvdISOGRAD(nbins1+1)+k*delta2
 !  END DO
-!
+
 !  DO k=1,nbinsout
 !     vdISOGRAD(k)=exp(logvdISOGRAD(k))
@@ -1040 +1036 @@
 
 ! Making dust size distribution (in um)
-!
+
       nbinsout=nbins+1
       b1=log(sizedustmin)
@@ -1053 +1049 @@
 !     binsHR(nb)=exp(b1+(nb-1)*stepbin)
 !  END DO
-!
+
 !  DO nb=1,nbinsHR
 !     binsHRcm(nb)=binsHR(nb)*1.e-4
@@ -1188 +1184 @@
 !  print
 !*,'zwstar=sqrt(2.*(',flag_wstarBL,'ale_bl+0.01*(',flag_wstar,'-100)*ale_wake))'
-  !
+
     DO i=1,klon  ! main loop
      zwstar(i)=sqrt(2.*(param_wstarBL(i)*ale_bl(i)+param_wstarWAKE(i)*ale_wake(i)))