@@ -460,9 +460,9 @@ contains
460460 !! @param fCson Speed of sound (EL)
461461 !! @param adap_dt_stop Fail- safe exit if max iteration count reached
462462 pure subroutine s_advance_step (fRho , fP , fR , fV , fR0 , fpb , fpbdot , alf , &
463- fntait , fBtait , f_bub_adv_src , f_divu , &
464- bub_id , fmass_v , fmass_n , fbeta_c , &
465- fbeta_t , fCson , adap_dt_stop )
463+ fntait , fBtait , f_bub_adv_src , f_divu , &
464+ bub_id , fmass_v , fmass_n , fbeta_c , &
465+ fbeta_t , fCson , adap_dt_stop )
466466#ifdef _CRAYFTN
467467 !DIR$ INLINEALWAYS s_advance_step
468468#else
@@ -679,10 +679,10 @@ contains
679679 !! @param myPb_tmp Internal bubble pressure at each stage (EL)
680680 !! @param myMv_tmp Mass of vapor in the bubble at each stage (EL)
681681 pure subroutine s_advance_substep (err , fRho , fP , fR , fV , fR0 , fpb , fpbdot , alf , &
682- fntait , fBtait , f_bub_adv_src , f_divu , &
683- bub_id , fmass_v , fmass_n , fbeta_c , &
684- fbeta_t , fCson , h , &
685- myR_tmp , myV_tmp , myPb_tmp , myMv_tmp )
682+ fntait , fBtait , f_bub_adv_src , f_divu , &
683+ bub_id , fmass_v , fmass_n , fbeta_c , &
684+ fbeta_t , fCson , h , &
685+ myR_tmp , myV_tmp , myPb_tmp , myMv_tmp )
686686#ifdef _CRAYFTN
687687 !DIR$ INLINEALWAYS s_advance_substep
688688#else
@@ -708,7 +708,7 @@ contains
708708 myPb_tmp(1 ) = fpb
709709 myMv_tmp(1 ) = fmass_v
710710 call s_advance_EL(myR_tmp(1 ), myV_tmp(1 ), myPb_tmp(1 ), myMv_tmp(1 ), bub_id, &
711- fmass_n, fbeta_c, fbeta_t, mydPbdt_tmp(1 ), mydMvdt_tmp(1 ))
711+ fmass_n, fbeta_c, fbeta_t, mydPbdt_tmp(1 ), mydMvdt_tmp(1 ))
712712 end if
713713 myA_tmp(1 ) = f_rddot(fRho, fP, myR_tmp(1 ), myV_tmp(1 ), fR0, &
714714 myPb_tmp(1 ), mydPbdt_tmp(1 ), alf, fntait, fBtait, &
@@ -722,7 +722,7 @@ contains
722722 myPb_tmp(2 ) = myPb_tmp(1 ) + h* mydPbdt_tmp(1 )
723723 myMv_tmp(2 ) = myMv_tmp(1 ) + h* mydMvdt_tmp(1 )
724724 call s_advance_EL(myR_tmp(2 ), myV_tmp(2 ), myPb_tmp(2 ), myMv_tmp(2 ), &
725- bub_id, fmass_n, fbeta_c, fbeta_t, mydPbdt_tmp(2 ), mydMvdt_tmp(2 ))
725+ bub_id, fmass_n, fbeta_c, fbeta_t, mydPbdt_tmp(2 ), mydMvdt_tmp(2 ))
726726 end if
727727 myA_tmp(2 ) = f_rddot(fRho, fP, myR_tmp(2 ), myV_tmp(2 ), fR0, &
728728 myPb_tmp(2 ), mydPbdt_tmp(2 ), alf, fntait, fBtait, &
@@ -736,7 +736,7 @@ contains
736736 myPb_tmp(3 ) = myPb_tmp(1 ) + (h/ 4._wp )* (mydPbdt_tmp(1 ) + mydPbdt_tmp(2 ))
737737 myMv_tmp(3 ) = myMv_tmp(1 ) + (h/ 4._wp )* (mydMvdt_tmp(1 ) + mydMvdt_tmp(2 ))
738738 call s_advance_EL(myR_tmp(3 ), myV_tmp(3 ), myPb_tmp(3 ), myMv_tmp(3 ), &
739- bub_id, fmass_n, fbeta_c, fbeta_t, mydPbdt_tmp(3 ), mydMvdt_tmp(3 ))
739+ bub_id, fmass_n, fbeta_c, fbeta_t, mydPbdt_tmp(3 ), mydMvdt_tmp(3 ))
740740 end if
741741 myA_tmp(3 ) = f_rddot(fRho, fP, myR_tmp(3 ), myV_tmp(3 ), fR0, &
742742 myPb_tmp(3 ), mydPbdt_tmp(3 ), alf, fntait, fBtait, &
@@ -750,7 +750,7 @@ contains
750750 myPb_tmp(4 ) = myPb_tmp(1 ) + (h/ 6._wp )* (mydPbdt_tmp(1 ) + mydPbdt_tmp(2 ) + 4._wp * mydPbdt_tmp(3 ))
751751 myMv_tmp(4 ) = myMv_tmp(1 ) + (h/ 6._wp )* (mydMvdt_tmp(1 ) + mydMvdt_tmp(2 ) + 4._wp * mydMvdt_tmp(3 ))
752752 call s_advance_EL(myR_tmp(4 ), myV_tmp(4 ), myPb_tmp(4 ), myMv_tmp(4 ), &
753- bub_id, fmass_n, fbeta_c, fbeta_t, mydPbdt_tmp(4 ), mydMvdt_tmp(4 ))
753+ bub_id, fmass_n, fbeta_c, fbeta_t, mydPbdt_tmp(4 ), mydMvdt_tmp(4 ))
754754 end if
755755 myA_tmp(4 ) = f_rddot(fRho, fP, myR_tmp(4 ), myV_tmp(4 ), fR0, &
756756 myPb_tmp(4 ), mydPbdt_tmp(4 ), alf, fntait, fBtait, &
@@ -783,7 +783,7 @@ contains
783783 !! @param fdPbdt_tmp Rate of change of the internal bubble pressure
784784 !! @param fdMvdt_tmp Rate of change of the mass of vapor in the bubble
785785 pure elemental subroutine s_advance_EL(fR_tmp, fV_tmp, fPb_tmp, fMv_tmp, bub_id, &
786- fmass_n, fbeta_c, fbeta_t, fdPbdt_tmp, advance_EL)
786+ fmass_n, fbeta_c, fbeta_t, fdPbdt_tmp, advance_EL)
787787 !$acc routine seq
788788 real (wp), intent (IN ) :: fR_tmp, fV_tmp, fPb_tmp, fMv_tmp
789789 real (wp), intent (IN ) :: fmass_n, fbeta_c, fbeta_t
0 commit comments